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Abstract 

Transmission capacity (TC) is a performance metric for wireless net- 
works that measures the spatial intensity of successful transmissions 
per unit area, subject to a constraint on the permissible outage prob- 
ability (where outage occurs when the SINR at a receiver is below a 
threshold) . This volume gives a unified treatment of the TC framework 
that has been developed by the authors and their collaborators over the 
past decade. The mathematical framework underlying the analysis (re- 
viewed in Ch. 2) is stochastic geometry: Poisson point processes model 
the locations of interferers, and (stable) shot noise processes represent 
the aggregate interference seen at a receiver. Ch. 3 presents TC results 
(exact, asymptotic, and bounds) on a simple model in order to illustrate 
a key strength of the framework: analytical tractability yields explicit 
performance dependence upon key model parameters. Ch. 4 presents 
enhancements to this basic model — channel fading, variable link dis- 
tances, and multi-hop. Ch. 5 presents four network design case studies 
well-suited to TC: i) spectrum management, ii) interference cancella- 
tion, in) signal threshold transmission scheduling, and iv) power con- 
trol. Ch. 6 studies the TC when nodes have multiple antennas, which 
provides a contrast vs. classical results that ignore interference. 



Contents 



1 Introduction and preliminaries [l| 

1.1 Motivation and assumptions 2 

1.2 Key definitions: PPP, OP, and TC [s 

1.3 Overview of tlie results 



2 Mathematical preliminaries |15 

2.1 Probability: notations, definitions, key inequalities IJ 

2.2 PPP void probabilities and distance mappings 2C 

2.3 Shot noise (SN) processes 22 

2.4 Stable distributions, Laplace transforms, and PGFL 3C 

2.5 Maximums and sums of RVs |35 

3 Basic model |39 

3.1 Exact OP and TC 4C 

3.2 Asymptotic OP and TC 43 

3.3 Upper bound on TC and lower bound on OP M 

3.4 Throughput (TP) and TC m 



1 



ii Contents 



3.5 Lower bounds on TC and upper bounds on OP 

4 Extensions to the basic model 

4.1 Channel fading 

4.2 Variable link distances (VLD) 

4.3 Multihop TC 

5 Design techniques for wireless networks 

5.1 Spectrum management 

5.2 Interference cancellation (IC) 

5.3 Fading threshold scheduling (FTS) 

5.4 Fractional power control (FPC) 

6 Multiple antennas 

6.1 MIMO with interference 

6.2 Categorizing MIMO in decentralized networks 

6.3 Single stream MIMO TC results 

6.4 Main results on multiple stream TC 

6.5 Practical issues and further research 

Notations and Acronyms 
A List of results by chapter 
Acknowledgements 



References 



1 



Introduction and preliminaries 



Wireless networks are becoming ever more pervasive, and the cor- 
respondingly denser deployments make interference management and 
spatial reuse of spectrum defining aspects of wireless network design. 
Understanding the fundamentals of the performance and behavior of 
such networks is an important theoretical endeavor, but one with 
only limited success to date. Information theoretic approaches, well- 
summarized by [24J, have been most successfull when applied to small 
isolated networks, where background interference and spatial reuse are 
not considered. Large network approaches, typified by transport capac- 
ity scaling laws [82], have given considerable insight into scaling laws, 
but are generally unable to quantify the relative merits of candidate 
design choices or provide a tractable approach to analysis for spatial 
reuse or the SINR statistics. Our hope for the transmission capacity 
framework has been to develop a tractable approach to large network 
throughput analysis, that while falling short of information theory's 
ideals of inviolate upper bounds, nevertheless provides a rigorous and 
flexible approach to the same sort of questions, and ultimately pro- 
vides the types of broad design insights that information theory has 
been able to achieve for small networks. 
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2 Introduction and preliminaries 



1.1 Motivation and assumptions 

This monograph presents a framework for computing the outage prob- 
abihty (OP) and transmission capacity (TC) ^80^ i79j in a wireless net- 
work. The OP is defined as the probabiUty that a "typical" transmis- 
sion attempt fails (is in outage) at the intended receiver, where outage 
occurs when the signal to interference plus noise ratio (SINR) at the 
receiver is below a threshold. Basing outage on the SINR, it is assumed 
that interference is treated as noise. The TC is defined as the maxi- 
mum average number of concurrent successful transmissions per unit 
area taking place in the network, subject to a constraint on OP. The 
OP constraint may be thought of as a reliability and/or quality of ser- 
vice (QoS) parameter — strict requirements on the fraction of failed 
transmissions result in low spatial reuse, low area spectral efficiency 
(ASE, measured in bps/Hz per unit area), and thus lower TC, while 
relaxing the outage requirement improves, up to a point, the ASE and 
thus TC. Viewing the OP as a (strictly increasing) function of the in- 
tensity of attempted transmissions, the TC is computed by inverting 
this function for the transmission intensity at the target OP. 

Note we use the word capacity in a distinctly different manner from 
its information-theoretic sense, i.e., Shannon capacity: the TC frame- 
work typically treats interference as nois^ while Shannon theory does 
not, and TC measures capacity in a spatial sense, while Shannon the- 
ory does not. The capacity in TC is also distinct from the transport 
capacity of |34] , defined as the maximum weighted sum rate of commu- 
nication over all pairs of nodes, where each pair's communication rate 
is weighted by the distance separating them. The transport capacity 
optimizes over all scheduling and routing algorithms and the focus is 
on the asymptotic rate of growth of the sum rate in the number of 
nodes n, either keeping the network area fixed or letting the network 
area grow linearly with n. TC, on the other hand, is a medium access 
control (MAC) layer metric that neither precludes nor addresses rout- 
ing. Although transport capacity is more general in that it optimizes 



^ see i)5.2l and the results of Chapter 6 as an exception: even here though the background 

(uncancelled) interference is then treated as noise. 
^ see i)4.3l for an exception, where a simple multihop model is added. 
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scheduling and routing, the cost of this generaUty is that typically the 
transport capacity results are less specific than those obtainable under 
the TC framework. The results are less specific in the sense that re- 
sults on the asymptotic rate of growth of the transport capacity as a 
function of n often do not specify the pre-constant. 

The advantages of using TC as a metric for wireless network per- 
formance are: i) it can be exactly derived in some important cases, and 
tightly bounded in many others, ii) performance dependencies upon 
fundamental network parameters are thereby illuminated, and in) de- 
sign insights are obtainable from these performance expressions. More 
fundamentally, the TC captures in a natural way essential performance 
indicators like network efficiency (ASE), reliability (OP), and through- 
put (TP). In fact, TC is precisely maximization of TP under an OP 
constraint, as discussed in §3.41 and is proportional to ASE, as discussed 
in gSH 

One limitation of the TC framework, at least as in this monograph, 
is the implicit assumption that the network employs the simplistic and 
sub-optimal slotted Aloha protocol at the MAC layer. The TC can 
also be extended to model other contention based MAC protocols at 
the cost of some tractability \26\ |25J , but we elect to stick to the sim- 
ple slotted Aloha protocol, where each transmitter (Tx) independently 
elects whether or not to transmit to its receiver (Rx) in each time slot 
by flipping an independent biased coin [1]. If the point process describ- 
ing the locations of contending transmitters at a snapshot in time form 
a Poisson point process (PPP), which we assume is the case, then un- 
der the Aloha protocol the locations of the active transmitters at some 
point in time also form a PPP, obtained by independent sampling of 
the node location PPP. The PPP model is necessary for preserving 
the highest level of analytical tractability of the TC framework, but of 
course it means that the computed TC is sub-optimal. The difficulty 
in relaxing the Aloha assumption lies in the fact that any realistic and 
useful randomized MAC protocol involves coordination among com- 
peting transmitters, which necessarily spoils the crucial independence 
property of the PPP. One's valuation of the TC framework typically 
rests on weighing the advantage of having an explicit expression for an 
insightful network performance indicator with the disadvantage of that 
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performance corresponding to a suboptimal control law. Throughout 
this monograph, we generally adopt the following assumptions in order 
for the OP and TC to be computable. See Fig. 11.11 



Assumption 1.1. The following assumptions are made: 

(1) The network is viewed at a single snapshot in time for the 
purpose of characterizing its spatial statistics. 

(2) Every potential Tx is matched with a prearranged intended 
Rx at a fixed distance u (meters) awajo: these Tx-Rx pairs 
are in one to one correspondence. 

(3) When mapping our results to a specific bit rate, we assume 
each Rx treats (uncancelled) interference as noise, and the 
rate from a particular Tx to its Rx at location o is given by 
the Shannon capacity c(o) = |log2(l +sinr(o)). 

(4) The potential transmitters form a homogeneous PPP on the 
network arena, taken to be M'^, for d £ {1, 2, 3}. This implies 
i) the number of nodes in the network is countably infinite, 
and a) the number of potential transmitters in two disjoint 
bounded sets of the plane are independent Poisson random 
variables (RV). See Fig. O 

(5) Every potential Tx decides independently whether or not to 
transmit with a common probability ptx- It follows that the 
set of actual transmitters is also a (thinned) PPP. 



A few remarks are in order: 

(1) The TC computes the maximum spatial reuse which is com- 
putable by looking at the network at a single snapshot in 
time. This perspective neither addresses nor precludes multi- 
hop or routing considerations. 

(2) Ass. (3) can be easily softened to account for any modulation 
and coding type that is characterized by a SINR "gap" from 
capacity. Typically, we directly utilize SINR for computing 



^ The extension to random distances is straightforward and given in 114.21 
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outage probability and TC and do not include the per-link 
rate in the results. 
(3) Ass. (4) makes clear our focus is on networks whose arena 

is the entire plane M?, and which have a countably infinite 
number of nodes. This along with Ass. (5) removes any con- 
cern about boundary effects and makes each node "typical" 
in a sense described below. 
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Fig. 1.1 A reference (typical) Rx located at the origin o G K'' is paired with a reference Tx 
at distance u, and is subject to interference (dashed lines) from a PPP \ of interfering 
Tx's (ea<;h of which has a unique associated Rx). 

1.2 Key definitions: PPP, OP, and TC 

Ass. (l)-(3) allow us to formally define the OP. 

Definition 1.1. Outage probability (OP). Define the constant R 

to be the spectral efficiency (in bits per channel use per Hz) of the 
channel code employed by each Tx-Rx pair in the network. Define 
the SINK threshold r = 22-^ - 1 so that R = |log2(l + r). For an 
arbitrary Tx-Rx pair with the Rx positioned at the origin o G M*^, let 
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c(o) = 2 log2(l + sinr(o)) be the random Shannon spectral efficiency of 
the channel connecting them when interference is treated as noise. The 
OP is the probability the random spectral efficiency of the channel falls 
below the spectral efficiency of the code, or equivalently, the probability 
that the random SINR at the Rx is below the threshold r: 

q{o) = P(c(o) < R) 

= pQlog2(l + sinr(o)) < ^log2(l + r)^ 

= P(sinr(o) < r). (1.1) 



It is worth emphasizing that the RV in P(c(o) < R) is the capacity c(o) 
of the channel connecting the Tx-Rx pair, computed at the snapshot 
in time at which we observe the network, and not the rate R, which is 
assumed fixed. In particular, c(o) is a function of the RV sinr(o), which 
is quite sensitive to the distances between the Rx at o and the random 
set of interfering transmitters at the observation instant. The OP is the 
cumulative distribution function (CDF) of the RVs c(o) and sinr(o). 

By the assumption that the transmitters (and receivers) form a 
PPP, it follows that all Tx-Rx pairs are typical, hence q{o) = q. More 
formally, we can condition on the presence of a test Tx-Rx pair where, 
without loss of generality, we assume the test Rx to be located at 
the origin o. The distribution of the PPP of potential transmitters is 
unaffected by the addition of this test pairQ 

q{o) = P(c(o) < R\ Rx at o) = P(c(o) < R). (1.2) 

Ass. (4)-(6) and the definition of OP allow us to formally define the TC. 
We first define a homogeneous PPP on MJ^. Fig. 11.21 shows a portion of 
a sample PPP on and illustrates the fact that the number of points 
in each compact set is a Poisson RV. 

Definition 1.2. Homogeneous Poisson point process (PPP). A 

PPP with intensity A > in d-dimensions is a random countable col- 
lection of points Ild,x = {xi)X25 . . .} C M'^ such that 



■*This result is due to Slivnyak [67]. See, e.g., g] Thm. 1.13 (p. 30), [36] Thm. A.5 (p.ll3), 
[To] p.41 and Example 4.3 (p.l21). 
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• For two disjoint subsets A, B C the number of points 
from n in these sets are independent RVs: n(A) _L 11(B), 
where n(A) = |n n A| is the number of points in 11 in A. 

• The number of points in any compact set ^ C M"^ is a Poisson 
RV with parameter X\A\, for |^| the volume of A. That is: 

n(yl)~Po(A|A|), (1.3) 

or equivalently, 

p(n(^) = k) = ^e-^l^l(A|^|)^ k e Z+. (1.4) 



Suppose Ilpot = ^d,Xpot is the PPP of potential transmitters with 
spatial intensity Apot discussed in Ass. 11.11 (4) and (5). Let ptx G (0, 1) 
be the common transmission probability employed by each node in Ass. 
11.11 (5). It follows that the PPP of actual transmitters at the observation 
instant (denoted 11^^^) is a thinned version of Ilpot with corresponding 
thinned intensity A = ApotPtx- 

It is intuitive that the OP is increasing in A: a higher spatial in- 
tensity of transmission attempts yields larger interference at each Rx, 
which decreases the SINR. We emphasize this dependence by writing 
the OP as q{X) and thereby view q : M+ — )■ [0, 1] as a map from the 
spatial intensity of transmission attempts to the corresponding OPifl 



Fact 1.1. The OP g(A) is continuous, strictly increasing, and onto 
[g(0), 1), where q{0) is the OP in the absence of interference. 



Because of this fact, the inverse q~^ : [q{0),l) — )• M+ is well-defined. 
For an outage constraint q* E [(/(O), 1), the inverse OP q~^{q*) is the 
(unique) intensity of transmission attempts associated with an outage 
probability of q* . Each such transmission succeeds with probability 
I — q*, and as such q~^{q*){l — q*) is the spatial intensity of successful 
transmissions. This is what we call TC. 



^ This redefines the OP (?(A) as a function of the intensity A £ R_|_ of the PPP — in HI. Il l 
and ljl.2|l q{x) denoted the OP at location a; S R**. 
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Fig. 1.2 Top: a realization of a homogeneous Poisson point process (PPP) of intensity 
A = 1 on R^. Bottom: the histogram of the number of points in each of the 10 X 10 unit 
area squares and the corresponding Poisson PMF for A = 1. 
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Definition 1.3. Transmission capacity (TC). Fix a maximum per- 
missible OP q* G [9(0), 1). The TC is the maximum spatial intensity of 
successful transmissions subject to an OP of q*: 

X(q*)^q-\q*){l-q*). (1.5) 



The intensity of failed transmissions is q^^{q*)q*, and the summed 
intensity of successful and failed transmissions is naturally q~^{q*). The 
TC (and all spatial intensities) are measured in units of (meters"''), i.e., 
an "average" number of nodes per unit area. 

Remark 1.1. TC and slotted Aloha. The TC X{q*) has operational 
sigificance for a wireless network of potential transmitters positioned 
according to a PPP of intensity Apot and employing the slotted Aloha 
MAC protocol with transmission probability ptx- Namely, if (g*, Apot) 
are such that X{q*) < Apot x (1 — <?*) then select 

Ptx= , (1.6) 
Apot X (l-q*) 

The resulting intensity of attempted transmissions ptxApot = q~^iq*) 
will be such that the OP is l-q*. li X{q*) > Apot x (1 - Q'*) then the 
network does not need an Aloha MAC throttling transmission attempts 
to achieve an OP of q*: setting ptx = 1 will result in an OP (/(Apot) < Q* ■ 



1.3 Overview of the results 

The results presented in this volume are listed in Tables lA.ll (Ch. [T]) 
through lA. 61 (Ch. [6|). We briefly discuss each chapter. 

Ch. [1] (Table IA.1[> . The key concepts are in §1.2( specifically, 
Def. Il.ll of the outage probability (OP), Def. 11.21 of the (homogeneous) 
Poisson point process (PPP), and Def. 11.31 of the transmission capacity 
(TC). 

Ch. [2] (Table [AT2]) . We first define the bah and annulus in M'^. 
(Def. 12. ip and gives their volumes (Prop. [27T]) . All results are given for 
arbitrary dimension d, where {1,2,3} are the three relevant values. 
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Throughout the volume we denote RVs in sans-serif font (Rem. 12.11 
in §2.ip . e.g., x. Note the acronyms and notation for standard proba- 
bihstic concepts in Def . 12.31 ^2.21 gives a short but essential coverage of 
the void probability (Prop. [T6|) . the mapping theorem (Thm. [27T]) and 
a derivative result on mapping distances (Prop. 12. 7p . The void proba- 
bility underlies most performance bounds derived in this volume, and 
the mapping result allow translation from a PPP on M'^ of intensity 
A G R+ (XI^^a) to an "equivalent" unit intensity PPP on M} (Ui^i). 

We cover (spatial) shot noise (SN) processes in §2.3l (Def. l2.4p . which 
are used to model the aggregate interference experienced by a reference 
Rx at the origin. We focus on power law SN (Def. 12. 5p by assuming 
the impulse response function in the SN definition is taken to be the 
standard pathloss attenuation |x|~" where a is the pathloss exponent. 
We also introduce here the characteristic exponent 6 = d/a, where to 
avoid trivialities we assume 5 £ (0, 1) {i.e., a > d) throughout. The sum 
SN process (Z) adds the interference contributions while the max SN 
(M) takes the largest contribution. The simple inequality M < Z forms 
the basis of most of the bounds in this volume in that the distribution 
of M (Cor. 12. ip is Frechet (Def. 12. 6p . and also can be derived directly 
from the void probability in Prop. 12.61 The Campbell-Mecke result 
(Thm. [23]) allows computation of moments (Prop. [2^91) of SN RVs. More 
important for us will be the series expansions of the SN distribution 
(Prop. I2.10p as these directly yield the asymptotic (tail) distributions 
(Cor. 12. 2p . which yield all the asymptotic performance results in this 
volume. 

A critical observation is that the SN is a stable RV, this is the focus 
of §2.41 We define this class (Def. 12.71 and 12. 8p , and introduce the Levy 
distribution (Def. 12. 9p which is the only stable distribution of relevance 
to us with a closed form CDF, and corresponding to S = ^. This allows 
the exact performance results in Ch. [3j We introduce the probability 
generating functional (PGFL) (Def. I2.10p . and identify its connection 
with the Laplace transform, the moment generating function, and the 
characteristic function of the SN RV. 

The results in this chapter are tied together in ^2.51 where we de- 
mostrate the key property that the the simple bound M < H is tight in 
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the sense that the ratio of the CCDFs for these RVs approaches unity in 
the hmit (Prop. I2.13p . We derive a similar result using subexponential 
distributions (Def. [TTT]) for a binomial point process (BPP). 

Ch. M (Table [X3]) . This chapter presents the main results on 
OP and TC in their barest, simplest form, so as to achieve maximum 
clarity. Exact OP and TC results are in gSH SINK is defined (Def. [3T]) 
and it is observed that the OP is the CCDF of the SN evaluated at a 
certain value. An explicit expression for the OP and TC (for 6 = ^) is 
given (Cor. O and E^D • 

Asymptotic OP and TC results are in §3.21 The asymptotic CCDF 
of the SN (Cor. [221) yields the asymptotic OP (as X ^ 0) and TC 
(as q* —7- 0) in Prop. 13. 3[ The asymptotic TC is interpreted as sphere 
packing in R"', where the sphere radius depends upon the key model 
parameters 6, u, r, d (Rem. 13. 1|) . 

The M < 51 SN inequality forms the basis for the OP lower bound 
(LB) and TC upper bound (UB) in §3.31 We adopt the language of 
dominant interferers (Def. 13. 3p to describe interferers capable by them- 
selves of reducing the SINR seen at the origin below its threshold r, 
but observe this concept is equivalent to taking the maximum interferer 
(Rem. 13. 2p . The main result is the bound on OP and TC in Prop. 13.41 

In ^3.41 we turn our attention to a third performance metric, the 
MAC layer throughput (TP), A(A), defined (Def. [33D as the spatial 
intensity of succesful transmissions. A TP UB is obtained from the 
OP LB (Prop. 13. 6p . We make the key observation that "blind" maxi- 
mization of TP leads to an associated OP of 67%. The natural design 
objective of maximizing TP subject to an OP constraint is shown to be 
precisely the TC, giving a more natural justification for this quantity as 
a meaningful performance measure (Prop. [3T7I) . In fact the TP and TC 
have the same unconstrained maximum and we relate their maximizers 
(Prop.EHD. 

Finally, ^3.51 gives an UB on OP and a LB on TC. A useful expres- 
sion for the OP in terms of its LB is derived (Prop. [3T9l) . which the OP 
LB and the three basic inequalities in §2.11 (Markov, Chebychev, and 
Chernoff) are combined to give three OP UBs. These are observed to 
vary both in terms of their tightness and their simplicity. 

Ch. [4] (Table IA.4P extends the basic model in three ways: fading 
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( P?T]l . variable link distances (P^. and multi-hop (P^. 

The bulk of this chapter is on fading ( §4.ip : the SINR under fading 
is defined (Def. 14. ^4.11 is split into three subsections: exact results 
( ^4.1.ip . asymptotic results ( ^4.1.2p . and bounds ( ^4.1.3p . The main 
result in §4.1.11 is Prop. 14.21 which gives the exact OP and TC under 
the assumption that the signal fading is Rayleigh (exponential). Note 
this exact result holds for all 5, while the only exact result available 
under the basic model in Ch. [3]is for (5 = ^ (Cor. [3T] and 13. 2p . For the 
asymptotic results in §4.1.21 we introduce the formalism of the marked 
PPP (MPPP) and exploit the important marking theorem (Thm. 14. ip 
which allows us to extend the distance and interference mapping re- 
sults for PPPs from ^2.21 to the MPPP case. The series expansions of 
the interference under fading (Prop. 14. 5p is used to derive the asymp- 
totic OP and TC (Prop. An important observation is that fading 
in general degrades performance relative to the non-fading case (Cor. 
14. 3p . In ^4.1.3l the concept of dominant interferers used in Def. I3.3l is ex- 
tended to incorporate fading (Def. 14. 2p . but under fading the strongest 
interferer need not be the nearest interferer to the origin. The main 
result is the OP LB (Prop. H77|l . where we observe the LB is in fact the 
MGF of a certain function of the signal fading RV. 

^4.21 addresses variable link distances, i.e., the Tx-Rx distance is a 
RV. The SINR and OP for this model are defined in Def . US] and [Ol 
respectively, and we present asymptotic results (Prop. 14. 8p and exact 
results (Cor.Ol). 

§331 extends the TC framework to a multihop scenario where sources 
send packets to destinations M hops away over a total distance R. Mul- 
tihop TC is defined in Def. 14.51 Although some fairly strong assump- 
tions must be made to preserve tractability, plausible insights can be 
drawn about the optimum hop count (given in Prop. I4.13P and end-to- 
end TC in terms of all the network parameters. 

Ch. [5] (Table [A.5p . The chapter on design techniques studies four 
natural approaches to improve the performance of a wireless network: 
§5. II studies the performance when the spectrum is split into a number 
of channels, ^5.21 considers performance when receivers are equipped 
with interference cancellation capabilities, ^5.31 evaluates the perfor- 
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mance when nodes only transmit when their signal fade is above a 
specified threshold, and §5.41 considers power control. 

In ^5.11 the design objective is to optimize the number of bands 
to form from the available spectrum, where each Tx selects a band 
uniformly at random. The intuitive tradeoff is that more bands gives 
fewer interferers but this also means the bandwidth per band is smaller, 
and thus a higher SINR threshold is required to achieve a given data 
rate. We define the model in Def. 15.11 and 15. 2[ The spectral efficiency 
optimization problem is formalized in Prop. 15.11 and we characterize 
the solution in Prop. 15.31 and then specialize the result to both the 
high (Cor. 15. 2p and low (Prop. 15. 4p SNR regimes. 

In ^5.21 the usual limitations of interference cancellation (IC) are 
captured through the (k, K, Pmm) model (Def. 15. 4p where k is can- 
cellation effectiveness, K is the maximum number of cancellable nodes, 
and Pmin is the minimum received power. The SINR is defined in Def. 
15. 5t and the main result is the OP LB (Prop. [53]) . 

In ^5.31 the fading coefficient threshold (Def. 15. 6p used to throt- 
tle transmission attempts naturally trades off between the quality and 
quantity of transmission attempts, and the TP metric A (Def. 15. 7p il- 
lustrates this tradeoff. Asymptotic results are given in Prop. 15.71 and a 
LB on OP is given in Prop. 15.91 

In ^5.41 the notion of fractional power control (FPC) is introduced, 
where the power control exponent sweeps between fixed power and 
channel inversion (Def. ISTS]) . The asymptotic results (Prop. [5.10p yield 
the optimal exponent is 1/2 (Prop. 15. lip . The notion of dominant in- 
terferers is used once again (Def. 15. 9p to compute the OP LB (Prop. 
[532]). 

Ch. [6] (Table Fa. 6h . The final chapter introduces multiple antennas 
at both the Tx and Rx, resulting in some of the first analytical work 
on MIMO that properly accounts for background interference. The re- 
sults are broken into two main categories, which are defined along with 
basics of the models in §6.21 §6.31 considers the case where despite the 
multiple antennas, only a single data stream is sent, with the balance of 
the antennas being used for diversity and/or interference cancellation. 
^6.41 considers the more general multistream case, where transmitters 
send more than one simultaneous stream to either a single receiver 
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(spatial multiplexing) or to multiple users (space division multiple ac- 
cess). Finally, the practical implications and limitations of the results 
are discussed in ^6.51 

In ^Q.3\ the results are further categorized into diversity ( ^6.3.ip and 
interference cancellation ( §6.3.2p . For receive diversity, the OP of MRC 
is given in Prop. 16.11 and the corresponding TC in 16.21 The result is 
equivalent for MRT (transmit MRC) and the generalization to 
diversity beamforming is discussed in Rem. 16.41 In ^6.3.21 a TC lower 
bound is given on a suboptimal technique called partial ZF in Prop. 
16.61 and a TC upper bound for MMSE in Prop. \67l\ These respectively 
bound the TC of MMSE and we see linear scaling can be achieved with 
the number of antennas. 

In j ]6.4l we first consider a class of results for spatial multiplexing 
in §6.4. H where multiple streams are transmitted from a single Tx to a 
single Rx. Prop . ISTSl and [6^91 give the optimal number of streams K* and 
TC scaling in terms of tit < for MRC and ZF receivers, respectively. 
This is extended to a BLAST receiver in Prop. 16.121 Then in ^6.4.21 we 
turn our attention to streams being sent to multiple Rx's at the same 
time. The main result for MRC receivers is given in Prop. 16. 13^ with 
the appropriate scaling results given in Prop. 16.141 



2 



Mathematical preliminaries 



In this chapter we present some necesssary mathematical prehminar- 
ies, mostly related to probabilistic analysis of functionals of PPPs. The 
most important example of such a functional is the aggregate interfer- 
ence experienced by a typical Rx in a wireless network when the loca- 
tions of interfering nodes form a PPP and the channel is distance depen- 
dent. Many of the results in this chapter are also found in the excellent 
monographs by Haenggi and Ganti ^36j and Baccelli and Blaszczyszyn 
[U |9]. Our treatment of this large field is quite selective: we present 
only those results directly relevant to computing the OP and TC. We 
recommend both these monographs for a more in depth treatment of 
application of the mathematical field of stochastic geometry to the per- 
formance analysis of wireless networks. To the extent possible we have 
used notation consistent with that used in [36l[8l[9]. Moreover, when- 
ever possible we give references in [521 El E] to corresponding results 
presented in this chapter. 

Denote the reals by M, the natural numbers by N = {1, 2, 3, . . .}, 
the integers by Z, and the complex numbers by C (and \f—\ by i). 
We use = for equality that holds by definition. We work in M*^ where 
d is the spatial dimension of the wireless network. Our analysis holds 
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for general d G N, but d G {1,2,3} are the relevant cases. A point in 
R'' is denoted by x = (xi, . . . ,Xd)- The Euclidean (L2) norm ||a;||2 = 

is denoted by \x\, and the origin (0, . . . , 0) by o. We denote 
the natural log as \ogx. For a natural number G N, we write [N] for 
the set {1, . . . , N}. We use the shorthand aAb = min{a, b} and a\/ b = 
max{a, b}. We use standard asymptotic order notation O(-), &{•). 
We begin with the d-dim. ball and annulus and their volumes. 

Definition 2.1. Ball and annulus. The d-dim. ball {d G N) centered 
at c G M'^ with radius r G M+ is: 

hd{c,r) = {xGR'^:\x-c\<r}. (2.1) 

The d-dim. annulus centered at c G M*^ with radii < ri < r2 is: 

&d{c,ri,r2) = {x G M"' : n < |x - c| < rs}. (2.2) 



Volume of a set 5 C M'^ is denoted IS*]. The ball and annulus volumes 
are given below (f35] (3)). 

Proposition 2.1. Ball and annulus volume. The d-dim. ball 

brf(c, r) and annulus adic,ri,r2) have volume 

\hd{c,r)\=Cdr'^, |arf(c, n, ra)! = ^(r^ - rf), (2.3) 

where 

d 

Cd = i (2.4) 
^vr^2'^(^)!, dodd. 



The relevant values of Cd are: 

4 , , 

ci = 2, C2 = vr, C3 = -TT. (2.5) 

We will also have frequent use for the gamma function. 
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Definition 2.2. Gamma function. The gamma and incomplete 
gamma function are, respectively, 

j-co r-th 

T{z)= e~^e-*dt, r{z,ti,th)= t'-^e-*dt (2.6) 

^0 Jti 
for z G C and < ti < th < oo. 

Note r(z,0,oo) = T{z) and that r{z) = {z - 1)! for z G N. We wiU 
have use for the identity: 

r(i-5)r(i + 5) = -^. (2.7) 

and the fact that r(-l/2) = -2^/^. See Fig. Ofor (pTj) . 

2.1 Probability: notations, definitions, key inequalities 

The material in this section is quite standard and is available in most 
textbooks on probability. The following is a key notational convention. 

Remark 2.1. RV notation. We indicate random variables (RVs) 
with a sans-serif non-italic font, e.g., x, h,u,c, and their realizations 
(as well as other non-random quantities) with an italicized serif font, 
e.g., X, h, u, c. A notable exception is the use of 11 (Def. II. 2p to indicate 
a random point process. 

Standard probabilistic quantities are denoted as follows. 

Definition 2.3. Standard probability definitions. Let x denote a 
continuous real-valued RV, and let t G M, G and s G C. 

(1) The cumulative distribution function (CDF) is -Fx(i) = IF'(x ^ 
t) for t G M. Denote the CDF for random x by x ~ F^. 

(2) The complementary CDF (CCDF) is F^{t) = 1 - Fx(t) = 
P(X > t). 

(3) The inverse CDF and inverse CCDF are F~^{p) and F^^{p) 
for p G [0, 1]. 
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(4) The probability density function (PDF) is /x(t) = ^i^x(i)- 

(5) Expectation is denoted by E[x], variance is denoted Var(x). 

(6) The Laplace transform (LT) is £[x](s) = E[e~'*'*], note s £ C 

(7) The characteristic function (CF) is </'[x](t) = E[e'*'*], note 
t G M. 

(8) The moment generating function (MGF) is M[x]{9) = E[e^'*], 
note 61 e M+. 

(9) The hazard rate function (HRF) is7i[x\{x) = logFx(x) = 

AM 
Px(x) • 

(10) A normal RV with E[x] = /i and Var(x) = o"^ is denoted 
X ~ N(fi, a). A standard normal is denoted z ~ N{0, 1) with 
CDF F^{t) and CCDF F^{t). 

(11) Equality in distribution between RVs x,y is denoted x = y. 



The LT with argument s € C is more general than both the CF 
and the MGF, but the LT and MGF need not exist, while the CF 
is guaranteed to exist. When all three exist, the CF and MGF are 
obtainable from the LT: 

M{t) = C[x]{-it), MM{9) = C[x]{-9). (2.8) 

We will have use for Jensen's inequality {e.g., Cor. 14. 3p . 



Proposition 2.2. Jensen's inequality. For a RV x, if / is a convex 
function then E[/(x)] > /(E[x]), with equality holding for / affine. 



The three inequalities of Markov, Chebychev, and Chernoff are each 
UBs on tail probabilities. The three inequalities build upon one an- 
other. In general it is fair to say that Markov is simpler to apply than 
Chebychev, and in turn Chebychev is simpler to apply than Chernoff. 
This is on account of the fact that Markov relies only upon the mean 
E[x], while Chebychev depends upon the variance Var(x), and Chernoff 
is a function of the moment generating function A^[x](s) = E[e*'^]. In 
general (but not always) it is further the case that the Chernoff bound 
is tighter than the Chebychev bound, and the Chebychev bound is 
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tighter than the Markov bound. The tightness of the Chernoff bound 
also comes about through the flexibiUty to tune the free parameter s. 
These inequahties wih be appUed in ^3.51 to derive an UB on the OP. 
We begin with Markov's inequahty. 

Proposition 2.3. Markov's inequality. For a nonnegative RV x and 
t G M+: 

P(x > t) < (2.9) 



Proof. Define the BernoulU indicator RV lx>t and observe x > tlx>t 
for all t G M+. Taking expectations yields 

E[x] > tE[lx>t] = tP(x > t). (2.10) 

□ 

Chebychev's inequality is obtained by applying Markov's inequality to 
the nonnegative RV |x — E[x]|. 

Proposition 2.4. Chebychev's inequality. For a RV x and t G M+: 

P(|x-E[x]| >t) < (2.11) 



Proof. Apply Markov's inequality with |x — E[x]|: 

P(|x - E[x]| >t)= P((x - E[x])2 > t") < ^[(^~^^M)'] . (2.12) 



Finally, Chernoff 's inequality is obtained by applying Markov's inequal- 
ity to the nonnegative RV e^'^. 
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Proposition 2.5. Chernoff's inequality. For a nonnegative RV x 
and t G M+: 

P(x > < inf E[e^^e-^* (2.13) 

6»>0 



Proof. Observe the equality of the events {x > t} and {e^^ > e^*} for 
all > and apply Markov's inequality to the nonnegative RV e^'^ 

P(x >t)=F i^e^"" > e^*) < E[e^^e-^*. (2.14) 

The above inequality holds for all ^ > and hence in particular for 
that 6 that minimizes the UB. □ 

2.2 PPP void probabilities and distance mappings 

Recall Ild,\ = {xj} denotes a PPP with points {xj} C M.'^ of intensity A. 
The two most important examples for us are 112, a and IIi^i. We often 
will write Hi i = {tj} in accordance with the usual interpretation of 
the points in a one dimensional point process as times. 



Assumption 2.1. Labeling convention for PPP. All point pro- 
cesses are assumed to number points in order of increasing distance 
from o: Ild,\ = {xj} with |xi| < |x2| < • • • . 



We present two key facts about distances for PPPs in this section. 
First, the void probability P(|xi| > r) is the probability that there are 
no points from 11^^ in the ball bd(o, r), i.e., that the nearest neighbor 
to o in Ild,x is at least at distance r. 



Proposition 2.6. Void probability. The RV |xi| has distribution 

F{Ud,x{hdio, r)) = 0) = P(|xi| > r) = e-^'^'*^', r e M+ (2.15) 
where Cd is defined in Prop. 12.11 



Proof. In words, |xi| > r is the event that the nearest point in PPP 
Ild,x to o is at least a distance r away. This is the same as there being 
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no points in the PPP lying in the ball bd(o, r). Thus, |xi| > r 44> 
Ud,x{hd{o,r)) = 0. Recall from Def. Othat the RV Ud,x{hdio,r)) is 
Poisson with intensity A|bd(o, r)| = Xcdr'^, and thus: 

P(|xi| > r) = F{nd,x{hd{o,r)) = 0) = e"^^''^'. (2.16) 

□ 

For d = 1 we recover the elementary fact that |xi| ~ Exp(2A), i.e., 
P(|xi| > r) = e"'^^'' and for d = 2 we have P(|xi| > r) = e~'^'"'^ Prop. 
12.61 will i) form the basis for the max SN distribution (Cor. 12. ip which 
in turn will yield the LB on the OP in §3.3( and will be generalized 
to a) a non-homogeneous marked PPP (MPPP) in Prop. 14. 3( and in) 
distances to the k^^ nearest neighbor in Thm. 15.11 

The second result in this section is a special case of a more general 
mapping theorem given below ([36]. Thm. A.l, p. 107 and [49| §2.3): 



Theorem 2.1. Mapping theorem ([H^ Thm. A.l). Let ^ be an in- 
homogeneous PPP on with intensity function A, and let / : M*^ — )• 
be measurable and A{f~^{y}) = for all y G M*. Assume further that 
fi{B) = A{f-^{B)) satisfies n{B) < oo for all bounded B. Then /(^>) 
is a non- homogeneous PPP on with intensity measure fi. 



We refer the interested reader to |36j for the formal definition of inho- 
mogeneous PPP, intensity function, and measurability. The following 
proposition is a special case of Thm. 12.11 



Proposition 2.7. Distance mapping. Let 11^ a = {xj} be a PPP in 
M'^ of intensity A, and IIi^i = {tj} a PPP in M of intensity 1. Then: 

\cd\x,\'^ = 2\ti\, ien. (2.17) 



Proof. Consider Thm. 12.11 for s = 1, <I> = 11^^ a homogeneous with 
intensity A(^) = X\A\ for all compact ^ C M"^ for some A G M+, and 
f{x) = ^|x|'^sign(x^"^^), where x^^^ is the first component in vector 
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X = {x^^\ . . . jX^")) and sign(a) = la>o — la<o is the sign of a G 
Consider bounded symmetric intervals B of the form [—t, t] for t E B 

A{rH[-t,t])) = x\f-\[-t,t])\=x\{xeR'':f{x)e[-t,t]} 

xeM.'^: ^lxl'^sign(t) e [-t,t] 
X e M"' : l^l < (-^ 



A 



2t 



(2.18) 



By the mapping theorem t, i]) = 2t, which is to say that f{Ild,x) 
is a homogeneous PPP of unit intensity, i.e., Hi^i. □ 

In particular, for d = 2 this result states 7rA|xj|^ ~ 2|ti|. Prop. EB] and 
12.71 are easily seen to be consistent for z = 1 in that they both give: 



— |xi| >r 



|xi| > 




-2r 



|ti| > r) 



(2.19) 

Prop. 12.71 is somewhat analogous to the standardization of normal 
N{fi,a) RVs to A^(0, 1), i.e., for x ~ N{ij,a) and z ~ A^(0, 1) the 
standardization of x is (x — /i)/cj, which is equal in distribution to z. 
Prop. 12.71 is used below in Prop. 12.81 for mapping probabilities associ- 
ated with functionals of distances in 11^,;^ to probabilities associated 
with functionals of distances in IIi^i. 



2.3 Shot noise (SN) processes 

Consider a system given injections of energy or noise at a sequence of 
random times (shock times), where each energy injection attenuates in 
time according to an impulse response function, so that the random 
cumulative energy seen at any given time t is the superposition of 
attenuated shocks from all shock instances prior to t. Such a process 
is termed a (temporal) shot noise (SN) process, and was first used 
by Schottky in 1918 [65j to explain how transfers of charge at random 
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Fig. 2.1 Top: a temporal SN process for d = 1 is the superposition of appropriately at- 
tenuated discrete time injections of energy / noise into a system. Bottom: a spatial SN 
process for d = 2 in the context of a wireless network is the superposition of interferences, 
where each interferer is appropriate attenuated through the corresponding channel. 



discrete units in time in vaccuum tubes give rise to current fluctuations. 
Fig. O shows a sample SN process. See e.g., P El [Ml ES] , 09] (Ch. 3), 
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|21j (Vol. 1, Ex. 6.1(d)) for more information. The following is a general 
definition of a SN process for arbitrary dimension d. 



Definition 2.4. SN process. A (sum) SN process is a real-valued 
random process {5I(x)}, indexed by the continuous parameter x G M'^, 
that is a functional of an underlying (stationary) point process 11 = 
{xj} C R'^, where 

ien 

Here I : M4- — ?• M-j. is a linear time-invariant impulse response function 
and {hj} is a collection of i.i.d. nonnegative RVs. A max SN process 
{M(a;)} formed from H, / is 

M^(x) = maxhiZ(|xi - x|), x G R'^. (2.21) 

ien 

Note M\j{x) < "^{lix) a.s. for each x G M'^. 



Remark 2.2. SN index convention. Throughout this volume we 
write functionals of PPPs by summing over their indices rather than 
their points, e.g., X^ien /(^O instead of X^x^en /(^i)- Although the lat- 
ter is maybe clearer in this case, it becomes awkward for marked point 
processes, say $ = {(xj, mj),z G N}, with points {x,} and marks {rrii}. 
In this case, writing /(^i; ^i) is more clear and compact than 



The case d = 1 is most common in the stochastic process literature, but 
the case d = 2 is most relevant for spatial models of wireless networks. 
The interpretation of /(|xj — x|) for d = 1 is the energy injected at 
time Xj attenuated over the time interval [xj,x], and thus T.\i{x) is 
the superposition of all time-attenuated energy injections seen at time 
X. The interpretation of /(|xj — x|) in the context of a d-dimensional 
wireless network is the interference generated by the node at position 
Xj attenuates in space over the distance |xj — x| at position x, and 
thus ^\j{x) is the superposition of all distance-attenuated interferences 
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seen at position x. The simple LB M{x) < will be shown to 

be asymptotically tight in Prop. 12.131 and will form the basis for the 
various LBs on OP and the UBs on TC. 

Remark 2.3. Radial symmetry. We have restricted Def. 12.41 to ra- 
dially symmetric functions Z(r) for r = \x\, based on the assumption 
that the impact on x of a shock at Xj depends only on the distance 
|x — Xj|. This assumption allows all integrals of / over M*^ to be replaced 
by integrals over R_|_ (see Thm. [212]) . 

For our purposes it suffices to consider a rather specific case. 

Assumption 2.2. Power law impulse response. Assume the fol- 
lowing for the SN process {T'^ix)} in Def. El 

(1) The impulse response function is a power-law truncated 
around o 

W{r) = r-"l,>„ r £ R+, (2.22) 

for a > 0, e > 0; 

(2) The stationary point process 11 is a PPP Ild,\', 

(3) The amplitude RVs {hj} are all unity; 

(4) We restrict our attention to the origin T.\j{o). 



The assumption on the amplitudes {hj} will be relaxed in §4.11 We 
will employ a special notation for T.\j{x) under Ass. 12.21 

Definition 2.5. Power law SN and characteristic exponent. 

The SN RVs at o under Ass. 12.21 are denoted 

^n(o) = ^fc,(o) = Z^;^(o) 
Mii(o) = Mj,- (o) = M^;^(o) (2.23) 
The characteristic exponent of 5I^'^(o) is defined as 

6=-. (2.24) 

a 
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Remark 2.4. Pathloss attenuation and the singularity at the 
origin. The impulse response function lafli''') = is a common 
choice to model the attenuation due to pathloss in wireless commu- 
nication but suffers the drawback of modeling amplification \x\~°' > 1 
rather than attenuation of received energy at distances |x| < 1, and in 
fact this amplification grows without bound as |x| — )• 0; this is further 
discussed in Prop. 12.91 below, in [SBJ (p. 24), and in [M]. Generaliz- 
ing lafi by truncating around the origin la^e removes this singularity. 
Note la,€ is used in [36] (§3.7.1) to model carrier sense multiple access 
(CSMA). 



The max SN RV M^'^(o) will be shown to obey the Frechet distribution, 



defined below. 



Definition 2.6. The Frechet distribution with parameters 7 > 0, 
(7 > 0, and /Li G M has CDF 

Fx(x) = exp |- "I , X > /i, (2.25) 

and for ^ = has moments up to order 7: 

E[xn = | ^'"^^'-^Z^^' ,f^ = 0. (2.26) 

[ 00, else 

The Frechet is one of three extreme value distributions 1501. 



The max SN RV M^'^(o) CDF is immediate from Prop. 12.61 



Corollary 2.1. The max SN RV CDF is 



For e = the max SN RV has the Frechet distribution in Def . 12.61 with 
7 = (5, o" = (Xcd) and fi = 0. 
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Having characterized M^'^(o), we now focus on characterizing the RV 
Z^'^(o), as it represents the aggregate interference seen at a typical loca- 
tion when interferers are positioned according to IT^^a and the pathloss 
attenuation function l^^^ is assumed. Cases of particular interest are 

Z^;°(o),Z^;^(o),Z?;°(o),Z?;?(o). (2.28) 

Many results will hold provided the characteristic exponent 5 < 1. For 
the important case d = 2 this translates to a > 2. Prop. 12.71 is used 
below to show that it suffices to consider 5l"'^(o). 



Proposition 2.8. Interference mapping. The following RVs are 
equal in distribution 

^d,xio) = I — 1 (o), and Z^\(o) = I — 1 (o). 

(2.29) 



Proof. Using Prop. 12.71 gives: 

(2.30) 

For T.^'^{o) simply observe {|xj| > e} = {|tj| > Acde'^/2}. □ 

Prop. 12.81 is important because it expresses the SN RV formed from 
n^^A with exponent a as a scaling of a SN RV formed from IIi^i with 
exponent a/d. In this sense d and A are inessential parameters. 

The next result is called the Campbell-Mecke Theorem; the version 
below is a special case of a much more general theory on moments of 
functionals of PPPs (see e.g., Thm. A. 2 and Lem. A. 3 in |36j). Our 
specialization is to homogeneous PPPs with measure A(dx) = Adx, 
and to radially symmetric functions /(|3;|). As mentioned in Rem. 12.31 
and 12.41 this assumption is natural for wireless networks, and has the 
advantage of allowing the d-dimensional integrals to be replaced with 
single dimensional integrals using the following theorem (from [TOj). 
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Theorem 2.2. Integration of radially symmetric functions 
(|10]). Let / : M+ — M be Riemann integrable on M_|_, and let 
l{x) = l{\x\) for X G M'^. Then / is Riemann integrable on R"' and 

r poo 

/ l{x)dx = dcd / /(r)r'^-Mr. (2.31) 
Jw^ Jo 



The above theorem is used to simphfy the integrals in the Campbell- 
Mecke Theorem below and in several other places throughout this 
monograph. 



Theorem 2.3. Campbell-Mecke ((HH] Thm. A.2, [8j Thm. 1.11, 
§3.2). The mean and variance of the RV T.\j{y) in (12:20]) for PPP U^^x 
and measurable / : M+ — )■ M+ are: 

Epii^Jo)] = AE[h] / /(|a:|)dx = A(icrfE[h] / /(r)r'^-Mr 
Var(Z[i^^(o)) = AE[h2] / l{\x\fdx = XdcdE[\^'^] ;(r)V~Mr 

jRd. Jo 

(2.32) 



The proof of Thm. [231 is essentially an exchange of the order of integra- 
tion and summation and is omitted. In particular, for la,oix) = 
we change variables from for x G to r E M+ to exploit the radial 
symmetry of the function /: 

E[Z"'°(o)] = X [ |x|-"dx = Xdcd r r~^r'^~^dr = ^^r'^"" . 

Jo d-a Q 

(2.33) 

As discussed in [36j ((3.4) p. 24) this integral diverges for a < d due to 
the upper limit of integration, and for a > d it diverges due to the lower 
limit of integration, which in turn is attributable to the singularity at 
the origin of the function r~". As stated in Ass. 12.21 we use la,t{f) = 
r~°'lr>e for r G M+, which can be interpreted as assuming a Rx has 
perfect interference cancellation within a ball of radius e. The following 
proposition summarizes this discussion. 
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Proposition 2.9. SN mean and variance. The means and vari- 
ances of the SN RVs (l2:23D are: 

M j oo, a < d or e = 

n^,,xio)] = i ^^^^ 



M w J Q<(i/2ore = 



The proof is a straightforward modification of (j2.33|) . Prop. [2^91 will be 
used along with the Markov and Chebychev inequalities in §2.1l to form 
the LBs on TC (UBs on OP) in g331 

The last SN specific result we will have use for is the series expansion 
of the PDF and CDF of a SN RV for d = 1. The following result is 
adaptecflfrom [57J (Eq. (29)). 

Proposition 2.10. SN series expansion (j57j). The series expan- 
sions of the PDF and CCDF of the RV for = ^ < 1 are: 

4.,o(y) = — V ^ > r(l + n5) sin(7rn,5)(2Ar(l - 6)y-'r 

■^i,A Try ^-^ n\ 

" n=l 

F^.,o(y) = --V ^ > T{l + n5)sm{^n5){2\T{l-5)y-'T 

n=\ 

(2.35) 



The asymptotic PDF and CCDF as y — )• oo is immediate upon taking 
the dominant n = 1 term from the above expansions. 

Corollary 2.2. The Asymptotic PDF and CCDF of the SN RV 

11°'^ as y — )• oo for 5 = ^ < 1 are: 

/^.,o(y) = 2\5y~^-^ + 0{y-^-^^), y ^ oo 
i,A 

F^.,o(2/) = 2Xy-^ + 0{y-^^), y ^ oo (2.36) 
i,A 

^ Eq. 1 12. 351 1 have a factor of 2 in front of A not present in [571 (29) due to the fact that their 
impulse response function (4) does not count contributions from t < 0. 
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Employ ()2.7p to simplify the n = 1 term. Cor. 12.21 gives the asymptotic 
approximations of OP (as A — ?■ 0) and TC (as q* — >■ 0) in Prop. 13.31 

2.4 Stable distributions, Laplace transforms, and PGFL 

A RV is said to be stable if iid sums of that RV are equal to an affine 
function of the original RV (closure under summation). Equivalently, a 
distribution is stable if convolutions of the distribution yield a trans- 
lation and/or scaling of the distribution (closure under convolution). 



Definition 2.7. Stable RV and distribution. Let x ~ F and 

(xi,...,Xn) be iid from F. Say x is a stable RV (F is a stable dis- 
tribution) if for each n G N there exists numbers (a„, bn) such that 

anX + 6„ = xiH hx^. (2.37) 

Moreover, if it exists, a„ = n^^^ for a characteristic exponent 5 G [0, 2]. 
If 6„ = for all n then x (F) is a strictly stable RV (distribution). 

See e.g., [59] Def. 1.5. Perhaps the most familiar example of a stable 
distribution is the normal: let (x, xi, . . . ,x„) ~ N{fi,a) be independent 
normal RVs with mean ^ and standard deviation a. Then xi -|- • • • -|- 
x„ ~ N{nfi, ^Jna^ and ax + 6 ~ N{aii + 6, aa). Choose an = ^/n and 
bn = (n — -y/n) fi to satisfy the requirement in Def. 12. 7[ 

Aside from a few special cases (the normal, Cauchy, and Levy dis- 
tributions), a stable RV x does not admit a closed form CDF. It does, 
however, have a special form for its CF 0[x](t) = E[e'*^] for t G R. 
We are interested in a specific sub-class of stable RVs appropriate for 
modeling SN RVs, and consequently the definitions in the remainder 
of this section are specialized to that class. 

Definition 2.8. Stable CF. The RV x is stable with characteristic 
exponent 5 G (0, 1), dispersion coefficient 7 > if it has CF 

(/)[x](i) = exp |-7^|t|'^(l - itan(7r5/2)sign(t))| , t G M. (2.38) 
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The above definition is adopted from Def. 1.8 from ^59j. It is special- 
ized in that we have fixed the location parameter at and the skewness 
parameter at its maximum value of 1 ("totally skewed to the right" in 
[59] p.l2). The support of this RV is M+ (j59] Lem. 1.10). As discussed 
in [59] §1.3, there is a wide variety of parameterizations for stable dis- 
tributions found in the technical literature, and their differences may 
easily lead to confusion. 

The case 5 = ^ corresponds to the Levy distribution ([59] §1.1). 

Definition 2.9. The Levy distribution with parameter 7 > and 
support E+ has PDF, CDF, and CF 

I — 7 

/ 7 e 2a; 

/x(a;) = \h —, X eR+ 

V ^ 2 

F,{x) = 2Fz (J^) , X GM+ 



</.[x](t) = e-^^\ t EM (2.39) 
where is the normal N{0, 1) CCDF. 

Observe from (j2.39p and Def. I2.8l that a Levy RV is stable with charac- 
teristic exponent 5 = ^ and dispersion coefficient 7. To see the equiva- 
lence between (lOHl) with 6 = 1/2 and (12:39]) note: 

tan(7r/4) = 1, 0[x](t) = j ' ^ ^ ^ , 1 ± i = V±2i. 

(2.40) 

Fig. [22] shows the Levy PDF and CCDF. Note the heavy tail in the 
right plot compared with the light tailed normal distribution. The char- 
acteristic exponent 6 fixes which moments of a stable RV are finite. 

Proposition 2.11. Stable moments. For x stable with characteris- 
tic exponent 6 and p G M+ : 

r <oo, p<6<2 _ 
[ = 00, else 
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Fig. 2.2 The PDF (left) and CCDF (right) for the Levy distribution with parameter 7 S 
{0.5, 1,2}. The CCDF is on log-log scale; the normal |Af(0, 1)| CCDF is shown as well. 



For example, a Levy RV x has E[x^] = oo for p > 5 = 

The notation for the characteristic exponent 6 in Def. 12.81 has been 
dehberately chosen to coincide with the notation for the characteristic 
exponent 6 of T.'^'^{o) in Def. 12.51 Our objective in this section is to 
characterize conditions under which Z^'^(o) in ()2.23p is stable. Cor. 12.61 
will show 5I^'^(o) is stable with characteristic exponent 5 = d/a < 1. 
To get to this result, we must find the CF of 5I^'^(o) and compare 
with Def. ESI We will show that the LT of a SN RV i:[^{o) in (|2:20]1 is 
expressible in terms of the probability generating functional (PGFL) 
of the underlying point process 11, and that this PGFL admits a closed 
form for PPP Ild,x and the truncated power law function /^.e (|2.22|) . 
We start by defining the PGFL ([36] Def. A. 5). Our methodology in the 
development that follows is to present results in a somewhat general 
form and then specialize as needed. Thus we distinguish between i) 
a stationary point process H and a PPP Ild,\, H) generic measurable 
functions l,iJ and the specific function la,e (|2.22|) . and SN RVs ^\j{o) 
dl^OD, Jo), and ^^'^(o) I^^M- We aim to clarify the impact of 
the assumptions of a PPP Ild,x and a particular pathloss function l^^^. 



Definition 2.10. Point process PGFL. The probability generating 
functional (PGFL) of a point process H and a measurable function 
v :R'^ ^ Ml is defined as 



(2.42) 
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If the point process n is a PPP Ud,x the PGFL simplifies (p] (A.3)). 
Proposition 2.12. PPP PGFL. For a PPP Ud,x the PGFL is 

gpd.A,!^] =exp|-A^^(l-z.(x))dx|. (2.43) 



The PGFL yields the LT of a SN RV. 



Corollary 2.3. Point process SN LT. The SN RV P^{o) in (IMj) 
for a stationary point process H with each hj = 1 has a LT expressible 
in terms of its PGFL: 

/:pUo)](s) = g[n,e-^'(i-i)], sec. 



(2.44) 



Proof. 



E 



s5:{j(o) 



E 



exp < -s^/(|xi|) 
I ien 



E 



.ien 



(2.45) 



When the point process H is a PPP Ild,\ we can combine Prop. 12.121 



and Cor. O to get the LT of I^^ ^ . 



Corollary 2.4. PPP SN LT. The SN RV ^(o) in (ICT]) for a 

PPP n^.A with each hj = 1 has a LT 



C[T^n,,M{s) = exp |-Adc, (l - e^'W) r'^-Mrj , (2. 



46) 



for all s G C for which the integral exists. 



Proof. Prop. EH and Cor. O yield: 

^Pn,,,(o)](s)=exp|-A^^ (l - e 
Now change variables from |x| to r using Thm 



Sim) 



(2.47) 
□ 



We next fix / to be la,e and obtain the MGF of 5:°'^(o). Recall 

that we can obtain the MGF from the LP SIM. 
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Corollary 2.5. Pathloss SN MGF. For (5 < 1 and e > the MGF 

of 51^'^ (o) in Def. Ois (for 9 eR+): 

-MK:A(o)]W=exp|^^' {e'y-l)y-^-^dy\. (2.48) 



Proof. Substituting la,e (l2:22]l into (f2:i6D gives 

C[Tl{{o)]{s) = exp |-Adc, (l - e-^^"") r^-'dr^ . (2.49) 
The change in variable: 

y = r^'^^r = y-^^r'^'^ = y'^,dr = -iy^^^My, (2.50) 

a 

yields 

C[TZiom=e^j>i^-^iy {l-e-^y)y-'-'dyy (2.51) 

Specializing to the assumed s = —9 for 9 G M+ yields the proposition. 
For (5 > 1 or e = this quantity diverges. □ 

Finally we fix / to be 1^,0 (IT^ and obtain the CF of ^|J'°(o) (fT^ . 
Recall that we can obtain the CF from the LP (j2.8p . 



Corollary 2.6. Pathloss SN CF. For 5 < 1 the CF of ^^');(o) in 
Def.ESlis (for t G M): 

(/.p^'°(o)](t) = exp ^\5cdT{-5) cos(7r5/2)|t|^ (1 - itan(^(5/2)sign(t))} . 

(2.52) 

In particular, 51^'^ (o) is stable as in Def. 12.81 with stability coefficient 
6 <1 and dispersion coefficient 

/ = -\ScdY{-6) cos(7r(5/2). (2.53) 
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Proof. Fix e = and s = —it for t G R in (I2.5ip . and assume 6 < 1: 



^■0(o)](t) = exp|-A5cd|t|^-^r(-(5) (-|t|cos(7rV2) +itsin(7r5/2))} . 

(2.54) 

Straightforward manipulations give (|2.52p . □ 
For the special case of 5 = ^ we can apply Def. 12.91 

Corollary 2.7. Pathloss SN for 6 = ^. Fox 6 = i, (f233|) gives 7 = 
|(Ac^)^ and the RV Z^'^(o) is Levy with parameter 7. In particular, 
7 = 2^ for Y.lj{o). 

2.5 Maximums and sums of RVs 

To finish this chapter, we combine several previous results to illustrate 
the asymptotic tightness of the LB M|^'^(o) < 5I|^'^(o). 

Proposition 2.13. Sum and max SN CCDF ratio. For e = and 

(5 < 1 the CCDFs of the RVs M;J'°(o) and T.'^''^{o) have a ratio that 
converges to unity: 

P(Z"f (o) > y) 
lim ^ '^"^^ ^ ^' = 1. (2.55) 
y^oop(M^-0(o)>2/) 



Proof. First apply Prop. [2^ and then apply Cor. 12.21 to Y.'^'^{o): 

V ^ 

= Ac^y--^ + ©(y-^-^) (2.56) 

Using Cor. 12. 1^ the first order series expansion of the CDF of M^'^(o) 
is: 

^ {K:xio) > y) = 1 - e-^'"''' = ^c,y'' + 0{y-^'). (2.57) 



K:iio) > ( ^ ) y 
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The limit of the ratio of the CCDFs ()2.56p and (I2.57|) as y — )• oo is one. 

□ 

Fig. ED shows the (exact) CCDF for Z (for 5 = 1/2) from Cor. 
the asymptotic CCDF for 5:^'°(o) from Cor. EZl and the CCDF for 
M^'^(o) from Cor. EH Observe the ratio of the CCDFs appears to 
converge to one as y — >• oo. This convergence will be used to establish 
the asymptotic tightness of the OP LBs and TC UBs in what follows. 

CCDF 




0.2 0.4 0.6 0.8 1.0^ 

Fig. 2.3 The exact and asymptotic CCDFs for the RV 51^'^ (o) and the CCDF for the RV 
M^'°(o). The parameter values are e = 0, A = 1/10, d = 2, and q = 4 (5 = 1/2). The 
M|^'°{o) CCDF is asymptotically equal to the 5I^'°{o) CCDF as j/ -> oo. 

The relationship between the max and sum of a sequence of RVs 
has a long history in the literature on probability. Levy (1935) |54j . 
Darhng (1952) [22], Chistyakov (1964) [TT], and Chow and Teugels 
(1978) |19j are important early works. Somewhat more recently, Goldie 
and Kliippelberg (1997) [29] characterize the class of "subexponential 
distributions" (first introduced in [17J). 



Definition 2.11. Subexponential distribution ([29j, Def. 1.1). Let 

{xj} for i G N be iid positive RVs with CDF F such that F{x) < 1 for 
all X > 0. F is a subexponential CDF if one of the following equivalent 
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conditions holds: 



lim ^("^ + • • • > = n, Vn>2 (2.58) 

x-5>oo F[X) 

lim "\ = 1, Vn>2. (2.59) 



Note in particular ()2.59p states sub exponential distributions have the 
ratio of the CCDF of the sum over the CCDF of the max approaching 
one asymptotically, which is precisely result in Prop. 12.131 For our 
purposes we require only one of their results, a sufficient condition 
for a distribution to be subexponential, which we condense and adapt 
below. Recall the hazard rate function T-L(x) in Def. 12. 3[ 

Proposition 2.14. SufBcient subexponential condition ([29], 
Prop. 3.8). 

limsupy'H(y) < oo =^ F is subexponential . (2.60) 



A simple and natural way to apply Prop. [2TT^ to our case is to condition 
on the number of nodes N from n^^A within a bounded domain, say 
B C M.'^ . More formally, suppose Ild,x{S) = In this case the A'' 
points are independent and uniformly distributed on B, and form a 
so-called binomial point process (BPP) [36j (§A.1.1). 

Definition 2.12. Fix A^ € N and bounded B C W^. The binomial 
point process (BPP) IIb,n = {xii • • • ,^n} consists of A^ points in- 
dependently and distributed uniformly at random in B. 

We fix S = brf(o, R) for ReR+, and derive the CDF and HRF for the 
interference contributions from each of the nodes in a BPP seen at o. 

Lemma 2.1. BPP distances and interference. Let Hr^n = 
{xi, . . . ,XAr} be the BPP on bd(o, R). The CDF for the RV |xi| is 

- ' ^<r<R ^2.61) 
1, r > R 



P(|x,| < r) 
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Assuming e = 0, the interference contribution RV Ij = -P|xj| " from 
each node i has support Ij G [ymin, 

Oo) for 2/niin — 

Pi?-" and CCDF: 



1 ' 



P(|, > y) =p f |xi| < f-V 1 = <! ^' y-y^^^ (2.62) 

\ \y J J [1, y< ymin 

The hazard rate function of I is 

n[\]{y) = S^ = l f' y>y-- . (2.63) 

' F{y) \ 0, else ^ ^ 

The fractional order moments are 

E[\P] = { 5-p ' (2.64) 

oo, else 



Applying Prop. [27T^ to the HRF (I2.63P gives that I is subexponential. 

Corollary 2.8. Subexponential BPP interferences. The RV I is 

subexponential as the HRF (j2.63p for I in Lem. [2TT] satisfies Prop. [2TT^ 

lim yn[\]{y) = S < oo. (2.65) 



Prop. [2l3] (for the PPP) and Cor. EE (for the BPP) both demonstrate 
that the sum and max RVs of the interference contributions under 
the pathloss model /(|x|) = \x\~°' for 6 £ (0,1) have CCDFs that are 
asymptotically equal. More succinctly, the probability of the sum be- 
ing large is roughly the same as the probability of the max being large. 
Large sums occur due to a small number of large individual contri- 
butions; they do not occur due to a large number of small individual 
contributions. This intuition helps explain why the LB on the OP and 
UB on the TC, which are ultimately derived from the simple bound on 
the interference RVs M < Z (Def. [23]) . are asymptotically tight. 



3 



Basic model 



In this chapter we consider the most basic model for computing the 
OP and TC. The wireless channel between any two nodes consists of 
pathloss attenuation with no fading. As indicated in Def. 11.11 and II. 3^ 
the key quantity is the SINR, defined below. 

Definition 3.1. Basic model SINR. Under the basic model, the 
SINR seen by a reference Rx located at o when all nodes use constant 
power P, the interferers form a PPP 11^ ,^, the noise power is A^, the 
channel model is la,t{T~) as in Ass. \2.2\ and each Tx is positioned at a 
fixed distance u> e from its Rx is 

sinr(o) ^ (3.1) 
where the received signal and interference powers are 

S = PLA^) = Pu-'^ and no) = PW{\^i\)- (3-2) 



We emphasize the only random quantity in sinr(o) in ()3.ip is 51 (o), 
and the only random quantity in 51 (o) in (13. 2p is the PPP 11^^ a- The 
following quantity will be used frequently throughout this chapter. 
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Definition 3.2. Rx SNR. Define ^ = - f j ° . Define the Rx 

SNR snr = 



Considering the = case = ut~) makes plain that ^ has units of 
meters. To simpHfy the analysis that follows it is convenient to make 
the following assumptions. 



Assumption 3.1. SNR LB. The Rx SNR obeys snr > r. Moreover, 
assume e < 



The first assumption states that the received SNR exceeds the SINR 
threshold, i.e., in the absence of interference a transmission attempt is 
successful. As will be clear below, the second assumption states that 
dominant interferers are possible. We compute the exact OP and TC 
in §3.11 asymptotic exact OP and TC for A — )• and — )• in §3.2^ 
UBs on TC (LBs on OP) in g331 and LBs on TC (UBs on OP) in M 
Several extensions on this basic model are presented in Ch. HI 

3.1 Exact OP and TC 

The next result gives the OP in terms of the CCDF of the SN RV 
representing the aggregate interference seen at o under the PPP IIi^i. 

Proposition 3.1. OP is SN CCDF. The OP for the SINR in Def. 
13.11 is expressible as the tail probability of a SN RV ^i^i^ (o) on IIi^i 
evaluated at y = (Acrf/2)~^^~": 

g(A) = P (rlf-'io) > y) , 6-' = ^> = Xc,e''/2. (3.3) 



Proof. By manipulation of the outage event {sinr(o) < r} and employ- 
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ing Def. EH and Prop. E 

Z(o) >--N 

T 



2^ Wdxii) > — -- 



= {^'I(«)>r"} 

= {(Acd/2)i/^zi/f'"^''^'/'(o) > r"} • (3.4) 

□ 

For no receiver guard zone (e = 0) and a path loss exponent of a = 4 
{6 = |) we use Cor. l2.7l to express the OP m a more shnple and exphcit 
form. 

Corollary 3.1. Explicit OP for 5 = i. For e = and (5 = i the OP 
in Prop. 13.11 is 



q{\) = 2F, l^^ ^^\ c,Xj - 1. (3.5) 
where -Fz is the standard normal A^(0, 1) CDF. 

Proof. Write Z = Z^''j'(o) for this proof. For the assumed 5 = 1/2 we 
have from Cor. 12.71 that Z is a Levy RV (Def. 12. 9p with j = 2tt and 

CCDF F^{x) = 1 - 2Fz (\f^) = 2Fz - 1. Using Prop. lO 



gives 

'2t-2d 



q{\) = ( {\c,/2r'Cn = 2Fz ( WtT^^I^CT " 1- (3-6) 



(Ac,/2)-2e- 



□ 



For e = we use Prop. 13.11 and Def. 11.31 to express the TC in terms of 
the inverse of the CCDF of the RV zf f (o). 
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Proposition 3.2. TC (e = 0). For e = the TC equals 



5 



M^-f) (1-9*) 

Cd[F^ {q*)r 

where F''^ is the inverse CCDF of the RV ^l^/'°(o). 



Proof. Write Y. = Y.\^^''^{o) for this proof. Let F-£{y) be its CCDF and 
F~'^{q) the inverse CCDF. Equate the OP q{X) with the target OP q*: 

q* = Ft ((Ac,/2)-i/^r") ^ F^\ql = {Xc,/2)-'/' C'' ■ (3.8) 
Now solve for A: 

2 

^ ^ cd'{F^Hq*)y ^^-^^ 

Finally, multiply by 1 — g* as in Def. 11.31 □ 

The condition e = is necessary in order to decouple the RV T.yf'^{o) 
from the parameter A, thereby enabling solving for A. An analogous 
result for the OP may be derived for e > in terms of the inverse 
CCDF of the RV Z|^'^(o), but this expression cannot be solved for A, 
and therefore there is no explicit expression for the TC. For e = and 
S = ^ we can use Cor. l3.1l to get a more explicit expression for the TC. 



Corollary 3.2. TC (e = and 5 = i). For e = and 5 = i the TC 
in Prop. 13.21 is 

where is the inverse of the normal A^(0, 1) CDF F^. 



The corollary is immediate from Cor. 13. H and Prop. [321 The expressions 
for exact TC and OP for 5 = i {d = 2 and a = 4) in Cor. OandEl 
are shown in Fig. 13.11 for varying SINR thresholds r, with u = 1, N = 
0,P = 1. Additional exact OP and TC expressions will be given for 
the case of Rayleigh fading in §4.11 
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if 




3.2 Asymptotic OP and TC 



A(q*) 
I.IO 

— T-- 




I (d = 2 and 



Fig. 3.1 The OP q(X) vs. A (left) and the TC A((?*) vs. q* (right) for S = 
a = 4) and SINR threshold r e {1, 5, 10}. 



In this section we obtain the asymptotic OP in the hmit as A — t- and 
the asymptotic TC in the limit as g* — by applying Cor. l2.2l to Prop. 
13.11 and Prop. \3.2\ both valid for the special case of e = 0. 

Proposition 3.3. Asymptotic OP and TC. For e = the asymp- 
totic OP as A — )• is: 

9W = , "'^ ,s + ^ ^ 0- (3-11) 

For e = the asymptotic TC as g* — )• is: 

Hql = TJ Q*+ 0{q*)\ q* ^ 0. (3.12) 



Thus the OP is linear in A for small A and the TC is linear in q* 
for small q* . These asymptotic approximations will be used in ^4.21 
on variable link distances. The following remark gives the asymptotic 
approximations for OP and TC for the special case 5 = 1/2 (and = 0) 
in Cor. lOl 

Remark 3.1. TC as sphere packing. The first order Taylor series 
expansion of (1 - q*)F-^{{l + q*)/2) in (fXTOD around = is 

(1 - q*)F-' (^) = ^q* + 0{q*f. (3.13) 
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Using this in Cor. 13.21 for no thermal noise {N = 0) and rearranging 
gives a low OP approximation for the TC for 5 = 1/2: 

A(,*) = + ^ 0, uiq*) = n ( jj) ' • (3.14) 

In particular, (j3.14p may be interpreted as the number of d-dim. spheres 
per unit area, each with radius u{q*). Observe that u{q*)/u is the guard 
zone factor by which each Tx-Rx distance u must be expanded to ac- 
count for the required SINR threshold r, the required outage proba- 
bility q* , and the stability exponent 5. Also note the asymptotically 
tight UB on TC in Prop. 13.41 has the same expansion (for e = 
and iV = 0) as Prop. 13.31 In this case the series expansion being 
(1 — g*)log(l — q*) = q* + 0{q*)'^ as q* — >• 0. These expansions are 
seen to be accurate over a reasonable range of q* in Fig. 13.21 




Fig. 3.2 The first order Taylor series expansion {^Tr/2q*) of (1 — q*)F^^{{l + g*)/2) for 
Fz the standard normal CDF used in the exact TC for S = 1/2 in Cor. 13.21 (left), and the 
expansion (5*) of -(1 - i}*)log(l - q") used in the TC UB (for 5 e (0,1)) in Prop. IXil 
(right). 

3.3 Upper bound on TC and lower bound on OP 

In this section we obtain an UB on TC (LB on OP). The bound is 
based on considering only "dominant" interferers and interference. 



Definition 3.3. Dominant interferers and interference. An in- 

terferer i £ Ild,x is dominant at o under threshold r if its interference 
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contribution is sufficiently strong to cause an outage for the reference 
Rx at o: 

W(|xil) >r"^e< |xi| <e. (3.15) 
Else i is non-dominant. The set of dominant and non-dominant inter- 
ferers at o under r is 

t^dAo) = {ie Ud,x : e < |x»| < e} , f^dAo) = Ud,x \ tldAo). (3.16) 
The dominant and non-dominant interference at o under r 

«Gnd,A(o) «end,A{o) 
are the interference generated by the dominant and non-dominant 
nodes. Note Z;j'^(o) = Z^'^(o) + t^Ko). 

The LB on OP is obtained by observing the aggregate interference 
exceeds the dominant interference, and thus the probability of the ag- 
gregate interference exceeding some value exceeds the probability of 
the dominant interference exceeding that value. 

Proposition 3.4. OP LB and TC UB. The OP has a LB 

gib(A) = l-e-^"''(«'-^'). (3.18) 

The TC has an UB 

A-N,-) = (3.19) 

When e = 0, the bounds are tight for A — )• and q* — )• 0, respectively. 

Proof. The key observation is the equivalence of the events {£^'^(0) > 
and {Ild^\{o) / 0}, where we observe nodes in the annulus ad{o, e, ^) 
are dominant interferers. Prom here we compute the corresponding void 
probability for Ild,\ using Prop. 12.61 

q'\X) = P (z^;^(o) > r") 

= p(nrf,A(o)/0) =i-p(nd,A(o) = 0) 

= l-FiUd,xMo,e,0)=0) 

= i_ 6-^-^(5'*-'^). (3.20) 
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Set this last equation equal to q* , solve for A, and multiply by 1 — 
as in Def. Oto get the TC UB. 

Finally, the tightness of the bounds can be verified by comparison 
with Prop. [3^ with Taylor series expansions of e"^'^''^^ ^ and log(l — 
q*) around A and q* . □ 



Remark 3.2. Dominant and maximum interferers. The LB on 

OP obtained via dominant interferers is exactly equivalent to the LB 
on OP obtained by retaining only the largest interferer: 

<J4> min > <^ 

«end,A:|xi|>e 

<^ niax |xi|~"l| |>, < ^ 

^ M^;^(o)<e, (3.21) 

for M^'^(o) in (12:211) and Def. [231 In our extensions of the OP LB 
(c./. Def. ia Prop. M\ Lem. [521 Prop. EM Def. \5M and Prop. [5T2]) 

we won't make this correspondence with the largest interferer explicit, 
although adapting the above derivation to establish this relationship 
in those cases is straightforward. 



Specializing Prop. [33] to the case e = 0, = and S = and compar- 
ing with Cor. 13.11 and 13.21 (with = 0) gives the following corollary. 



Corollary 3.3. OP and TC bounds {e = 0, N = 0, 6 = The 

LB on the OP and the exact OP are: 

qib(^) = I _ Q-Xc^u^V^ < ^ (y^Acdu'^^/f) - 1. (3.22) 
The UB on the TC and the exact TC are: 



(3.23) 
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These expressions are plotted in Fig. 13.31 for d = 2 and a = 4 (with 
u = 1 and r = 5). Note the OP LB appears to be asymptotically 
exact as A — )• 0, and the TC UB appears to be asymptotically exact as 
— )• 0. The expressions in (j3.22p may be simplified to 

F,{z) > F^\z) = 1 - ^e-^^ z > 0, (3.24) 

by defining z = Xcdu'^^/r. This inequality is shown in Fig. 13. 4( again 
note the bound appears to be asymptotically exact as z — )■ 0. 




3.4 Throughput (TP) and TC 

One justification for claiming that TC is a natural performance metric 
is obtained by comparison with MAC layer TP, as discussed below. 



48 Basic model 



Definition 3.4. The MAC layer TP of a wireless network employing 
the slotted Aloha MAC protocol, where the active transmitters form a 
PPP Ud,x, is 

A(A) = A(l-g(A)), (3.25) 
where q{X) is the OP in Def. 11.11 

The TP has units of successful transmissions per unit area, and (|3.25p 
may be read as saying the TP is the intensity of attempted transmis- 
sions per unit area (A) thinned by the success probability (1 — q{X)) of 
each transmission. In light of Rem. 11.11 the design question for the TP 
is: given Apot how to select ptx so as to maximize A(ptxApot)- Before 
considering this question, we first recall the following basic facts about 
(saturated) slotted Aloha in a wireless uplink setting under the colli- 
sion channel model with a single collision domain (no spatial reuse). 



Proposition 3.5. Slotted Aloha TP and OP. For slotted Aloha 
in a single collision domain under the collision channel model with N 
saturated (backlogged) users transmitting to a common base station, 
employing a common transmission probability p, the TP and OP are 

A{N,p) = Np{l - p)^'^, q{N,p) = 1 - (1 - p)^'\ (3.26) 

The TP optimal p is p* (N) = 1 /N with associated TP and OP 

A(N,p*{N))=(^l-^y \ q(N,p*{N)) = 1-(^1-^ 

(3.27) 

For large N and p = the asymptotic TP and OP are 
A(A) = lim A(iV,A/A^) = Ae-^ 

q{X) = lim q{N,X/N) = l-e-^. (3.28) 
Af— s>oo 

The asymptotic TP optimal choice for A is A* = 1 with associated 
asymptotic TP and OP 

A(A*) = le-i f« 0.367879, q{X*) = 1 - e"^ w 0.632121. (3.29) 
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The statements in Prop. 13.51 are simple to prove. Note that achieving 
the maximum asyptotic TP of 37% requires an incurred OP of 63%. 
These relationships are illustrated in Fig. 13.51 We now return to the 
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Fig. 3.5 The TP (left) and OP (right) of slotted Aloha in a single collision domain (no 
spatial reuse) for N backlogged users employing transmission probability p = A/Af, for 
N S {5,10} Also shown is the TP and OP as Af — ^ oo. Achieving a high TP requires 
incurring a high OP. 

context of wireless networks with spatial reuse. The LB on the OP in 
Prop. [33] leads directly to an UB on the TP in Def. [331 This UB is 
equivalent to the TP in the non-spatial context of Prop. 13.51 

Proposition 3.6. MAC layer TP UB. The TP in Def. [331 has UB 

A(A) < A"b(A) = A(l - gi^(A)) = Ae-^^'^(«'-^'). 
The TP bound optimal A is 

A* 



1 



and the associated TP UB and OP LB are 
1 



(3.30) 
(3.31) 



A"^(A*) 



, gib(A*) 



ecd(e'^ 



1 



0.632121. 



(3.32) 



The main point of Prop. 13.61 is that achieving the optimal TP requires 
incurring an OP of 63%. Given that wireless devices are energy con- 
strained and that failed attempted transmissions are wasted energy, it 
is natural to question if unconstrained TP maximization is the right 
design objective. TC is in fact the TP of a wireless network under Aloha 
subject to an OP constraint, as shown below. 
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Proposition 3.7. TC is constrained TP maximization. The op- 
timization problem of maximizing TP subject to a constraint on the 
outage probabihty q* G (0, 1) 

max{A(A) : q{X) < q*} (3.33) 
A 

has solution A* = q~^{q*) and TP equal to the TC in Def. Ol 

A(A*) = A*(l - g(A*)) = q-\q*)il " Q*) = A*(g*). (3.34) 



To clarify, A* = q'^{q*) is the solution of (iSlSSjl . and X*{q*) is the TC as 
defined in Def. ll.3l In summary, TC is TP under an OP constraint. Note 
that maximization of TP A(A) over A G M is equivalent to maximization 
of the TC A*(g*) over the target OP q* G [0, 1], as made precise below. 



Proposition 3.8. Maximum TP equals maximum TC. The TP 

optimization problem 

max A(A) (3.35) 



has a unique maximizer Aopt and an associated maximum value Aopt = 
A(Aopt)- The TC optimization problem 

max \*(q*) (3.36) 
<?*e[o,i] 

has a unique maximizer q*^^ and an associated maximum value A*pt = 
X*{q*p^). Furthermore, the maximum values are equal and the maxi- 
mizers are related through the OP function q{X): 

^opt = A*pt and ^(Aopt) = ^opt- (3-37) 



Proof. We first prove (/(Aopt) = ^opt- The monotonicity of q{X) guaran- 
tees that A(A) has a unique stationary point on (0, oo) and that this is 
the unique maximizer Aopt > 0. Taking the derivative 



A'(A) = l-g(A)-Ag'(A) 



(3.38) 
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and equating with zero at the stationary point gives 

A'(Aopt) = 1 - g(Aopt) - Aoptg'(Aopt) = 0, (3.39) 
or equivalently, after rearranging, 

^'(Aopt) = 1^1^^. (3.40) 

It is hkewise straightforward to show X*{q*) is concave on [0,1] and 
therefore the unique optimizer is the unique stationary point on (0, 1). 
Taking the derivative and applying the inverse function theorem gives: 

X*'{q*) = (1 - q*)q-''{ql " g-^?*) = ^7(^^^ " 1~Hq*)- (3-41) 



Equating with zero at the stationary point gives 
or equivalently, after rearranging. 



Q'iQ-HQ*)) = ^ZTT^y (3.43) 

There is a unique solution Aopt for (I3.40p , and likewise there is a unique 
solution g*pt for ()3.43p . The two optimizers are related by g(Aopt) = Qopt 
(equivalently, Aopt = Q'^ilopt)) since under this relationship these two 
equations are the same (in that Aopt solves (|3.40p iff q*^^ solves (|3.43p ) . 
We next prove Amax = Aj^^x- Rearranging (I3.40p and (I3.43P gives 



Amax — Aopt <?'( Aopt), A^ax — ,/ _i/ * • (3.44) 

y W {Qopt)) 

The square root of their ratio is 

-g'(Aopt) = l, (3.45) 



' Amax Aopt 



V Amax 1 ~ Q'(Aopt) 
and thus Amax = Kva.yL- ^ 



The corollary below shows that Prop. [3^ holds for the TP UB in Prop. 
Mland the TC UB in Prop. [M 
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Corollary 3.4. Maximizing TP and TC UBs. Denote Cd(^'^ - e'') 
found in Prop. [T^ bv a. The TP UB optimization problem 



max A"^(A) 



(3.46) 



for A"'^(A) in ()3.30p has a unique maximizer Aopt = 1/a and an associ- 
ated maximum value AJ^^ = m- The TC UB optimization problem 



max A"^(g*) 

9*e[o,i] 



(3.47) 



1 — i and an asso- 



for A"^(g*) in (l3T9]) has a unique maximizer Qop^ 
ciated maximum value X"^^^ = ^. Like Prop.EHl the maximum values 
are equal and the optimizers obey q*^^ = (?'^(Aopt) for q^^ in (j3.18p . 

The corollary is obtained by simple calculus on the functions A"'^(A) 
and A*^*^ ((?*). Fig. 13.61 illustrates the quantities in Cor. 13.41 

A(A) 
0.05 




Fig. 3.6 The spatial TP UB A"'^(A) ITSOt vs. the spatial intensity of attempted trans- 
missions A (left), and the TC UB A"''((j*) vs. the target OP q* (right). For a = 4,d = 
2,u = 1,T = 5, N = we have a = c^^'^ = ^/Eir 7.02 and optimal A* = 1/a 0.14 and 



maximum TP An 



l/(ca) Si 0.052, g^pt = 1 - 1/c 0.63, and g*^^ = g"^(Aopt). 



3.5 Lower bounds on TC and upper bounds on OP 

In this section we obtain three LBs on TC (UBs on OP). The bounds 
are obtained by using the bound from §3.31 along with three UBs on 
the tail probability for the non-dominant interference. The tail UBs are 
given in ^2.1l and these in turn employ the moments from the Campbell- 
Mecke theorem given in Thm. 12.31 We express the OP in terms of the 
distributions of the dominant and non-dominant interference. 
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Proposition 3.9. Exact OP in terms of OP LB. Define the 

CCDFs of tlie independent RVs t^'^(o), t'^'^io) from Def. [331as: 

F^iy) = P(r^;^(o) > y), F^{y) = P(Z^;^(o) > y). (3.48) 
The OP equals: 

q{x) = + i^£(r")i^f (r"). (3.49) 

Moreover, the OP equals: 

q{X) = q'\X) + (1 - q'\X))F^iCn (3.50) 
for q^^{X) given in Prop. 13.41 



Proof. Write 51 = 5I^'^(o),II = Y.'^'^{o),Y. = Z^'^(o) for this proof. 

Recall Z = Z + H. Express the outage event for Y. in terms of (Z, Z) 
and decompose it into three (overlapping) regions: 

{Z>r"} = {(Z,Z):Z + Z>r"} (3.51) 

= {z > r"} u {z > r"} u {z < ^ < r", ^ > r"}- 

Recall the equivalence of the events {Z < ^^"} and {Z = 0}, and 
observe this implies the third event becomes null: {Z < n Z > 
^-"} = 0. The OP is therefore 

q{x) = p(z > c") + ip(^ > r") - p(^ > r" n z > r")- (3.52) 

Apply the independence of (Z,Z) to the third term and group terms 
to get ([339]). Finally, recognize g'^(A) = F^(C"°) to get (f330]l . □ 

In ([T^ it is clear that we can UB q{X) by an UB on the CCDF of 
the non-dominant interference -F^(^~"). We now apply the Markov, 
Chebychev, and Chernoff bounds to the RV Z|^'^(o). 



Proposition 3.10. Markov inequality OP UB for a > d is 

^ub,Mar(^) = q'^{X) + (1 - q'^iX))^^^''. (3.53) 
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Proof. Markov inequality (Prop. [T3|) for Z|J'^(o) in Def. 13.31 gives 

p(z^;^(o)>r")<riEK:I(o)]. (3.54) 

Now apply Campbell's Theorem (Tlim. [2^ to compute E[Z^'^(o)]: 



E 



oo 



Xdcd I r-V-Mr = ^^r^+'i. (3.55) 
^ a — a 



The integral is finite for a > (i. Substitution yields the proposition. □ 
Second, apply the Chebychev inequality. 

Proposition 3.11. Chebychev inequality OP UB for a > d is 

^.b,Cheb(^) = q^\X) + (1 - ^^^A))-^- (3.56) 

The bound is trivial for 

a — d , . 



Proof. Denote Y. = Y.'^'^(o) for this proof. The Chebychev inequality 
(Prop. EID applied to the RV t in Def. Ogives 

P(Z>r") = P(^-E[Z] >r°-Ep]) 

< P(|Z-E[Z]| > -Ep]) 

Var(Z) 

< — . (3.58) 

Now apply Campbell's Theorem (Thm. [2^31) to compute Var(Z): 



Var \xi\ "l|x,|>5 = Adcrf / 



2a -d 



(3.59) 

Evaluating the integral for a > d/2, substituting E[Z] = ^^C"""^"*, 
and cancelling the common term yields the proposition. The first 
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inequality used above yields a trivial bound for ^ " < which may 
be expressed as bound on A: 

r" < Ep] ^ e"" < ^^^-"+'^. (3.60) 

a — a 

□ 

The threshold ()3.57p is sufficient but not necessary for the bound to be 
trivial. Third, apply the Chernoff inequality. 

Proposition 3.12. ChernofF inequality OP UB ioi a > d is 

^ub.Cher^^-) = q^^{X) + (1 - g"^(A))e-"(^) , (3.61) 

where 



c(A) = sup I (e'y- l) y-'-'dy ] . (3.62) 

6»>0 \ CH Jo 



Proof. Denote Y. = Z^'^(o) for this proof. The Chernoff inequality 
(Prop. [23]) applied to the RV t in Def. gives 

P(Z > r°) < inf E[e^^]e-^«"". (3.63) 

e>o 

The MGF of 51 is obtained by selecting e = ^ in Cor. [23] yielding: 

^ a Jo 

Simple manipulations yield the proposition. □ 



M[t]{e)=e^pl^l^^ \e<^y-l)y-'-'dy\. (3.64) 



Finding the optimal 9* in (I3.62p must in general be done numerically, 
although certain simplifications hold for 5 = 1/2. 

The Markov, Chebychev, and Chernoff bounds on the OP and TC 
are shown in Fig. 13.71 The top two plots are OP vs. A and the bottom 
two plots are TC vs. q* . The right plots are an inset of the left plots. 
Consider the OP plots. The OP LB is seen to be tighter than each 
of the three OP UBs. For A small (corresponding to small q{X)) the 
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Chernoff and Chebychev bounds are nearly equivalent and are better 
than the Markov bound. For moderate to large A (corresponding to 
moderate to large q{X)) the Markov and ChernofF are nearly equivalent 
and better than the (trivial) Chebychev bound. All three bounds thus 
have their value: i) Markov is bad for small A but good for larger A and 
is simple, ii) Chebychev is good for small A but bad for larger A and 
is intermediate in simplicity between Markov and Chernoff, and in) 
ChernofF is as good as Markov and Chebychev for all A, but is much 
more complicated than the other two. Similar trends naturally follow 
for the three TC plots as for the OP plots. 
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Fig. 3.7 Top: the OP q(A), the LB, and the three UBs (Markov, Chebychev, Chernoff) vs. 
A e [0,0.1] (left) and A e [0,0.01] (right). Bottom: the TC A(g*), the UB, and the three 
LBs (Markov, Chebychev, Chernoff) vs. q* e [0,1] (left) and q* e [0,0.1] (right). 
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Extensions to the basic model 



In this chapter we study three extensions of the basic model in Ch. O 

(1) Channel fading: allow for iid channel fading on top of the 
pathloss attenuation, with Rayleigh fading a special focus. 

(2) Variable link distances: allow each Tx~Rx pair to be sepa- 
rated by a random distance u, iid across pairs. 

(3) Multi-hop networks: measure performance for a multi-hop 
extension of the model using suitably modified OP and TC 
metrics. 

These three extensions are chosen because they are among the most 
obvious steps towards a more realistic decentralized network model. We 
shall see that fading can be added to the model without any major dif- 
ficulties, and in fact when all fading is Rayleigh, exact results are easier 
to compute than without fading. Variable link distances are straight- 
forward to include, and result only in a multiplicative constant, which 
justifies the use of the less realistic but simpler fixed distance model 
adopted in the rest of the monograph (and most of the literature) . Mul- 
tihop is a nontrivial extension and requires end-to-end definitions of OP 
and TC, but under a simple model we are able to preserve tractability 
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and determine quantities like the optimum number of hops. 
4.1 Channel fading 

In this section we modify the basic model of Ch. [3] to include channel 
fading. Recall in Ch. [2] and [3] we operated under Ass. 12.21 where in 
particular the amplitudes {hj} in Def. 12.41 were assumed to be unity. 
We now relax that assumption. 

Definition 4.1. SINR under fading. Fix e = in this section. Let 
ho be the fading coefficient on the signal channel between the reference 
Tx and the reference Rx at o, and let {hj} = hi, h2, . . . be iid RVs rep- 
resenting the fading coefficients on the channels between each interferer 
and the reference Rx at o (as in Def. 12. 4p . Let S(o),5I(o) denote the 
random signal and interference powers seen at o each normalized by 
the transmission power P: 

S(o) = hon-", and Z(o) = Z;^;j;(o) = Yl ^^l^^l""" (^'l) 

The SINR at o is as in the basic model in Def. 13.11 with S(o) and 51 (o) 
updated as above: 

Note we have normalized the noise by P since we defined S, 51 to be 
normalized by P as well. 



Remark 4.1. Signal and interference fading coefficients. Under 
Ass. 14.1 1 both the received signal S(o) and interference power 5I(o) are 
random, where randomness in S(o) is due to hg, and randomness in 51 (o) 
is due to both the random positions in II^^a and the fading coefficients 
{hj}. Note the convention that {hj} = hi, h2, . . ., and in particular the 
signal channel fade ho is not in the collection of interference channel 
fades {hi}. To be clear, h denotes a generic interferer fading coefficient, 
hj denotes the fading coefficient for interferer i, and ho denotes the 
signal fading coefficient. Throughout this section it is important to 
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observe the distinct impacts of h vs. ho on performance, and the distinct 
requirements for analytic tractability. 

This section is divided into three subsection: exact OP and TC ( §4.1.1|) . 
asymptotic OP (as A — 0) and TC (as q* — )• 0) ( §4.1.2p . and a lower 
(upper) bound on OP (TC) ( UTT^ . 

4.1.1 Exact OP and TC with fading 

The Laplace transfom of the interference Z^'^(o) is given in the follow- 
ing proposition ([36] (3.20)). 

Proposition 4.1. LT of the interference. The Laplace transform 
of the interference (o) under Def. 14.11 and for 5 < 1 is 

C[T2^^{o)]{s) = exp {-AQE[h^]r(l - 6)s'] , seC. (4.3) 

The RV Z^']^(o) is stable: the characteristic function (f>[T.'^'^{o)]{t) is 
given by Def. 12.81 where the characteristic exponent is 5 < 1 and the 
dispersion coefficient is 

7 = (AcdE[h^]r(l - 6) cos(7rV2)) * • (4.4) 

For the special case 6 = 1/2 the RV Z^']^(o) is Levy as defined in Def. 
12.91 with parameter 

7 = |(Ac,E[A])'. (4.5) 



Proof. The proof of (|4.3|) is given in [36] (3.20). Write Y. = Z^'t'(o) for 



this proof. We find the CF for Z from (j4.3p . Using the expression for 
7 in (14. 4p we write the Laplace transform as 

In general we can obtain the CF from the LT via (j2.8p . 
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To put this in the form of Def. 12.81 we must estabhsh: 

' it)'' = -7'^|t|'^(l -itan(7r5/2)sign(t)). (4.8) 



cos(7r(5/2) 

We consider the case t > 0, the case t < is similar: 

-i = e-i-/2 

(-i)<5 = cos(7r(5/2) -isin(7r(5/2) 

= l-itan(7r5/2). (4.9) 

cos(7ro/z) 

□ 

We obtain an exphcit expression for the OP and TC when the signal 
fade ho is assumed to be exponentially distributed (usually interpreted 
as modeling Rayleigh fading). The following result is found in [6] and 
is discussed in |Ml ^3.3. 



Proposition 4.2. OP and TC under Rayleigh signal fading, 
general interference fading ([6j). The OP and TC under Def. 14.11 
for signal fading coefficient ho ~ Exp(l) (Rayleigh signal fading) are: 



/(A) = 1 - exp I -AcdE[h^]r(l - 6)t^u' 



5„.d '_ 

snr . 



(-log(l-g*)--^)(l-g*) 



Proof. Write 51 = 51^'^ (o) for this proof. Solve the outage event for ho, 
condition on Z, and apply the exponential CCDF: 

1-qiX) = P(sinr(o)>r) 

= P(ho > ™°(Z + iV/P)) 

= E[P(ho > Tn"(Z + iV/P) I Z)] 

= E[exp{-r'u"(Z + iV/P)}] 



(4.11) 
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Now recognize the Laplace transform £[5I](s) at s = tu'^, use the defi- 
nition of Rx SNR in Ass. 13. H and apply Prop. 14.11 

l-g(A) = e-s5?£p](s) 

= e-i^exp|-AcdE[h'5]r(l . (4.12) 

L J S=TU°' 

Simplification gives the proposition. □ 

The following corollary gives the OP and TC for the special case when 
the interference coefficients are also unit exponentials (i.e., Rayleigh 
fading for both signal and interference channels). The key step to obtain 
the corollary is the following expression for the fractional order moment 
of the exponential distribution. 

Lemma 4.1. Moments of exponential RV. For a unit rate expo- 
nential RV h ~ Exp(l) and 6 G (0, 1): 

E[h^] = r(l + 5), E[h-^] = r(l-(5). (4.13) 



Lem. 14.11 and the identity (|2.7p yield the following corollary, which is 
one of the most widely used results appearing in this monograph, and 
due to 0. 

Corollary 4.1. The OP and TC under Rayleigh fading 

(ho,hi,h2,... ~Exp(l)) are: 

q[X) = 1 - exp <^ -A . r u 

1^ sm(7roj snr 



sin(7r5) (-lo g(l-g--)- ^) (1 - q*) 



K<1*) = " l.^y. (4.14) 



Notice that this is an exact closed-form result for OP and TC that does 
not require bounds, asymptotics, or approximations. We can specialize 
it even further by fixing 6=1/2 and = (snr = oo) in the above 
corollary. 
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Corollary 4.2. The OP and TC under Rayleigh fading {6 = ^, 

N = 0) are 



^ -2Mi-p(i-,-) 



We compare the impact of fading on the OP and TC in Fig. 14.11 by 
plotting tlie OP and TC with Rayleigh fading for S = 1/2 (from Cor. 
14. 2p with the non-fading case for 5 = 1/2 (from Cor. I3.1l and l3.2( ). The 
plots illustrate that fading degrades performance in the low outage 
probability regime. This may be roughly understood by arguing that 
fading is a source of variability (and hence diversity) but our slotted 
Aloha MAC protocol does not exploit that diversity, and therefore per- 
formance suffers under the variability. This degradation is quantified 
in the asymptotic (A — >• and q* — >• 0) regime in Cor. 14.31 and Fig. 14.21 
We will see how this diversity may be exploited via scheduling ( §5.3p 
and power control ( ^5.4p leading to OP and TC that outperform the 
non- fading counterpart. 




Fig. 4.1 The OP (left) and TC (right) for both no fading and Rayleigh fading at 5 = 1/2 
and N = 0. The other parameters are u = I, t = 5, and d = 2 (a = 4). 



4.1.2 Asymptotic OP and TC with fading 

The asymptotic OP and TC in the basic model of Ch. [3] were given 
in ^3.21 (Prop. [3T3j) . These results were obtained by i) mapping 5I|^'^(o) 
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to (o) using Prop. 12. 8( ii) applying the series expansions of the 

CCDF of ^I'x in Prop-ETTUl and in) using this in Prop . \3A] which gives 

the OP in terms of the CCDF of Y.Y^'^{o) evaluated at a certain y. We 
now extend this same procedure to the fading model in Def. 14.11 

Remark 4.2. Marked PPP (MPPP). The appropriate formalism 
to incorporate channel fading under Def. 14. II is that of the marked Pois- 
son point process (MPPP), denoted by = {(xj; hj)} where the RV 
hj G is the "mark" associated with point x, G M*^. In general the dis- 
tribution of the mark is allowed to depend upon the point, but marks 
must otherwise be independent of other marks and other points. Let 
hj E M_|_ be the fading coefficient for the channel between the interferer 
at Xj and the reference Rx at o, and assume that hj is independent of 
Xj and the other points and marks. The key result is that the MPPP 
^d,\ with pairs (xj, hj) taking values in x is in fact a (in general 
non-homogeneous) PPP with "points" (xj, hj) taking values in R'^ xR_|_. 
Thus the addition of marks does not spoil the tractability of the un- 
marked PPP framework. See |49] Ch. 5, and in particular the marking 
theorem in §5.2, given below. 



Theorem 4.1. PPP marking theorem (09]). The MPPP $d,A = 
{(xj,hj)}, where {xj} is a (in general, non-homogeneous) PPP on 
of intensity A(x) and the marks admit a conditional PDF f^y^{h\x) for 
each X G M"^, is a (non-homogeneous) PPP on x with intensity 
measure 

KC) = [ A(x)A|,(/i|x)dxd/i, (4.16) 
J{x,h)ec 

for all Borel sets C C M'^' x M+. 

Prop. 12.61 gives the void probability for a homogeneous PPP; this is 
generalized for non-homogeneous MPPPs below. 



Proposition 4.3. The void probability of the non-homogeneous 
MPPP ^d,x with intensity measure ^(C) in Thm. 14.11 is 

P($rf,A(C) = 0) = e-^(^), (4.17) 
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for all Borel sets C C M'^ x E. 



Thm. allows for marked versions of Prop. [2771 and [THl 

Proposition 4.4. MPPP distance and interference mapping. 

Let ^d,x = {(xj^hj)} be a MPPP where {xj} is a homogeneous PPP 
on of intensity A and the marks {hj} are independent of the points, 
and let $i^i = {(tj, hj)} be a MPPP where {u} is a homogeneous PPP 
in M of intensity 1. Then: 

{Xcd\xi\'^,hi) = {2\ti\,hi), ieN. (4.18) 

Further, the following RVs are equal in distribution: 

^s;(")^(^)'^5^"). (4.19) 



Thus, as with the non-fading case, the interference may be treated as 
a scaled version of that generated by a unit rate PPP on M. Mapping 
from M'^ to M. is important because it allows us to apply a series repre- 
sentation of the CCDF (as in Prop. I2.10p and take the dominant term 
to get the asymptotic CCDF as y — t- oo (as in Corr. 12. 2|) . Prop. 12.101 
(taken from [57j (29)) ignored the marks since in that section we as- 
sumed hj = 1 (Ass. 12. 2p . but in fact ^57j gives the series expansion for 
given below. Recall the footnote under Prop. \27l0[ 



Proposition 4.5. SN series expansion ([57]). The series expansions 
of the PDF and CCDF of the SN RV T."''? for 5 = ^ < 1 are: 

f 4y) = — V ^ \ r(l + n5) sin(7rn^)(2Ar(l - J)E[h^]y-^)" 
■^i,A Try ^ — ' n! 

^ n=l 

F ,(y) = — ^r(l + n(5)sin(7rn5)(2Ar(l -5)E[h'5]y-^)'^ 

^i,A Trd nn! 

71=1 

(4.20) 
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The asymptotic PDF and CCDF as y — )• oo are: 

= 2XSE[h%-'-' + 0{y-^-^'), y^^ 
F^.^y) = 2XE[h%-' + Oiy'^'), y^oo (4.21) 

The following theorem establishes the stability (in the sense of Def. 
EH]) of T.l^J'^{o) (adapted from [6l] Thm. 1.4.5 and 02] Thm. 3). 

Theorem 4.2. Interference with fading is stable ([6i]). For 6 < 1, 
liy^''"{o) is stable (with CF in Def. 12. Sp . with characteristic exponent 
5 and dispersion coefficient 

7 = (r^r{2 - 6) cos(7r5/2)E[h^] j ' . (4.22) 



Remark 4.3. Interference and fading moments. In both Prop. 
14.51 and Thm. IT2] the impact of the random marks (fading coefficients) 
is restricted to the fractional order moment E[h'^]. See [36] (p. 33) for 
an extended discussion of this fact. 



The asymptotic CCDF of Z"']^ in Prop. [33] leads directly to the asymp- 
totic OP (as A 0) and TC (as q* -?- 0) under fading, just as Cor. [Q 
led directly to Prop. 13.31 A key difference is that random signal fading 
invalidates Ass. 13.11 

Remark 4.4. Fading and outage with no interference. With ran- 
dom signal fading ho there is the possibility of a bad fade causing outage 
even in the absence of interference. This outage event 

for snr defined in Ass. 13.11 has probability 

g(0) = P(£:o)=Fh„(^). (4.24) 
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Note q{0) is the OP evaluated at A = 0. Denote the complement of Sq 
by £o, and its probability by 

m = 1 - q{0) = F{£o) = Fho (— ) • (4.25) 

All analysis of OP (and hence TC) must therefore condition on hg 
being above or below r/snr to distinguish between the case of outage 
being possible vs. outage being guaranteed. The range of the OP q{X) 
is [g(0), 1], and the domain of q* in the TC X{q*) is [q{0), 1]. 



Proposition 4.6. The asymptotic OP under fading as A — t- is: 



liX) = 1 - 1 - CdE[h^]E 



N 



TU'^ P 



£o 



Aj m + o{x^) 

(4.26) 



The asymptotic TC under fading as q* — t- g(0) is 



1-9(0) 



CrfE[h^]E 



TU°' P 



So 



+ 0{q*-q{Q)f. (4.27) 



In the no noise case (A^ = 0, g(0) = 0) these expressions become 
g(A) = CrfrVE[h^]E[ho^]A + 0(A2) 

K<f) = T'f + Wf 

^ ' Qr'5n'^E[h5]E[hoT 



(4.28) 



In the no noise and Rayleigh fading case (Hq and {hj} exponential RVs): 

,d 



sin(7r(5) 



A + 0(A' 



sin(7r(5) ^ *n2 
jq* + 0{q*Y 



(4.29) 



Proof. The proof is analogous to that of Prop. [3T| with a key difference 
being the requirement to condition on possible vs. guaranteed outage, 
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c.f. Rem. 14.41 Applying Def. 14.11 to Def. II. 1^ and conditioning on hg 
leaves as the sole source of randomness in sinr(o): 



P(sinr(o) < r) 

P (sinr(o) < t\£o) ^(0) + 1(1 - g(0)) 
1 - (1 -P(sinr(o) < t\So)) g(0) 
1- (l-E[P(sinr(o) <T|ho)|^o])g(0) 
/ 

P 



1 



1 -E 







ho) 


So 



\ 



E[.: 



Note the probability is a RV that is a function of the RV ho, and the 
expectation is with respect to the conditional PDF for hg: 

f^,^{h\£o) = { ^w,(iS7)' (4.31) 
[ 0, else 

Now apply Prop. and 1^31 to the interference probability noting that 
ho and Z are independent. 

/ 



EM 



E 

E 
E 



^11 (o) > Hr 



V 

T' (o) > y 



So 



E[h^]y-^ + 0(y 



y > 
y > 



(4.32) 



Note So = {ho > ^} is equivalent to {y > 0}. Simplification yields 
Solving q{X) = q* for A yields dOZ])- Expressions are im- 

mediate upon substituting = 0, and expressions (|4. 29p are immediate 
upon applying Lem. 14.11 and ()2.7p . □ 



Remark 4.5. OP and fading moments. In Rem. 14.31 we noted 
the distribution of the interference depended upon the fading coeffi- 
cients only through E[h'^]. In the no noise case of Prop. 14.61 we see the 
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asymptotic OP and TC depend upon the signal fading coefficient ho 
and the interference fading coefficients {hj} only through the product 



Comparing the no-noise asymptotic OP and TC with fading in ()4.28p 
in Prop. 14.61 with the analogous results without fading in Prop. 13.31 
yields the following corollary. 



Corollary 4.3. Fading degrades performance. Let q{X),X{q*) de- 
note the OP and TC without fading, and q{X), A(g*) denote the OP and 
TC with fading. Let the signal and interference fading distributions be 
equal (h = hg). Under no noise (iV = 0) the asymptotic OP (as A — >■ 0) 
and TC (as q* — >■ 0) with and without fading have ratio 

|W=^=E|h'lE[h^'l>l. (433) 
For Rayleigh fading this ratio is given by Lem. 14.11 and (j2.7p : 

E[h^]E[ho ^] = r(l + 6)T{1 -6) = > 1. (4.34) 

sm(7roj 

Fading degrades asymptotic performance relative to non-fading. 



Proof. Jensen's inequality (Prop. [221) asserts E[/(x)] > /(E[x]) for con- 
vex /(•) and RV x. Use convex function f{x) = 1/x and RV h''. □ 

Fig. 112] shows vs. 6. For 5 = 1/2 {e.g., d = 2 and a = 4) this 

quantity is 7r/2 w 1.57, i.e., the asymptotic OP / TC is 57% worse in 
the presence of Rayleigh fading than without fading. 

4.1.3 LB on OP (UB on TC) with fading 

The OP LB (and TC UB) was established for the basic model through 
the use of dominant interferers, defined in Def. l3.3l as those nodes whose 
interference strength (under pathloss without fading) was sufficient to 
individually cause outage at the reference Rx. The following definition 
extends the concept to allow for fading. 




Definition 4.2. Dominant interferers and interference. An in- 
terfere! i in the MPPP ^d,x defined in Prop. 14.41 is dominant at o 
under threshold r and signal fade ho if its interference contribution is 
sufficiently strong to cause an outage for the reference Rx at o: 

< T 4^ h,|xr" > ^ - TT- (4.35) 



hi|xj|-" + iV/P ' ' ™" P 

Else i is non-dominant. The set of dominant and non-dominant inter- 
ferers at o under (r, ho) is 

^d,x{o) = |i G $d,A : h,|xi|"" > ^ - ^} ' ^d,x{o) = <^d,x \ ^d,x{o). 

(4.36) 

The dominant and non-dominant interference at o under (r, ho) 

^d'>) - E h.ixr", ^t(^) - E h.ixr- (4.37) 



are the interference generated by the dominant and non-dominant 
nodes. Note T.°'xio) = ^2a(o) + a(o). 



This definition leads directly to a LB on the OP and an (in general, 
numerically computed) UB on the TC. 
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Proposition 4.7. OP LB. Under Def-Othe OP has a LB 



glb(^) = 1 - E 



exp < —XcdMl 



ho N 



TU"^ P 



So 



m (4.38) 



where the outer expectation is w.r.t. the random signal fade hg. In the 
case of no noise (A^ = 0, q{0) = 0) the LB is: 



,1b 



(A) = 1 - E exp I 



exp<| -AcdT^n'^E[h>o^} 



(4.39) 



In particular, with no noise the LB is expressible in terms of the MGF 
of the RV — Hq ^ at a certain 6: 



(4.40) 



Proof. Repeat the proof of Prop. 14.61 up to (|4.30|) . LB in terms of 
o), and then express tl 

(A) = P(sinr(o) < r) 



^^(o), and then express the LB outage event in terms of ^d,x{o): 



= 1 
> 1 



1 

1 -E 



1 -E 

1 -E 
1 -E 



■r-a,h / \ ^ 



ho 


N 


rn" 


~ P 


ho 


N 




~ P 



IP ^d,x(o) = 



"0 



ho 

^0 



m 



ho 
ho 

m 



m 
m 



(4.41) 



The PPP ^d,x is a homogeneous PPP with intensity measure given by 
Thm. BTT] with A(x) = A and fh\x{h\x) = fh{h). The key observation 
is that the probability that ^d,x{o) = equals the void probability for 
^d,x on the set Co = {(x, h) : > wq} for wq = — ^. Note ho 

is random and hence so is Co and wq. Using Prop. 14.31 and simplifying 
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gives: 



■(^d,A(Co) = 0| he 



exp ^ —X I fh{h)dxdh > 

'{x,h)eCo J 



exp < —A 



'-fc|x|-">wo 



/h(/i)dxd/i 



= exp |-A^^E[lh|a.|-c>w„| wo]dx| 

= exp I P(h|xr" > wo| wo)dx|(4.42) 

Express the CCDF for h in terms of and apply Thm. l2.2l to reduce 
the integral over to an integral over R: 



exp < —A 



h 

Wo 



> \x\ 



Wo dx ^ (4.43) 



exp { -Xdcd I F\ [ — ] >r 



Exchange the order of integration as follows: 



exp |— Adcd 
exp < —Xdcd 



Wo r'^ ^dr 



'h > wor"| wo)r'^ Mr 



= exp Adcd 
Evaluate the inner integral and rearrange: 



Uh)dh\r''-^dr 

(/i/wo)!/" \ 

r'^'^drj fh{h)dh 



exp < — AdCfl 



1 / h 



k{h)dh 



expl-AcrfWo-^y h^Mh)dh^ 
expl-AcrfWo-'Efh^]} (4.44) 
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Substituting this last expression into (I4.4ip yields ()4.38p . □ 

Note the MGF Ai[—\r\Q^](6) is not in general expressible in closed form, 
and consequently the OP LB in Prop. 14.71 must be numerically com- 
puted, see Fig. 14.31 Further, the MGF is not in general analytically 
invertible, and thus the corresponding TC UB must be numerically 
computed. Fig. 14.41 shows several OP and TC expressions from ^4. II for 
the specific case of Rayleigh fading (for both signal and interference 
channels), d = 2 and a = 4 {5 = 1/2), and no noise {N = 0). The other 
parameters include r = 5 and u = 1. The plots include exact results 
from ^4.1.11 (Cor. 14. 2|) . asymptotic results from ^4.1.21 (Prop. 14. 6p . and 
bound results from ^4.1.21 (Prop. 14. 7p . 




Fig. 4.3 The MGF for the RV -h"^ in Prop. ET] (A4[-h-*] (6)) vs. 9, for h ~ Exp(l) and 
d = 2 and « e {2, 4, 6} (5 e {1, 1/2, 1/3}). 

4.2 Variable link distances (VLD) 

In this section we make perhaps the simplest enhancement to the basic 
model: the distance separating each Tx-Rx pair is allowed to be a RV 
u ~ -Fu iid across pairs, which we call variable link distances (VLD), as 
opposed to fixed link distances (FLD). For simplicity of exposition we 
fix e = throughout this section. We retain the assumption of constant 
power P for each node. There are now two sources of randomness: the 
PPP Ild,\ yielding a random interference 51 , and the Tx-Rx separation 
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Fig. 4.4 The OP (left) and TC (right) for Rayleigh fading (signal and interference), d = 2, 
a = 4 {5 = 1/2) T = 5, u = 1, and N = 0. Shown are the exact, asymptotic, and bound 
results from i|4.1l 



distance u, yielding a random signal power S = Pu requiring an 
updated definition of SINR. 



Definition 4.3. SINR for VLD. Under VLD u ~ F^, and e = 0, the 
SINR in Def. 13.11 holds but with signal power a RV S = Pu~". 



The total law of probability gives an updated definition for the OP. 

Definition 4.4. The OP for VLD with u ~ is 

qiX) = E[P(sinr(o) < t|u)]. (4.45) 



The impact of VLD on the OP and TC is most easily seen by applying 
the asymptotic OP and TC from Prop. 13.31 for the case of no noise 
{N = 0). Denote the OP and TC with VLD q{X), X{q*) and denote the 
OP and TC with FLD as q{X), X{q*). 



Proposition 4.8. Asymptotic OP (e = 0, N = 0). For e = and 

iV = the asymptotic OP with FLD n = E[u] and VLD u ~ are 

q{X) = CdT%ufX + O{X'^),X^0 

q(X) = CdT%u'^]X + 0{X'^), X^O (4.46) 
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For e = and iV = the asymptotic TC with FLD u = E[u] and VLD 
u ~ Fu are 

1 



q* + 0{q*f, q*^q{0) 



1 



q*+Oiq*f, q*^qiO) 



(4.47) 



Note for small A, q* the ratio of OP and TC under VLD vs. FLD is: 
g(A) X{q*) E[u'^i 



q{X) X{q*) E[u]d 



> 1. 



(4.48) 



The inequality is shown by applying Jensen's inequality to the convex 
function n'^. We conclude that VLD degrades performance for d > 1 in 
the regime of small A, q* relative to FLD with the same mean. 

We next apply the OP LB in Prop. 13.41 to Def. 14.41 to express the 
OP LB for VLD in terms of the MGF of -u'^. 



Proposition 4.9. OP LB as an MGF. For e = and iV = the 

OP LB for VLD is the MGF of the RV -u'^ evaluated at 6' = XcdT^: 



q^^{X) = 1 - E 



1- M[-u'^]{e) 



(4.49) 



We specialize this example to a particular choice of for which we 
can explicitly compute the MGF A^[— u'^]. Fix u to be the contact 
distribution for a PPP, i.e., the random distance from any point x G M!^ 
to the nearest point in H^^^. As motivation, suppose the PPP n^;^^ 
represents locations of base stations, these locations induce a Voronoi 
tesselation of W^, and a Tx in the PPP Ild,x is assigned to the base 
station for the transmitter's cell. The CCDF is immediate from the 
void probability in Prop. 12.61 



Proposition 4.10. Nearest neighbor RV characteristics. The 

nearest neighbor contact RV u for a PPP H^^^ has CCDF: 

Fuiu) = P(u >u)= ¥{Ud,^{hd{o, u)) = 0) = e-^^-^"'. (4.50) 
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The RVs u and u'^ have means and ratio 

which evaluates to 1 (d = 1), ^ 1.273 (d = 2), and 1.404 (d = 3). 
The RV -u'^ has MGF: 

M[-u'm = -^. (4.52) 
The OP LB in Prop. for this choice of U is 

q{X) > g'^A) = 1 - = (4-53) 



We may further specialize this example to the case of 5 = 1/2 and 
employ the exact OP in Cor. 13.11 in Def. 14.41 

For u in Prop. 
1. (4.54) 



Corollary 4.4. Exact OP (e 
liTOl the OP is 



q{X) = 2E 



0, N = 0, 6 



Fig. 14.51 shows six curves of OP vs. A for the case of r = 1, n = 
1/100, d = 2, a = 4: (5 = 1/2). The six curves are the exact, LB, and 
asymptotic OP for both VLD using the nearest neighbor distance u in 
Prop. lilOl and the FLD OP (from Ch. ED with u = E[u] = 5. Observe 
the LBs lie below the exact OP as expected and are asymptotically 
tight as A — )• 0, and that the asymptotic approximations as A — )• 
are valid. Observe the VLD OP is slightly higher than the FLD OP 
for A — )• 0, as expected by ()4.48p . Note however, that the ordering for 
A — )• does not hold for larger A, as we in fact observe FLD OP exceeds 
VLD OP in this regime. Although the example illustrates there may 
be some significant dependence of the OP (and hence the TC) on the 
link distance variability for A (and hence the OP) large, we also observe 
the impact of this variability to be minimal in the small A (and hence 
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small OP) regime, which is typically of more practical interest. In short, 
this section shows that a variable Tx-Rx link distance can be included 
in the OP and TC results without any major modifications; the link 
distance is simply conditioned on and averaged which, for a given link 
distribution, causes a fixed reduction in TC. 




(5*001 0.005 0.01 0.05 0.1 (?.001 0.005 0.01 0.05 0.1 



Fig. 4.5 The OP (exELct, LB, and asymptotic approximation for A — > 0) vs. A for both FLD 
and VLD. The top left plot shows all six curves, various subsets of these are shown in the 
remaining five plots. 



4.3 Multihop TC 



Since an ad hoc network - or any wireless network - may wish to use 
intermediate relays to transfer packets from source to destination, an 
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extension of the TC metric to multiple hops is very desirable. This is 
particularly the case in "power-limited" networks; i.e., those where the 
source-destination distance U (say) is sufficiently large that it cannot be 
bridged with a single transmission, but instead must be broken into M 
shorter hops using other nodes in the network as relays. This is a chal- 
lenging extension, and in this section we overview one recent approach 
that retains fairly good tractability. Building on the TC framework so 
far, in this section we assume: 

(1) All links experience unit mean Rayleigh fading. 

(2) All hops are equidistant of length u = U /M, and on a straight 
line between the source and destination. This can easily be 
shown to be "best case" in terms of success probability. 

(3) A packet is repeatedly transmitted on each hop until it is 
successfully received by the next hop, or until a timeout cor- 
responding to a maximum total number of end-to-end (e2e) 
attempts A occurs. 

(4) Each transmission attempt experiences iid fading and inter- 
ference, and hence has the same OP which is independent of 
all other attempts. 

(5) A packet must reach the final destination before a new one 
is injected by the source; i.e., there is no intra-route spatial 
reuse. The spatial intensity of attemped transmissions at a 
point in time is A (nodes per square meter); this is also the 
spatial intensity of source nodes. 

This model can be conceptualized by Fig. 14.61 Based on these assump- 
tions, we introduce multihop TC, originally defined in [sj^. We acknowl- 
edge prior work [69] that proposed a similar model but focused more 
on conditions for queue stability over the multihop routes. 

Definition 4.5. Multihop TC. Define Amh(A, U, M, A) as the spatial 
intensity of successfully delivered packets (packets per square meter) 
when the intensity of source nodes is A, each source destination pair is 
separated by U meters, each such pair employs M — 1 relays positioned 



^ This paper introduces a slightly different metric called random access transport capacity. 
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Hh = 3 



Fig. 4.6 The multihop TC model with M = 3. 



equidistant on the straight line connecting them, and each packet is 
permitted at most A transmission attempts (end to end) before time- 
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out. Let Tm ~ Pascal(M, g(A) jl be the RV denoting the total number 
of independent transmission attempts required to achieve M successes 
when each trial fails according to the OP ^(A). Then 

where we thin by the probability of successful end to end delivery 
^{Tm < A) and divide by the average number of transmissions re- 
quired for end to end delivery, E[Tm A A]. The multihop TC is defined 
as the maximum of ()4.55p over the number of hops M: 

Amh = A max — -— — . (4.56) 

It has units of packets per unit area. 

Recall \A\ = {1,...,^}. Neither the numerator P(Tj\,f < A) nor the 
denominator E[Ta// A A\ in (14.550 can be computed directly, however, 
a useful inequality is given in the following proposition. 

Proposition 4.11. Multihop TC inequality. For all A G N and all 
MG {!,..., yl}: 

P(TAf < A) 1 

e[Tmaa] - e[Tm]" ^ • ^ 



The proof is given in [5]. The main idea is that as A increases the 
numerator P(TAf < yl) — t- 1 and the denominator E[Tjvf A ^] — t- E[Ta,/] 
both increase, but the numerator does so slightly more quickly, resulting 
in an upper bound. 

With Ass. (2) above (equally spaced hops), the per-hop per-attempt 
OP is simply the standard Rayleigh fading OP (Cor. WA\ with per-hop 
distance u = U/M and snr = P{U /M)-'^ /N: 



special case of the negative binomial distribution. Sometimes these distributions are 
defined as the number of failures (instead of trials) until M successes are achieved. 
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Amh 
0.04 




Fig. 4.7 The multihop TC Xmh from (|4.56| l and its upper bound A^*;^ from jfiSO^i versus 
the allowed number of transmission attempts A. 



Since each transmission attempt is iid, the RV Tm is in fact the sum 
of M independent geometric RVs, Jm = T^^'^'' + • • • T^^\ where T^] is 
the number of independent trials required until success is achieved on 
hop i, and each trial fails with OP q{X,M), i.e., 

P(Tg=t) = g(A,M)*-i(l-g(A,M)), ten. (4.59) 

Therefore, the average number of attempts per hop is E[Tj^] = 
1/(1— g(A, M)) and the expected number of total transmissions required 
to move a packet from source to destination is E[Tm] = -^IE[T^]] = 
M/(l — g(A, M)). Combining all these observations yields the following 
UB on the multihop TC. 



Proposition 4.12. Multihop TC UB. An UB on the multihop TC 
from ()4.56p is 

Amh < AlS, = A max ^ ~ "^[^^ ^\ (4.60) 
for q{\M) in (ITOD . 



The upper bound is somewhat loose for small A and is fairly tight for 
large A, as can be seen in Fig. 14.71 The natural next step is to find 
the optimal hopcount M*{A), i.e., the value of M that maximizes the 
multihop TC. 
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Definition 4.6. Optimal number of hops. The (bound-) optimal 
number of hops using the multihop TC UB in Prop. [^T2] is 

M = are max . (4.61) 

Me[A] M ^ ' 

The following proposition characterizes M* for the special case d = 2. 

Proposition 4.13. Optimal number of hops. The optimal number 
of hops M* in Def. 14.61 for d = 2 is the solution to the equation 

M° - 2Ar2/'"K„C/2M"-2 = Q. (4.62) 

This results in closed-form expressions for M only when a £ {3, 4, 6, 8}, 
in which case M* is the largest positive root of (j4.62p . The constant 
Ka = 7r^(5 csc(-7r5) with 5 = ^. 

Proof. The objective is 

i^^^ = ^ eM-hM-'^ - k,M-^) (4.63) 

where = ^ and k2 = A7r^<5csc(7r(5)r'^C/^, which follows from (14.580 . 
Setting the derivative equal to zero gives 

exp(A;iM-" + k2M~'^) (l - fciaM"" - 2^2^-^) = (4.64) 

which results in 

Af" - 2A;2M°-2 _ j^^f^ = 0. (4.65) 

By the Abel-Ruffini theorem, a formula solution to a polynomial equa- 
tion only exists when the degree of the polynomial is 4 or less, therefore 
M can be found in closed-form only for a € {2, 3, 4, 6}, although it can 
also be found in principle for a € {1,8}. The solutions for a = 6 and 
a = 8 follow similarly to the a = 3 and a = 4 solutions due to the 
sparsity of the polynomial. 

□ 

As noted previously, the expected interference in a 2D PPP is infinite 
for a = 2 (recall d = 2 in Prop. I4.13|) . Hence we give the results 
for a G {3, 4} in the following two corollaries. Proofs follow standard 
algebraic arguments (although perhaps unfamiliar ones for a = 3) and 
are given in [5]. 
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Corollary 4.5. Optimal number of hops for a = 3. Solving (|4.65p 
with a = 3 yields two possible solutions, depending on the polarity of 
the equation's discriminant D 

M*(3) = r^C/ 



4 



27 



. When L» > 0, 



3/3iV 3/3iV 

2P ■' y 2P 



(4.66) 



where 



and Ks = ^VStt^ 
When L> < 



/ 



3iV 



8K^ 



A3 



(4.67) 



2P ) 27 
7.6 is the -K^b csc(7r5) term evaluated at a = 3. 



M*(3) 



2\\ — - — T3 (7cos 



1 



-arc cos 



3\/3A^ 



4\/2P(Aif3)2 



(4.68) 



Although the expressions (j4.66|) and (|4.68|) at first appear to be quite 
different, in fact they have a quite similar dependence in terms of the 
main parameters of interest. 



Corollary 4.6. Optimal number of hops for a = 4. Solving (j4.65p 
with a = 4 yields a unique maximum positive real solution for M*: 

1 / 7r2 \ vr^ iiV 
M*(4) = t^C/WA- + yA2- + — . (4.69) 



The tightness of the bound is shown in Fig. 14.81 where it is seen that 
the quality of the bound improves as A increases. 

We conclude by commenting on the five enumerated assumptions 
that were made at the beginning of this section in order to get to this 
result. The first assumption, of Rayleigh fading, is fairly innocuous and 
used just because it allows a precise equality for per-hop outage prob- 
ability. Any of the other channel models discussed in this monograph 
could of course be used, resulting in appropriately modified bounds. 
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Amh(M) Amh(M) 
0.04 0.04 




Fig. 4.8 The quantity AmhCA, U,M, A) in 1)4. 55p and its upper bound versus the number of 
hops M for each M = 1, . . . , A when A = 6 (left) and A = 12 (right). Note the maximizing 
M* for the exact solution and the UB are close to one another in both cases. 

The second assumption of equi-distant relays on a straight hne is per- 
haps the most physically questionable assumption, especially since the 
network interference is drawn from a PPP. This model is best case, 
so it preserves the upper bound, but it can be legitimately asked how 
loose the bound might be compared to a network with relays that 
must be chosen from a random distribution of nodes. This question 
has been explored in p6], but deserves further study. The third as- 
sumption, of transmission until success or repeated failure (timout) is 
reasonable and does not require modification. The fourth assumption, 
of each slot being iid, could be approximated in a frequency hopping 
or other diversity-harvesting system; in which each slot could truly 
see iid interference and fading, and thus have nearly independent suc- 
cess probability. This can be relaxed, see e.g., [Sj which builds upon 
[5]. And finally, the fifth assumption regarding no intra-route reuse is 
pessimistic for a given source-destination pair, especially for large M. 
However, in a network-wide sense sparse reuse of a given route allows 
other source-destination pairs to transmit and so there is no net change 
to the multihop TC, since this is a spatially averaged metric. Formally, 
the loss in multihop TC from M in the denominator (which assumes no 
intra-route reuse) can be exactly balanced by a corresponding increase 
in A in the numerator, resulting in the same interference intensity. 



5 



Design techniques for wireless networks 



In this chapter we consider four design techniques/issues for a decen- 
trahzed wireless network based on the models of Ch. [3] and HI 

(1) Spectrum management: the network bandwidth W (Hz) is 
to be broken up into B bands each of size W/B (Hz). This 
allows the density of interferers (per band) to be controlled, 
since the density will then be X/B. 

(2) Interference cancellation: allow a Rx to cancel a fraction k 
of the interference generated by the K strongest interferers. 

(3) Threshold scheduling: exploit fading by scheduling transmis- 
sions for Tx-Rx pairs with a strong channel. 

(4) Power control: select the transmission power to partially 
compensate for the channel to the Rx. 

Although this list of enhancements and design issues is by no means 
exhaustive, these four topics are major issues that any systems engi- 
neer designing a centralized network protocol would be faced with. We 
will see that the TC framework is able to provide insight and quantita- 
tive network-level performance analysis regarding these issues, whereas 
other approaches typically have been unable to achieve this, and must 
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either use simulations or greatly simplified analytical models. Chap- 
ter 6 continues this direction, focusing exclusively on multi-antenna 
transmission and reception. 

5.1 Spectrum management 

In this section we study the impact of multiple frequency bands on OP 
and TC [l5]. We consider a bandwidth of W (Hz) that is divided into 
B uniform bands, each with bandwidth W/B (Hz). The objective is to 
select B to maximize the TC. The solution of this optimization is non- 
trivial on account of the following two dependencies: as B increases 
there is simultaneously less interference on each band and a higher 
SINR required to achieve a fixed rate R. The lower interference is on 
account of there being more bands from which to choose, while the 
higher SINR requirement is due to Shannon's formula c = Tyiog2(l + 
sinr); the nature of this dependence will be made clear in what follows. 



Assumption 5.1. Random band selection. Throughout this sec- 
tion fix e = and set each hj = 1 (no fading) in Ass. 12.21 Further, 
assume each Tx will independently and uniformly at random select a 
single band in [B] = {1, . . . ,B} on which to operate. Let bo be the 
random band selected by the reference Tx, and {hi} be the random 
bands selected by each interferer i £ Ild,x- 



Remark 5.1. Thinned interference seen by reference Rx. Hav- 
ing bo and {hi} be uniform and independent implies that the SINRs 
on each band are iid, and hence, without loss of generality, we may 
assume the reference Tx selects bo = 1 (say). It follows that only in- 
terferers that select band 1 are of relevance to the reference Rx at o. 
More formally, the MPPP ^d,x = {(xi,bi)} induces B iid PPPs, each 
of intensity X/B, and the PPP of interferers on band 1 is 



Hd.A/B = {xi : (xi,l) e <^d,x}- 



(5.1) 
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It furthermore follows that the (normalized) interference SN RV seen 
at o on band 1 is 



The random channel capacity and channel spectral efficiency are de- 
fined below, along with equivalent definitions of outage event (|l.ip . 



Definition 5.1. Noise, SINR, SNR, capacity, spectral effi- 
ciency. 

(1) The noise power spectral density is uniform at r] (W/Hz). 
The noise power over a band is N{B) = rj^ (W) and the 
noise power over the full spectrum is = rjW (W). 

(2) The SINR seen at o on band 1 is 

Pu~" 1 
sinr(o) = ^ = ^ , (5.3) 

where the full-spectrum Rx SNR is 

(3) The Shannon channel capacity at o on a channel of W/B 
(Hz), treating interference as noise (Ass. fTT]) . is the RV 

W 

c(o) = — log2(l + sinr(o)), (bps). (5.5) 

(4) The corresponding Shannon spectral efficiency is the RV 

^ = log2(l + sinr(o)), (b/s/Hz). (5.6) 

(5) The rate and spectral efficiency requirements are defined as 

W R 
^=^log2(l+r) (bps), = — S = log2(l + r) (b/s/Hz). 

(5.7) 
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(6) The outage event (II. ip may be expressed in terms of both 
the rate and spectral efficiency: 

sinr(o) <T^ c(o) < R < v. (5.8) 

W/B 



The above definitions and discussion motivate the following modi- 
fications to the definitions of OP and TC. 



Definition 5.2. The OP and TC under multiple bands are de- 
fined as follows: 

(1) The OP at o is a function of the per-band spatial intensity 
of interferers \/B: 

q{X/B, B) = P(sinr(o) < r). (5.9) 

(2) The TC under outage constraint q* G (0, 1) is defined as 

X{q\B*) EE ra^^Bq-\q\B){l - q*), (5.10) 

where the maximization is over all possible number of bands 
B £ N. Here q~^{-,B) is the inverse of the monotone increas- 
ing OP q{-,B) in (I5.9p . The multiplication by B is on account 
of the fact that q~^{q* , B){1 — q*) is the average number of 
successful transmissions per unit area on band 1. 



The above definitions make clear that the TC in Def. 15.11 is the 
following modification of the TC (Prop. 13. 2p in Ch. [3l 



Proposition 5.1. Using Def. 15. li the TC under multiple bands is: 

X{q* , B*) = K{q*)uj{B*) (5.11) 

where 

2(1 -a*) . . 
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and is the inverse CCDF of the RV ^}''/'°(o). Moreover, 

s 



u{B) = B (—^ , (5.13) 

V2#^-l snrS; ^ ' 



with 



B* G arg max wfS), (5.14) 

Be[Li3maxJ] 



for Smax satisfying 



W 

-Bmax = log2(l + Snr^max)- (5.15) 



Proof. Multiplying the TC in Prop. 13.21 bv B as in (jS.lUp gives: 

= ^ L ' (5-16) 

where we write t{B),N{B) to emphasize their dependence upon B. 
Substitution of ()5.4p and ()5.7p and algebra yields the proposition. □ 



Remark 5.2. Optimal number of bands independent of target 
OP. The form (15. lip highlights the fact that the optimal number of 
bands B* is independent of the target OP q* . In fact K{q*) captures 
the spatial components of the network through its dependence upon 
u, d, 5, q* , while uj{B) captures the spectral components of the network 
through its dependence upon R, W, snr, 6. 



It is natural to change the design variable from the number of bands B 
to the spectral efficiency u = ^B. The following definition is central 
to what follows. 

Definition 5.3. Energy per bit and receiver SNR. The energy 
per bit is 



rj rjR 



(J/bit), (5.17) 
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and relates to the Rx SNR as 



R — = Wsnr = 
V 



Pu- 



(5.18) 



Remark 5.3. Relaxation of integrality constraint. Although B 
is naturally restricted to be integral, taking values in [[-BmaxJ] for -Bmax 
in (I5.14|) . by continuity of the objective oj{B) we can safely relax the 
domain to the continuous interval [0, -Bmax] and then take the near- 
est integer. It follows that the corresponding domain for u is [0, fmax] 
defined below. 



Simple algebra gives the following corollary of Prop. 15. li 



Corollary 5.1. Using Def. 15. li the TC under multiple bands is: 



\{q\v*) = —K{q*)u{u*) 



(5.19) 



where Fj- ^ and ^(g*) are as in Prop. 15.11 




(5.20) 



with 



y* £ arg max 



max 



(5.21) 



and t'max determined by ^ and satisfying the equation 




(5.22) 



Note ()5.22p does not have a solution for all ^. 
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Proposition 5.2. Minimum energy per bit required for solu- 
tion. Equation ()5.22p has a solution precisely for 

— > log 2 f« 0.693 f« -1.59 (dB). (5.23) 



Proof. Application of x > log(l + x) for x G M+ to (I5.22p yields: 

(log 2)fmax = log ( 1 + — Z^max ) > — l^max- (5.24) 
\ V J V 

□ 



Remark 5.4. Low SNR regime. The minimum energy per bit 
threshold of —1.59 (dB) is a fundamental quantity in the wideband 
regime. The interested reader is referred to [75] for a general discussion 
and to [35] (footnote 4) for an elaboration in this context. 



Fig. 15.11 shows uj{B) vs. B (I5.13p . uj{iy) vs. u (I5.20p . and two plots 
illustrating z^max(^) (I5.22P using parameters: 

d = 2 n = l(m) aG{3,4} PG{3,6}(dBW) 

R=l (Mbps) W = 10 (MHz) r] = 10"^ (W/Hz) 

(5.25) 

These are the default parameters throughout this section except when 
indicated otherwise. Observe that each of the four (a, P) pairs in Fig. 
15.11 has a distinct optimal B* (resp. u*) that maximizes uj{B) (resp. 
ti;(i/)). The B* maximizing uj{B) in (I5.14|) is seen to depend upon three 
parameters {R/W, snr, 6) while the u* maximizing oj{i') in ()5.2ip is 
seen to only depend upon two parameters (^ and S). For this reason 
we restrict our attention in what follows to optimization over u instead 
of over B, with the understanding that B* = ^i^* (recall Rem. 15. 3p . 
The following result characterizes u* ( |45j Thm. 1). 
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ijj(M) LjUvj 
10 1.0 




Fig. 5.1 Top: the spectral component tu(B) (left) and lo{i/) (right) vs. the number of bands 
B (left) and the spectral efficiency i/ (right) under Prop. [STTI (left) and Cor. lS.lK right) using 
parameters 1)5. 25p . Bottom: the value I'max obeys I'max = log2(l + ^^^max) and thus is 
the intersection of the u line with the log2(l + ^i') curve (left). Equivalently, !^max(^) vs. 
^ (in dB) is shown on the right. No solution t'max exists for ^ < log 2 (si —1.59 dB). 



Proposition 5.3. The optimal spectral efficiency u* in (j5.2ip is 
the unique positive solution of: 

-1/(2'^ - 1) - (1 - 6){2'' - if - 5(log 2)-u^2'' = 0. (5.26) 
r] 7] 

Furthermore, v* is increasing in ^ and decreasing in 5. 



The proof is in |45) and requires only simple calculus. Fig. 15.21 shows 
V* vs. ^ for various 5 as well as the TC A(g*, v*) vs. ^ for various q* . 
Observe that v* is increasing in ^ and decreasing in 5 as asserted in 
Prop. 15.31 Observe that the TC is increasing in both ^ and q* . 

We next study the asymptotic behavior of v* in the high SNR — )• 
oo) and low SNR — ?• log 2) regimes. 
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Fig. 5.2 Left: the optimal u* satisfying l[06)l vs. for <5 e {2/3,2/4,2/5}. Right: the 



TC A((j*,i/*) vs. 



^ for 5 : 



1/2 and q* e {0.01,0.05,0.10}. 



Corollary 5.2. The asymptotic optimal spectral efficiency v* in 

the high SNR regime — )• oo) is the unique positive solution of: 



1 - 2"'^ = (log2)(5i/. (5.27) 



The corollary is immediate upon observing the second term in ()5.26p 
is negligibly small as ^ — )• oo. This function i'*{6) is shown in Fig. 15.31 
Observe the convergence of i^* to that predicted by Cor. 15.21 is slower 
for small 6 (equivalently, for high a). 




Fig. 5.3 Left: the optimal high SNR spectral efficiency w* (S) in Cor. 15.21 Right: the optimal 
spectral efficiency u* vs. — (in dB) from Fig. l5.2l (left) along with the high SNR asymptotes. 



Proposition 5.4. The maximum possible spectral efficiency 

i^max (obeying (j5.22p ) in the low SNR regime (^ — )■ log 2) to first 
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order in — log 2) and second order in v is 

2(^-log2) , /£, \2 



^""'^ (log2)(2£^-log2) ' ' ' \v 

The optimal spectral efficiency (obeying (15.260 ) in the low SNR regime 
(ff ~^ log 2) to first order in — log 2) and second order in u is 

.+e.(.Vc.f--iog2 



(log2)((log2)-(l-25)(f -log2)) 



(5.29) 



Proof. The series expansion of log(l + ^i/) — (viewed as a function 
of {^,1^)) around the point (log 2,0) to first order in (^ — log 2) and 



to second order in is: 



log (l + „ = (zi5i^^2 + 0(^3) 

+o(^|-log2^ . (5.30) 

Cancelling the common factor of from the first two terms, equating 
with zero, and solving for i/ gives ()5.28p . Note the second order expan- 
sion for v is required to be able to solve for z/ — a first order expansion 
in z/ is inadequate in this regard. The series expansion of (|5.26|) to first 
order in ^ and third order in v around {^,1^) = (log 2,0) is 

dos]) = (^-l(iog2)V + o(z/4; 

+ ((iog2)(i - sy + ^(log 2)2(1 - 26y + oy: 

|-log2^ +o(^^-log2y. (5.31) 

Cancelling the common factor of z^^ form the first two terms, equat- 
ing with zero, and solving for v gives (I5.29p . Again, the third order 
expansion for v is required to be able to solve for v. □ 
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Fig. 15.41 (left) shows the exact value of fmax (obeying (I5.22p ) and the 
low SNR approximation from (I5.28P vs. ^ (in dB). Similarly, Fig. 15.41 
(right) shows the exact value of (obeying (|5.26|) ) and the low SNR 
approximation from (j5.29p for 6 S {2/3,2/4,2/5}. Note the validity 
of both the max and optimal approximations quickly degrades as ^ 
increases, and that the optimal approximation is more accurate for 
high 5 than for small 6. 




Fig. 5.4 Exact and approximate values for /''max (left) and u* (right) vs. ^ (in dB) using 
the low SNR approximations in Prop. [53] The plot shows exact (thick) and approximate 
(thin) values for 5 of 2/3 (solid), 2/4 (dashed), 2/5 (dotted). 



5.2 Interference cancellation (IC) 

In this section we study the impact of interference cancellation (IC) on 
the OP and TC [77]. There are (at least) three major constraints in 
designing a Rx capable of cancellation. 



Definition 5.4. The (k, iT, Pmin) IC model is defined as follows. 

(1) Cancellation is inevitably imperfect, due to imperfect chan- 
nel estimation and signal reconstruction. To model this, we 
multiply interference from cancellable nodes by k S [0,1]. 
Note K = is perfect cancellation and /c = 1 is no cancella- 
tion. 

(2) Processing delays and computational complexity constraints 
limit the number of cancellable interferers. Therefore, we can- 
cel at most the nearest K £'N interferers. 

(3) Received power limits the decodability of an interferer, i.e., 
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it is difficult to accurately estimate and subtract low SNR 
interferers, due for example to finite resolution receivers. 
Thus, we only cancel interferers with a received power above 

Pmin £ 



These constraints may be interdependent, e.g., there may be natural 
design tradeoff in (k, K) in that high quality cancellation (small k) may 
require larger delays and thereby UB feasible K. We make the following 
assumptions throughout this section. Recall Ass. I2.1l in §2.21 We retain 
Ass. 12.21 in §2.31 and in particular assume e = 0. We next define the 
SINR seen at o under the {k, K, Pmm) IC model. 



Definition 5.5. The SINR at a (k, K, Pjam) IC capable reference 

Rx located at o is as in Def. 13. H except interference is split into par- 
tially cancelled and uncancelled components: 

" .PZPc(o) + PZ^c(,)+^ ^ (5-32) 
where the partially cancelled / uncancelled interference components are 
ZP^(o)^ |xi|-", ^'^^(0)= Yl l^^^l'" (5.33) 

Here, the set of interferers Ild,x = {xi} is partitioned into the sets of 
partially cancelled and uncancelled interferers based on {K,P^[j^): 

nP3^(o) = {iG Ud,x : i e [K] and P|xr" > ^'min} 

u^%{o) = n,,,\nP3^(o). (5.34) 



Note the i G [K] requirement and recall the labeling convention in 
Ass. 12.11 The following lemma translates the received power constraint 
(-fmin) and cancellation quantity constraint (K) for interferers in Ild,x 
(from Def. l5.4p to a maximum distance from o under Hi i, and the dom- 
inant vs. non-dominant criterion in Ild^x (from Def. 13. 3p to a maximum 
distance from o under IIi^i. 
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Lemma 5.1. Constraint mapping. Satisfying both the received 
power constraint P|xj|^" > Pmin and the cancellation quantity con- 
straint i G [K] for i £ Tid,x is equivalent to |tj| < tP'^ for i G IIi^i 
where 

tP^-|tx|A^f-^)'. (5.35) 

^ \ ^rnin / 

The requirement for a partially cancellable interferer i £ li^'^^{o) to be 
dominant as in Def. 13.31 is equivalent to |ti| < t^^^^^ for i G IIi^i where 

Cn. = ^-'e^ (5.36) 

for ^ in Def. 13.21 The requirement for an uncancellable interferer i G 
ir^\{o) to be dominant is equivalent to |tj| < t'^'^^ for i £ IIi^i where 

t = ^e'- (5.37) 



Proof. All three quantities follow from Prop. [27F1 First, (j5.35p : 



P|xi|-°>P^in4^|xi|^< ^\t,\<^[-—\ . (5.38) 

\ -'min / \ ^min / 

Take the minimum with {txl to also ensure i G [K]. Next, ()5.36p : 

^""'^ >r«^|xi|^<Ac^e'^|ti|<^^¥. (5.39) 



The constraint in (j5.37p is obtained by an identical development but 
without the k on the interference term. □ 



Note tP'^ is a function of the RV |t/^| and is therefore a RV. The 
next lemma applies Lem. 15.11 to define the sets of dominant and non- 
dominant nodes (again in the sense of Def. 13. 3p among the sets of 
cancellable and uncancellable nodes using thresholds i"^^ , t^o^, t^om- 
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Lemma 5.2. The set of partially cancellable / uncancellable 
nodes in Hi i is: 

nP'i(o) = {ie ni,i : |t,| < tP^}, n^^i(o) = ni,i \ ul';,{o), (5.40) 

for tmax in Lem. 15.11 The set of dominant and non-dominant partiahy 
canceUable and uncancehable nodes is: 

<,(o) = {^E^l,l:|t.|<min{tP^^P^J} 
nP'l(o) = {^eU,,,:tZ^<\u\<t^^} 

n?^i(o) = {i e ni,i : tp= < |t,| < Cm} 

ny^i(o) = {iG^l,l:max{tP^Cm}<|t.|} (5-41) 

for imax°™ , imax°™ in Lem. 15.11 The aggregate dominant and non- 
dominant partially cancellable and uncancellable interference is: 
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(5.42) 



Proof. The sets in (j5.40p and ()5.4ip are immediate from Lem. 15.11 The 
aggregate interferences are simply the summed interference contribu- 
tions under the corresponding interferer sets, using the SN RV mapping 
in Prop. EEl □ 

These sets are illustrated in Fig. 15.51 We see that the RV t^^ is a function 
of |tj<-|, the distance of the Kth nearest Rx from o in IIi^i. The next two 
results characterize this RV. The following theorem (from |35j . Thm. 
1) extends Prop. 12.61 
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UC UC PC PC PC uc uc 
non dom non dom non dom non 



lUC ,pc ,pc /UC 

^''doin "''doni ''dom ''dom 
\ \ \ \ \ \ \ t 

-4 -< 



Uncancellable Partially cancellable Uncancellable 

Fig. 5.5 The thresholds t^'^ partitions R into the region of partially cancellable (PC) nodes 
and uncancellable (UC) nodes. The thresholds t?'^ , t further divide each of these regions 

^ dom ' dom ^ 

into dominant and non-dominant regions. 



Theorem 5.1. Ordered distances marginal distributions ([35]). 

The successive (ordered) random distances from o of points in Ild,\ = 
{xfc} with |xi| < |x2| < • • • have marginal (generahzed Gamma) PDFs: 

^|><^|^*) = friiir^"'''*'' * e 1^+' ^ ^ (5.43) 



By the mapping theorem we will have need only for d = 1 and A = 1. 

Corollary 5.3. Ordered distances in Hi i marginal distribu- 
tions. The successive (ordered) random distances from o of points in 
IIi^i = {tfc} with |ti| < |t2| < • • • have marginal (Gamma) PDFs: 

= ^(Hl)!^''*' i e M+,fc G N. (5.44) 



Moreover, 



r(fc+p) 



, p > —k 



E[\tk\P] = { 2nfc-i)!' (5.45) 
oo, else 



In particular IE[|tfc|] = |. 



The PDF is immediate from Thm. 15.11 The moments are obtained by 
calculus. The PDFs for k G {1, . . . , 5} are shown in Fig. 15.61 The main 
result of this section is the following LB on the OP. 
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ftk(t) 




Fig. 5.6 The PDFs f\t^\it) for fc e {1, . . . ,5} from Cor.ES] Note E[|tfc|] = |. 



Proposition 5.5. The OP LB under (k, A', Pmin) IC Rx model is: 

g^A) = 1 - E [exp { -2 min {tP^ } exp {-2( t - t^T }] 

(5.46) 

where = max{x,0}, and tP'^, ijom' *dom given in Lem. [STTl tP*^ 
is a RV that is a function of the RV {txl (from Prop. 15. 3p . and the 
expectation is with respect to |tx|- 

Proof. We first apply the SINR definition from Def. 15.51 and solve for 
the interference: 

q{X) = P(sinr(o) < r) 

< T 



PZ(o) + N 
^ o > — - ^ 



czp;=f'°(o) + z:;X'°(o)>r") (5.47) 

for ^ in Def. 13.21 We apply the law of total probability by conditioning 
on Iti^l which effectively partitions M into the regions in Fig. 15.51 With 
those regions fixed for each possible value of {txl we achieve a LB by 
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only considering the dominant interferers (both partiahy cancellable 
and uncancellable) . 



> 



E 
E 



ypc,a,0 



o) 



^^Ix (o) + Z,^ • (o) > C 



|tA' 



(5.48) 



We then take the complement of the outage event, yielding the event 
that the sum dominant partially cancellable and dominant uncan- 
cellable interference is less than By construction this event holds 
iff there are no dominant interferers of any type. 



,1b 



(A) 



E 
E 



ypC.O.O, 



o) + z-r'°(o)<r" 



0^^1(0) = 0nnri(o) 



(5.49) 



Because the dominant partially cancellable interferers and dominant 
uncancellable interferers are defined on disjoint regions of R, it fol- 
lows by independence property of the PPP that their occupancies are 
independent RVs. 



(5.50) 



Prop. 



with d = X = 1 applied to the two dominant regions gives 



□ 



Fig.E71and[53]plot the OP LB q^^{X) from Prop.ESlfor various choices 
of the (k, K, Pmin) SIC Rx model. The default parameters are 



rE {1,5} u = l a = 4 d = 
A = 0.025 k; = 0.05 K = 3 Pmin 



(5.51) 



2 P = 1 iV = O 
= 1 = O(dBW) 

Several points merit comment. First, note that in all cases Fig. 15.71 
shows the corresponding performance without SIC using the OP LB 
q^^{X) from Prop. 13.41 (with e = 0). Recall that this bound was shown 
to be tight. Second, in most (but not in all) cases the plots show greater 
sensitivity of the OP LB to the IC parameters for r = 1 (right) than for 
r = 5 (left). Third, the plots illustrate the benefit in reducing the OP 
from improved IC is quite parameter dependent; we now give several 
examples of this. In the top plots of Fig. 15.71 we see that foi K = 3 
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Fig. 5.7 The OP LB q^^{\) from Prop. [53] for the (K.^.P^in) IC Rx model for r = 5 
(left) and r = 1 (right). The OP is shown vs. A for varying ft (top), K (middle), and Pmin 
(bottom). Default parameters are in 115. 5H . In all cases we also show the OP LB g"'(A) 
from Prop. [T4l for no IC for the same parameters. 



and Pmin = 1 (0 dBW) the performance of IC with k = 0.25 or more is 
indistinguishable from no IC for r = 5, but this is not true for r = 5. 
In the middle plots of Fig. 15.71 we see the OP is insensitive to K for 
K > 1 when r = 5, and this also holds for r = 1 and A small, but not 
for A larger. In the bottom plots of Fig. 15.71 we see the OP for Pmin = 6 
or greater is the same as no SIC for r = 5, but that the OP under 
.fmin = —6 is the same as for Pmin = for r = 1 (all powers in dBW). 
In the top plots of Fig. 15.81 we observe that improving the cancellation 
quality k by an order of magnitude reduces the OP by only a few 
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Fig. 5.8 The OP LB g"'{A) from Prop. [53] for the (K,/f,Pmin) IC Rx model for r = 5 
(left) and t = 1 (right). The OP is shown vs. k for varying K (top), vs. k for varying Pmin 
(middle), and vs. Pmin for varying K (bottom). Default parameters are in 1)5. 5ip . In all 
cases we also show the OP LB q^^{X) from Prop. [X4l for no IC for the same parameters. 



percentage points in the given parameter regime. In the bottom plots 
of Fig. 15.81 we observe again limited dependence of OP upon K > 1, 
while the right plot (r = 1) shows insensitivity to Pmin below dBW. 



Remark 5.5. Other bounds on OP and TC under IC are found 
in the literature. We specifically discuss [58] centered upon the OP LB 
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for perfect (k = 0) IC: 



> y 



> 



\>^K+l\ 



> y) 



xi^+i| < 



(5.52) 

where Ass. 12. II dictates the labeUng of the points. The CDF from Thm. 
l5.1l for k = K + 1 may be used to give an exphcit OP LB. Note that this 
bound holds also for imperfect IC (k G (0, 1)), but is not as insightful as 
the bound depends upon K but not on k. Perfect cancellation naturally 
leads also to an UB on OP using the Markov inequality, as done in jl6] 
(Thm. 1 and Lem. 2), where the essential step is below: 



>y\ < -E 

y 



.i=K+l 



^ oo 

i y E 



X,; 



(5.53) 

The moment E[xj|~"] may be computed using the PDF in Thm. 15. H 
and this in turn be UBed using Kershaw's inequality on the Gamma 
function. 



5.3 Fading threshold scheduling (FTS) 

In §4.11 we extended the basic model of Ch. [3] to incorporate fading, 
and Cor. 14.31 showed that the asymptotic OP and TC (as A — )• and 
g* — )• 0, respectively) under fading were worse than without fading 
[75] . The assumed slotted Aloha MAC protocol does not exploit fading 
as a source of diversity, and this (partially) explains why OP and TC 
degrade under fading. In this section we seek to exploit fading through 
the use of scheduling; in ^5.41 we will exploit fading by selecting the 
transmission power. Only certain forms of scheduling will retain the 
critical features of the basic model that allow for analytical tractability, 
namely, the PPP independence property that the number of points in 
disjoint regions of are independent RVs. Scheduling in its usual 
usage refers to inducing a negative spatial and or temporal correlation 
for activity among adjacent nodes so as to minimize the occurence of 
collisions, and in this sense the set of active nodes under scheduling will 
not form a PPP. Instead, our usage is restricted to "fading threshold 
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scheduling" (FTS): a MAC protocol where potential transmitting nodes 
decide to transmit when the fading coefficient of the channel to their 
intended Rx exceeds a threshold. As we assume independent channels, 
it follows that the transmission decisions of two nodes are independent. 
The following definition formalizes the model for this section. 



Definition 5.6. Fading coefficients, signal interference, and 

SINR. Extend Def. 14.11 throughout this section as follows. Let 
{(hi,i, hj^o)} for i € {0, 1, 2, . . . , } be iid RVs indicating the fading coef- 
ficient for the channels between each Tx i and its Rx and between each 
Tx i and the reference Rx at o: 

(1) ho,o: fading coefficient from reference Tx to the reference Rx. 

(2) hj_0) « G N: fading coefficient from interferer i to reference Rx. 

(3) hi^f. fading coefficient from interferer i to its own Rx. 

Let Ild Xpai denote the set of poientia/ transmitters (c./. Rem. [TTT]) . and 
assume each Tx elects to transmit precisely when its fading coefficient 
to its intended Rx exceeds a threshold h £ M-|_. The FTS threshold 
h must lie in the support of the RV h, and moreover we assume the 
threshold to be such that the success is guaranteed in the absence of 
interference, i.e., {c.f. Rem. | 



h N ^ T , , 

— >0^h> — . (5.54) 



rti" P snr 

More formally, the MPPP $d,Apot = ^i^o, Ki)} induces a MPPP of 
actual interferers 

of homogeneous intensity A = ApotIP(h > h). Assume the reference Tx 
attempts transmission (since otherwise there is no chance of outage at 
the reference Rx), so that the channel coefficient between the reference 
Tx and Rx ho o has a CCDF 



Fr{h) = P(h >h\h>h) = { ^hW - (5.56) 
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The signal, interference, and SINR seen at the reference Rx at o are 
given by 

S ( 

sinr(o) = ^^^^ + jy/P ' ^^"^ ^ KoW", ^(o) = Yl Ko\^i\~'^^ 

(5.57) 

where {\r\ifl} are iid with distribution F^, ho,o is independent of {hj^o} 
and has distribution ()5.56p , and ^ is as in (I5.55P . 



In this section we modify the OP (Def. [TI]) and TP (Def. 
definitions as follows. 



Definition 5.7. OP and TP. The OP is the probability of success at 
o conditioned on having a channel fade above the threshold. Because 
the threshold parameter h directly controls the intensity of attempted 
transmissions (A = Apot-Fh(^)) as well as the distribution on the signal 
strength ho,0) it is more natural to define the OP as a function of h 
instead of Apot or A. 

qCh) = P(sinr(o) < T|ho,o > h), (5.58) 

where interference is from a PPP ^ of intensity A. Likewise, the 
MAC TP is the spatial intensity of successful transmissions 

A{h) = X{l-q{h)). (5.59) 



Remark 5.6. Quantity vs. quality of transmissions through 

FTS. The design objective in this section is to maximize the TP A(/i) 
over h. Increasing h has two effects: the signal power increases on av- 
erage, and the intensity of attempted transmissions decreases. In other 
words, increasing h increases the quality but decreases the quantity of 
attempted transmissions. Viewing A as the quantity and 1 — q{h) as 
the quality of transmissions, we see the TP is the product of quantity 
times quality, and thereby yields a natural optimization over h. There 
are also fairness and delay costs under h: nodes that happen to have a 
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fading coefficient above the threshold are effectively given priority over 
those with a fading coefficient below. The fairness and delay issues is of 
particular concern for channels with long coherence time where waiting 
for a "good fade" may be prohibitively costly. 

Recall ^4. II addressed OP and TC under fading from three perspectives: 
exact ( ^iXT]) . asymptotic ( WL2\i . and bound results. This 

section follows the same outline, but focuses on OP and TP. 

We first explain why FTS spoils the analytical tractability of §4.1.11 
Recall the derivation of the explicit expression for the OP and TC when 
the signal fading is Rayleigh (ho ~ Exp(l)) in Prop. 14.21 was obtained 
by conditioning on the interference, applying the exponential CCDF, 
and recognizing the resulting expression as the LT of the interference, 
given in Prop. 14.11 If hg ~ Exp(l) then 



for y = — N/P. The key difficulty is the indicator function within 
the expectation that precludes using the LT of 51 in Prop. 14.11 

We next turn to establishing the asymptotic OP (as A — 0) under 
FTS. The following result is the FTS analogue of Prop. WM 




h> h 
else 



(5.60) 



The analogous development to Prop. Br2] gives (writing Y. = Z'^'t'): 



q{X) = P(sinr(o) < r) 

= P(ho,o >Tn°(Z + iV/P)) 

= E[P(ho,o >™"(Z + iV/P)|z)] 

= E[exp{-(Tn"(Z + iV/P)-/i)}l5:>jj 



(5.61) 



Proposition 5.6. Asymptotic OP under FTS (/i — > oo, A — > 0): 



( 



) 



5' 



q{h) = CdE[h^]E 



ho,o 



N 



(5.62) 



TU' 



a 



P 
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In the no noise case {N = 0) this expression becomes 

q{h) = CdT^n'^E[h^]E[hoJ]ApotFh(/i) + 0{\^). (5.63) 
In the no noise and Rayleigh fading case (hg and {hj} exponential RVs): 
qCh) = Qr^n'^r(l + 5)r(l - 6, h, oo)Apot + ©(A^), (5.64) 
for T{z,ti,th) the incomplete Gamma function defined in Def. | 



Proof. The proof is analogous to that of Prop. IT6l with the notable ex- 
ception that the assumption on h in (|5.54|) ensures outage is impossible 
without interference; this leads to a simplified development below. 



qCh) 



E 
E 



sinr(o) < r|ho,o > h 

sinr(o) < T|ho,o 

ho,o N 
P 



d,x ru" 



'0,0 



(5.65) 



Now apply Prop. IXH and 113] noting that ho,o and ^ are independent. 



E 

E 
E 



^Vt" > 



> y "0,0 
E[h^]r^ + 0(y"^^) 



ho,o 


N 










" \ 2 





'0,0 



ho,o N 

TU'^ P 



(5.66) 



Simplification yields (|5.62|) . (|5.63|) is immediate upon substituting 
0. Assuming Rayleigh fading we obtain (j5.64p by observing 



E[hoo]^h(/i) = r(l-5,/i,oo). 



(5.67) 
□ 



This result is applied below to obtain the asymptotic TP under FTS 
and its optimization over thresholds h. For simplicity we restrict our 
attention to the no noise case in what follows. 
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Proposition 5.7. The asymptotic TP under FTS and no noise 
(iV = 0) as /i ^ oo (A ^ 0) is: 



Aih) = XpotFuCh) (l - 6ApotE[ho5Fh(/i)) + ^(A^) (5.68) 
for b = CdT^u'^'E[h^]. The asymptotic TP optimal h satisfies 

n^-\>h]+nh-\^',] = j;^. (5.69) 



Proof. Expression ()5.68p follows upon substitution of ()5.63p into (|5.59p . 
The derivative of the TP with respect to h is 



A'{h) = -/h(/i) l-6ApotE[ho^]Fh(/i) 



Apot 

-6ApotFh(A)4r |E[ho^]Fh(/i)| • (5.70) 
ah ^ ' ^ 

Observe that ()5.56p gives: 

nKfi]F^Ch) = I h''Mh)dh = E[h-\^fJ 

J h 

i{E[hoJ]Fh(A)} = -h~'Mh) (5.71) 

Substitution, equating with zero, and rearranging gives ()5.69p . □ 

Fig. 15.91 presents some numerical results for Prop. 15.71 the left plot 
shows the asymptotic TP (I5.68P vs. the fading threshold h, and the 
right plot shows the optimal threshold /lopt solving (|5.69p . The optimal 
h trades off quantity with quality of attempted transmissions — as 
expected both the optimal threshold and the associated optimum TP 
increase with Apot as the MAC can afford to be more selective on link 
quality without sacrificing the quantity of transmission attempts. We 
next show that the asymptotic expressions for TP without fading, with 
fading but without scheduling, and with fading under FTS are ordered 
when the signal and interference fading distribution are equal. This 
result extends Cor. 
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b Apot 



Fig. 5.9 Left: the asymptotic TP A(h) Il5.68| l vs. the fading threshold h for Aopt S 
{1/10, 1/5, 1/2}. Right: the TP asymptotic optimal threshold hopt that maximizes 1 15.681 1 
as a function of bAopt- All fading is Rayleigh; other parameter values are d = 2, a = 4 



(S = 1/2), 



5, ti = 1 (6 3.11). 



Proposition 5.8. FTS exploits fading to improve performance. 

Assume the signal and interference fading coefficients are equal in dis- 
tribution: ho,o = h. Assume no noise {N = 0) and define a = ^CdT^u'^. 
Then the asymptotic TPs (as A — )• and h — )• oo) without fading (NF, 
Prop. [3^ . with fading but without scheduling (FNS, Prop-HSj), and 
with FTS (Prop. EZD are: 

AN^(A) = X{l-aX + 0{X^)) 
A™S(X) = A(l-aE[h^]E[ho^^]A + 0(A2)) 

A^^^{X) = A(l-aE[h^]E[ho^^]A + 0(A2)), (5.72) 

where A in A^"""^ implies selecting h so that Apot-Fh(/i) = X, i.e., h = 
P-^ (A/Apot). Large h {h > E[h^]^ is sufficient) ensures 

E[h^]E[hoJ] < 1 < E[h^]E[hoJ], (5.73) 
which in turn ensures 

A™S(A) < aNF(a) < A^^^iX). (5.74) 



Proof. The asymptotic TP expressions for A^^(A) and A^^^(A) are 
immediate from the TP definition A(A) = A(l — q{X)), and the ex- 
pression for A^^^(A) is immediate from Prop. ISTTl using the assumed 
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XpotFh{h) = A. The second inequalities in (|5.73p and (j5.74p are not 
asymptotic, they hold for all A. To prove the first inequalities, note the 
assumption > E[h'^] ensures the second inequality below: 

E[h'5]E[ho ^] = E[h^] < E[h^] _ < 1. (5.75) 

-^ho,o(^) -^ho,o(/i) 

□ 



Fig. 15.101 illustrates the asymptotic TP under no fading, fading with- 
out scheduling, and fading with threshold scheduling. The figure makes 
clear that threshold scheduling transforms fading from an overall per- 
formance penalty (fading without scheduling being inferior to no fad- 
ing) to a performance enhancement (fading with threshold scheduling 
being superior to no fading). Threshold scheduling successfully exploits 
the available fading diversity. Note as a — )• 2 (5 — )• 1) that fading with 
threshold scheduling is equivalent (in terms of asymptotic TP) to no 
fading; one can also show that as a — )• oo ((5 — )• 0) fading without 
scheduling is equivalent to no fading. 

A(X) A(A) 



0.06 




Fig. 5.10 The asymptotic (A — > 0) spatial TP A(A) vs. the spatial intensity of attempted 
transmissions A under i) no fading, ii) fading without scheduhng, and Hi) fading with 
threshold scheduhng for o = 4 (left) and a = 2.5 (right). Other parameters are d = 2, 
« = 1, T = 5, Apot = 1. 



Finally, we establish a LB on the OP and an UB on TP under fading 
threshold scheduling. These bounds follow immediately by the natural 
extension of Def. 14.21 and its subsequent Prop. 14.71 in §4.1.31 replacing 
ho with ho,0) yielding the following result. 
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Proposition 5.9. OP LB and TP UB. Extend the concepts defined 
in Def. 14.21 to apply to threshold scheduling by replacing ho with ho,o- 
Then Prop. iZl yields the OP LB 

-6^ 



,1b 



1 - E 



exp<^ -ApotFh(/i)cdE[h^ 



np.o 



N 



(5.76) 



where the outer expectation is w.r.t. the random signal fade ho,o- In 
the case of no noise {N = 0) the OP LB is: 



Jb 



(/i) = 1 - E 



exp{-Apoti^h(/i)cdr^n'^IE[h^]ho,o}] , (5.77) 



In particular, with no noise the OP LB is expressible in terms of the 



MGF of the RV -h^ g at a certain 9: 



In all cases the TP UB is 

, ub 



A''''{h) = X^Mh)(l-q'\h) 



(5.78) 
(5.79) 



Fig. 15.111 compares the TP UB from Prop. [5^ with the asymptotic TP 
from Prop. [5771 In closing, we refer the interested reader to recent work 




Fig. 5.11 The TP UB and the asymptotic (A — 5> 0) spatial TP A(A) vs. the spatial intensity 
of attempted transmissions A for a = 4 (left) and a = 2.5 (right). Other parameters are 
d = 2, u = I, T = 5, Apot = 1. 



on threshold scheduling by Kim, Baccelli, and de Veciana [H], and Cho 
and Andrews [18]. 
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5.4 Fractional power control (FPC) 

In §5.31 we found that exploiting fading through threshold scheduling 
improved performance (measured as spatial TP) compared with no fad- 
ing; in this section we exploit fading through a specific form of power 
control (PC). In general, PC in decentralized wireless networks refers 
to an algorithm where each node i at each slot t updates its power Pi^t 
based on feedback so that the nodes jointly converge in a distributed 
manner to a globally (or at least locally) optimal power allocation vec- 
tor. As with scheduling in ^5.3|, analytical tractability requires we adopt 
a rather restricted definition of PC. Again, the key property we seek 
to maintain is that the set of interferers at some point in time form a 
PPP, and this requires each interferer act independently of the other 
interferers. Consequently, power control in this section refers to each 
interferer i S 11^ ,^ selecting a (random) power Pj that is a function of 
the fading coefficient to its intended Rx hj^j [Uj. The following defini- 
tion of fractional power control (FPC) makes these concepts precise. 



Definition 5.8. Transmission powers and SINR. Consider the 
fading model and MPPP used in Def. 15.61 in ^5.31 



$d,A = {(xi,hi,i,hi,o)} , (5.80) 

and the same reference pair fading coefficient ho,o- As in Def. 15.61 all 
fading coefficients are independent. Assume {hj^j} are iid with CDF 
-^hi.i' {hj,o} are iid with CDF Fh^^, and ho,o has CDF Fho^. Let / G M 
be the FPC exponent and let P be the average transmission power. 
The reference Tx (0) and each interferer i G Ild,\ select a random 
transmission power according to 



So^ = (5-81) 

The SINR at o is 
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where the received signal and interference powers at (o) are 

'^'^ " ^ h-/h,o|xr". (5.83) 



To clarify, S(o) is proportional to hg q'' due to FPC (Hqq) and fading 
itself (ho,o)- Interference from i seen at (o) is proportional to h~/ due 
to FPC and to hj^o and due to fading. Observe the transmission powers 
are selected so that E[Po] = IE[Pi] = P for all i G I^d,x- Two special 
cases for / merit mention: (i) / = is no PC, i.e., Pq = Pj = P, and 
we revert to the results up until now; (ii) / = 1 is channel inversion, 
i.e., the transmitted signal power inverts the fading coefficient so that 
the received signal power is not dependent upon the fading coefficient. 
Therefore, FPC can be viewed as a fairly general and flexible form of 
power control. The following lemma gives the moments and variance 
of P under Rayleigh fading. 

Lemma 5.3. Transmitted power RV moments under FPC. For 

Rayleigh fading (ho,o> ~ Exp(l)) the moments and variance of Pq, Pi 
under FPC in (I^IHTTl are 



p 



"rii-pf), pf<i 



E[PP] = \ Vr(i-/)y ^•'-^ , (5.84) 

[ oo, else 

Var(P) = ( (S^TF - ^) ' ^ < 

oo, else 



The variance Var(P) is plotted in Fig. 15.12] vs. the FPC exponent /. 

The design objective in this section is the selection of the FPC 
exponent / E M. The outline of this section mirrors that of §5.31 we 
first discuss asymptotic OP (as A — )• 0) and TC (as q* — )• 0), and then 
discuss a LB on OP (UB on TC). Although in ^5.3l we assumed h to be 
sufficiently large to preclude the possibility of a bad fade guaranteeing 
outage even in the absence of interference (c./. Rem. 14. 4p . under FPC 
there is again this possibility. We update Rem. 14.41 for FPC below. 
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Fig. 5.12 The variance of the transmitted power Var(P) vs. the FPC exponent / for P = 1. 
Note Var(P) = oo for / > 1/2. 



Remark 5.7. Fading and outage with no interference. With ran- 
dom signal fading under FPC there is the possibility of a bad fade 
causing outage even in the absence of interference. This outage event 
£qj has the three equivalent forms: 

— — -<T^hoo<{ — lEhnn ^ ? r: < 0, 

N Vsnr ^ O'OV E[h-/]™" P 

(5.86) 

for snr defined in Ass. 13.11 and has probability 

qf{0) ^F{£oj) = Fu,,, (^(-I-E[ho-^])^) . (5.87) 

Note g/(0) is the OP evaluated at A = 0. Denote the complement of 
£oj by £oj, and its probability by 

97(0) = 1 - 9/(0) = P(^o,/) = i^ho,o ((;^IE[hoj) ^) . (5.88) 

All analysis of OP (and hence TC) must therefore condition on ho,o 
being above or below ^i^IE[hQ q]^ ^ ^ to distinguish between the case 
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of outage being possible vs. outage being guaranteed. The range of the 
OP q{X) is [g/(0), 1], and the domain of q* in the TC X{q*) is [g/(0), 1]. 

The following result gives the asymptotic OP and TC under FPC. 

Proposition 5.10. The asymptotic OP under FPC as A — is: 

(?(A) = l-(l-aA)g/(0) + O(A2) 

for 



a = CrfE 



E 



'1,0 



E 



E 



N 



The asymptotic TC under fading as q* — )• qf{0) is 



(5.89) 



(5.90) 



In the no noise case (A^ = 0, (7/(0) = 0) these expressions become 



g(A) = XcdT^u'^E 
X{q*) = 



E 



'1,0 



E 



"i,i 



E 



'o,o 



E 



,-{i-/)5 
'o,o 



''1,1 



E 



''1,0 



E 



'^1,1 



E 



'o,o 



E 



,-(i-/)5 

'0,0 



+0{q*f (5.91) 
For no noise and Rayleigh fading (ho,0) ^ifl^ exponential RVs): 

qix) = AQrVr(i + 5)r(i - /5)r(i - (i - m + ©(a^) 



A(g*) 



QrMr(i + d)r(i - /<5)r(i - (1 - /)(^) 

provided 6, f6, (1 — f)6 are all in (0, 1). 



+ 0(^*)f5.92) 



Proof. The proof is analogous to that of Prop. 14.61 We condition on 
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ho,o and isolate the random normalized interference. 

q{X) = P(sinr(o)<r) 

= P (sinr(o) < t\£oj) qj{0) + 1(1 - qj{0)) 
= 1- (l-P(sinr(o) <T|4/))g/(0) 

= 1 - 1 - E[F(sinr(o) <r|ho,o)|goj]^ qf{0) (5.93) 

V m J 

<r 

P , -/, I ^ 

-f^ l^iella.x i,i "ifil^i 



E[-] = E 



E[h-{] 



E 



'0,0 



-a + N 



'0,0 



0,/ 



P 



10,0 



for 



Define the RV 



''o.o 



N 



as a function of the RV ho,o and apply Prop. to obtain 



(5.94) 
(5.95) 

(5.96) 
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Now apply Prop. [13] to obtain the asymptotic CCDF 
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EM = E 
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hiihi,o 
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(5.97) 



(5.98) 



Note the equivalence of the events £qj and y > 0. Substitution of (|5.96p 
into ([OH]) yields ([SlM]) . Solving for A gives ([OO]) . Expressions ([OT]) 
are immediate after substitution of = 0. Expressions (I5.92p follow 
after application of Lem. 14.11 □ 
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The next result gives the asymptotic optimal choice for / in the special 
case of no noise and when the signal fading coefficient ho,o is equal in 
distribution to the interference fading coefficients hi^i. 

Proposition 5.11. Asymptotic optimality of / = 1/2. Setting 
/ = 1/2 minimizes the asymptotic OP and maximizes the asymptotic 
TC under no noise {N = 0) and ho,o = hi,i- 



Proof. It is clear that the optimal / for minimizing the asymptotic 
OP equals the optimal / for maximizing the asymptotic TC, hence we 
restrict our attention to minimizing the asymptotic OP. The asymptotic 
OP under the assumption of no noise and ho,o = (which we denote 
in this proof by h) gives the optimization problem 



minE 

/6M 



(5.99) 



We will show the following more general result: for any RV x the func- 
tion g{f) = E[x~'^]E[x~-^] for / G R (and / = 1 - /) is convex in / 
and has a unique global minimizer at / = 1/2. The proposition then 
follows by using this result for the RV x = h'^. Recall that a function is 
said to be log-convex if it satisfies 

g{cfi + c/2) < gihfgihf. Vc g [0, i], /i, /2 e R, (5.100) 

where c = 1 — c. Establishing log-convexity of g is sufficient to estab- 
lish convexity of as log-convex functions are convex. Recall Holder's 
inequality asserts 

1 



E[xiX2] < E [x?]t E [x?]5 , V(p, g) : i + ^ = 1. 



(5.101) 



We apply Holder's inequality twice with pi = p2 
9{ch + c/2) = E[x-(^^i+^^2)]E[x-(i-(^^i+^-^2))] 

= E[x-^^l^ 



i and qi = q2 
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As g is convex, its global minimum (if it exists) is found by equating 
the derivative of g{f) with zero and solving for /. The derivative is: 



.'(/) = ^\ 
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(5.103) 
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It is clear that g'{f) = for / = 1/2, establishing the result. 



Remark 5.8. Optimal FPC may incur large power variance. 

Prop. Ism establishes / = 1/2 is the FPC exponent that optimizes the 

asymptotic OP and TC (under the assumptions = and ho,o = hi,i), 
while Lem. 15.31 establishes Var(P) is only finite under Rayleigh fading 
provided / < 1/2. Thus at / = 1/2 we optimize the asymptotic OP 
and TC but the variability of power required is large, and in practice 
the radio power constraint would limit the ability of nodes with poor 
fades to achieve (j5.8ip . Note that the asymptotic OP and TC under 
the above assumptions are symmetric in / around 1/2, i.e., the OP at 
/ equals the OP at 1 — / (c./. ()5.99p ). and in particular the asymptotic 
OP and TC under / = equal that under / = 1. Although no PC 
(/ = 0) has the same asymptotic performance to channel inversion 
(/ = 1), the former is greatly preferred to the latter on the basis of 
power variance (Lem. [53]) . 

Interestingly, recent work modeling and analyzing Qualcomm's 
FlashLINQ peer-to-peer protocol [7j, which has a CSMA-type MAC, 
found a similar result. Namely, with a different interference and trans- 
mission model (namely, SIR-based scheduling) , but also using stochas- 
tic geometric tools, they found that inverse square root power control 
is optimal for the FlashLINQ protocol, which they verified with ex- 
tensive simulations. This shows that this earlier result on FPC in a 
simpler Aloha-type MAC was fairly robust to the proposed model, and 
provides some degree of confidence that the other design-type results 
given in this chapter will be as well. 
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Our next result gives an OP LB and TC UB under FPC. We modify 
the definition of dominant interferer and interference (Def. 14. 2p for 
FPC. 

Definition 5.9. Dominant and maximum interferers. An inter- 
ferer i in the MPPP ^d,\ in Def. 15.81 is dominant at o under threshold 
r, signal fade ho,0; and FPC exponent / if its interference contribution 
is sufficiently strong to cause an outage for the reference Rx at o: 



(5.104) 

Else i is non-dominant. The set of dominant and non-dominant inter- 
ferers at o under (r, ho) is 

l>d,A = i G <5d,A : h,:;/h,,o|xr" > lE[h J] ( J"'/ - ^) I ' 

(5.105) 

and ^d,\ = ^d,\ \ ^d,x- The dominant and non-dominant interference 
at o under (r, hg) 

= E h-/h.o|xr", ^t'(o) = Yl K-/h.o|xr" (5.106) 

are the interference generated by the dominant and non-dominant 
nodes. Note ^^^(o) = ^dl^o) + 



Proposition 5.12. The OP LB under FPC is 
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where £qj and ^/(O) are defined in Rem. [^3] and 
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(5.108) 



120 Design techniques for wireless networks 



In the case of no noise {N = 0, (7/(0) = 0) the LB is: 
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In particular, with no noise the LB is expressible in terms of the MGF 
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0,0 



(0) 



e=bx 



(5.111) 



Proof. The proof requires only minor adaptation of that of Prop. 14.71 

q{X) = P(sinr(o) < r) 

= P (sinr(o) < t\£oj) g/(0) + 1(1 - qf{0)) 

1 - P (sinr(o) < t|^o,/) | 9/(0) (5.112) 
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(5.113) 



and T.'^'^{o) in (I5.95p . Lower bound in terms of T.'^'^{o) in ()5.106|) . and 
then express the LB outage event in terms of ^d,\ in (|5.105p : 



> E 



E 



Z^;j;(o)>E[h-(] 



''0,0 



E[hoi]rn" 



N 



'0,0 



'0,0 



0,/ 



1 -E 



'0,0 



Si 



0,/ 



9/(0) 



0,/ 



(5.114) 



In what follows we denote h = hi i and h = hi q. The PPP ^d,\ is 
a homogeneous PPP with intensity measure given by Thm. 14.11 with 
A(x) = A and a joint PDF on mark pairs (h, h) independent of x: 



(5.115) 
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The key observation is this: given /io,o, the probabihty that ^d,x is 
empty equals the void probability for x on the set 



Co 

Wo 



{{x, h,h) -.h ^h\x\ " > Wo] 



'-0,0 



(5.116) 



Note in what follows that ho,o is random and hence so is Cq and wq. 
Using Prop. 14.31 and simplifying gives: 



10,0 



>d,A(Co) =0| ho,o) 



(5.117) 



exp<[-A / _ f^{h)f~^{h)dxdhdh 



exp < —A 



h|a;| " > wqI Wq) dx 



where we have defined the RV h = h/h-^. Now echo the development in 
the proof of Prop. starting at (|4.43p . replacing h with h, yielding: 



'0,0 



= expl-AcrfWo-^Eih^]} 
= exp{-AcrfWo^E[hf^o]E[hr,f]} (5.118) 

where in the last step we have exploited the assumed independence of 
h and h. Substituting this last expression into (I5.114|) yields (15.1070 . □ 



Fig. 15.131 shows the asymptotic and LB OP vs. the PCE / (all for 
no noise {N = 0) and Rayleigh fading on all coefficients). We also 
show the no fading asymptotic and LB OP, and observe the OP under 
FPC is always larger than the OP without fading. Although threshold 
scheduling was able to exploit fading and improve the OP/TC to lie 
above that of no fading, this is seen to not be possible under FPC. All 
plots use d = 2, a = 4, ((5 = 1/2), n = 1, r = 5. Further, the asymptotic 
and LB OP are seen to predict the same optimal PCE /* = 1/2 for 
small A, while for larger A the LB optimal FPC is seen to be /* = 0. Fig. 
l5.14l shows the OP LB vs. A, the TC UB vs. q* , as well as the asymptotic 
OP and TC using the same parameters as Fig. 15.131 Similar comments 
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apply: FPC cannot exploit fading to improve performance relative to 
that of no fading, and the (bound) optimal PCE is / = 1/2 for small 
A and / = for large A. 
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Fig. 5.13 The OP q{X) vs. the FPC exponent / for A = 0.01 (left) and A = 0.05 (right). 
The top plots arc the OP LB, the bottom plots arc the asymptotic OP (as q* — > 0). The 
OP is shown for no fading, fading without PC (/ = 0), fading with PC at / = 1/2, and 
fading vs. PCE /. For small A the asymptotic and LB OP agree that the optimal PCE is 
/ = 1/2, while for larger A the LB optimal PCE is / = 0. 
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q(A) q(A) 




Fig. 5.14 The OP <j(A) vs. A, and the TO X{q*) vs. q* . The LB OP and UB TO are on the 
left, while the asymptotic OP and TC are on the right. Curves are shown for no fading, 

and for fading with PCE / = and / = 1/2. The performance without fading is superior 
to performance with fading, for both / = and / = 1/2. The bound on OP and TC is seen 
to be better under / = 1/2 (/ = 0) for small (large) A. 
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Multiple antennas 



6.1 Ml MO with interference 

In this chapter we consider the use of muhiple antennas at the transmit- 
ter (Tx) and/or receiver (Rx). Multiple antenna transmission and re- 
ception, often broadly referred to as "MIMO" (multiple-input multiple- 
output), has been one of the most extensively studied topics in physical 
layer communications over the past fifteen years. It is now a successful 
commercial technology, having strong support in modern cellular stan- 
dards, e.g., LTE and WiMAX, and more recent instantiations of WiFi, 
namely 802.1 In and its successors. Nevertheless, when considering the 
realities of both WiFi and cellular, not to mention ad hoc networks, 
one notices that nearly all theoretical work on MIMO ignores the role 
of interference. 

Understandably, the heralded early results on point-to-point MIMO 
neglected the role of interference or multiple concurrent users. From 
those results we learned that MIMO offered the possibility of not only 
increased reliability through diversity and combining gains, but also 
held out the promise of linear capacity scaling (at high SNR) with the 
number of antennas, assuming they were equally balanced at the Tx 
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and Rx |23t FfTj . This led naturally to considering multiple simultane- 
ous users, for example the downlink of a cellular system where the base 
station has (Tx) antennas, and each user has < (Rx) anten- 
nas. This situation is typically referred to as space division multiple 
access (SDMA) or multiuser MIMO. Here as well, the linear capacity 
(now summed over the users) scaling can still be maintained by using 
a pre-cancellation technique at the Tx (in the downlink), or Rx SIC 

(uplink) [laEHllHllEI]. 

Notably, however, these well-known results all ignore the spurious 
interference that occurs from "uncontrolled" interferers in the network. 
Only a very small body of work has considered the effect of interference 
on MIMO transmission, despite fairly clear warnings that the classic 
results may not hold for moderate to high levels of interference, and in 
fact contradictory results may hold instead [HI [l2l [Til H] • This is prob- 
ably because it appears to be quite difficult to analyze MIMO systems 
(which already have significant randomness to deal with from the ran- 
dom matrix channel) when considering non-Gaussian interference as 
well. The aforementioned cautionary works primarily used simulations, 
as did contemporaneous work on ad hoc networks [62l [SH [I5l [83] . Even 
today, the most promising theoretical MIMO techniques, namely spa- 
tial multiplexing (SM) and SDMA, have proved largely disappointing 
in the field, in large part due to the effect of interference. 

One of the most popular uses of the TC framework has been to 
attempt to better understand MIMO systems that are subject to in- 
terference. The results we present in this chapter are again for the de- 
centralized one- hop ad hoc networks with an Aloha- type MAC, but ex- 
tensions to other examples of interference-limited networks along these 
lines seem well within reach, and are discussed more in §6.51 

6.2 Categorizing MIMO in decentralized networks 

There are two main differences that emerge when considering MIMO 
in a large wireless network, versus an isolated point-to-point or down- 
link/uplink scenario. As we just emphasized, there is the random (non- 
Gaussian) interference aspect, which causes the tradeoffs that emerged 
from traditional MIMO analysis to no longer necessarily hold. Second, 
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there is the spatial reuse aspect to consider, since in a decentrahzed net- 
work the key metric is throughput per unit area. A scheme that gets 
high throughput for a typical Tx-Rx pair at the expense of poor spatial 
reuse may result in a poor TC, or any related network- wide through- 
put metric. Therefore, in this section we distinguish between "single 
stream" and "multi-stream" TC, dividing multiple antenna techniques 
accordingly. Of course, multi-stream techniques reduce to single stream 
as a special case, but for ease of exposition we consider them separately. 




Fig. 6.1 Simple configurations of 1 X 2 receive diversity, 1x2 spatial multiplexing (SM), 
and K = 2 user space division multiple access (SDMA). 



6.2.1 Single stream techniques 

We now provide a general model for a MIMO ad hoc network when 
only a single data stream is transmitted by any given node. Our focus 
is on linear transmitters and receivers but the approach here could be 
generalized to nonlinear systems without any major conceptual differ- 
ences. Among single stream techniques, there are two major competing 
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approaches for how to use the antenna arrays. The first we refer to 
broadly as diversity, where the spatial degrees of freedom (DoF) are 
used to strengthen the desired portion of the received signal. The other 
is interference cancellation (IC), whereby the DoF are used to suppress 
interference. To the best of our knowledge, the first paper on this topic 
was [59] (extended in [iQ]), with [31j shortly thereafter. It has since 
been extended by numerous other works to many scenarios. 

The model used in this section starts with the SINR given in ^4.H 
since these multi-antenna techniques are only relevant when the chan- 
nels have temporal and/or spatial selectivity. In this chapter, we as- 
sume that all channels experience iid Ray leigh fading, i.e., the channels 
are iid complex Gaussian with zero mean and unit variance. Further, 
we assume the dimension d = 2 throughout this chapter, and that there 
is no guard zone, i.e., e = 0. The relevant SINR is still 



and the diversity techniques attempt to increase S(o) (and reduce its 
variance) while IC reduces 51 (o) by removing a portion of one or more 
interfering nodes. 

Definition 6.1. MIMO single stream SINR. For both diversity 
and IC approaches, when linear Tx/Rx filters are used and the number 
of Tx and Rx antennas used per node are and n^, respectively; a 
common model can be adopted for SINR, which is: 



where channel fading is captured by x random matrices Hq for 
the reference Tx-Rx pair, Hj for the channel from interferer i to the 
reference Rx, and filtering is modeled by x 1 Tx filters {vq, vi, . . .} 
and X 1 Rx filters {wq, wi, . . .}. These filters are random in that they 
are in general chosen as functions of the random channels. 

The relevant design question is how to pick the various v and w to 
maximize SINR. For single stream MIMO approaches, the definitions 
of OP and TC developed in the preceding chapters still apply without 
modification. 






2, A 



(6.2) 
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Diversity 

Diversity teciiniques pick v and w to maximize the numerator of (j6.2|) . 
while ignoring the interference. Although some feedback may be nec- 
essary for the design of the Tx filter v, this class of techniques has the 
merit that interference does not need to be estimated, learned, or oth- 
erwise acquired; such acquisition can be expensive in terms of overhead 
and computation, and in many cases {e.g., mobility) not very robust. 

An important point is that in this model, the "diversity" affects the 
SINR in two ways. First, by sending a single stream over multiple ef- 
fectively uncorrelated channels, variations in the SINR are significantly 
reduced and the formerly random SINR hardens to a constant value as 
the antenna arrays grow large. However, the effect of hardening on TC 
is small, since the TC already considers a spatial throughput average. 
Second, the average SINR increases because of array gain, which aligns 
the transmitted energy towards the dominant eigenvalue of the matrix 
channel. This latter effect has a more significant effect on the TC. 

Interference cancellation (IC) 

IC is considered in this chapter for single stream techniques, i.e., we 
design w to be roughly orthogonal to a subset of nearby interferers, e.g., 
the na — 1 strongest interferers. Additional single stream IC gain can in 
theory be achieved if each Tx is able to learn the channels to its — 1 
closest active "victim" receivers, in which case v can be designed to be 
approximately orthogonal to those channels. However, this introduces 
significant overhead and is unlikely to be very robust or practical — it 
also complicates the exposition considerably. For the purposes of this 
chapter, we limit ourselves to Rx IC when considering single stream 
transmission. 

Three different linear Rx IC approaches are considered in the next 
section. Two of them also achieve a diversity gain. The three techniques 
are: 

(1) Zero-forcing (ZF), where w is chosen to be precisely orthogo- 
nal to the channels of the strongest — 1 interferers. There 
is no diversity gain. 
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(2) Partial zero- forcing (PZF), where w is chosen to be orthogo- 
nal to the z < — 1 strongest interferers, and the remaining 

— z dimensions are used for Maximal Ratio Combining 
(MRC). 

(3) Minimum Mean Square Error (MMSE), achieves both IC and 
diversity and maximizes the SINK (and hence TC), but the 
tradeoff between them is not as easy to observe as in PZF. 

Although it is impossible to eliminate all interference with a finite num- 
ber of antennas, for sufficiently large and for low to moderate SINR 
thresholds r, all appreciable interference can in theory be removed 
(since sufficiently weak interferers will not cause outages even in aggre- 
gate). 

6.2.2 Multi-stream models: spatial multiplexing and SDMA 

We would also like to allow the communication of multiple simultane- 
ous information-bearing streams that originate and/or terminate at a 
single user in the network. Although with single stream transmissions, 
multiple Tx-Rx pairs communicate one stream simultaneously, we refer 
here to the case where a single node transmits (or receives) multiple 
streams. Three cases are of interest. 

(1) Spatial Multiplexing (SM): Tx sends K streams to Rx 0. 
All other active transmitters follow suit to their respective 
receivers. 

(2) Tx SDMA: Tx sends K streams total, one each to Rx's 
0, 1, ...,K — 1. For simplicity and maximum contrast to SM, 
we limit our attention to this case, but it is straightforward 
conceptually to extend this to Kj streams sent to Rx i, where 

(3) Rx SDMA: Similarly, the simplest scenario is where Tx's 
0,1, K — 1 send one stream each to Rx 0. 

For such a setup, it is necessary to modify the definitions of OP and TC, 
given in Ch. [TJ Since multiple streams are sent over a matrix channel, 
different SINR statistics will be seen on different streams, motivating 
a per-stream outage constraint ql G (0, 1). 
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Definition 6.2. Multi-stream OP and OCD. Fix a common max- 
imum permissible per-stream OP for each stream = q* £ (0, 1). The 
maximum spatial intensity subject to an OP of q* is the same as in 
the single stream case, but must hold for all K streams: q~^{q*)- For 
the multi-stream case, we call q~^{q*) the optimal contention density 
(OCD). 

For such a definition, one stream will become the limiting factor in the 
OCD, and hence the TC. Because there are now potentially multiple 
streams per transmission, the OCD and TC are no longer interchange- 
able, leading to the following generalized definition for TC. 

Definition 6.3. Multi-stream TC. The TC with K < data 
streams per transmission is the optimal spatial density of concurrent 
multi-stream transmissions q~^{q*) per unit area allowed subject to an 
OP constraint q* , i.e., 

C{q*)=Kq-\q*){l-q*). (6.3) 



This definition generalizes our previous TC definition in Def. 11.31 which 
was simply for K = 1. 

6.3 Single stream MIMO TC results 

We now summarize a subset of instructive results on single stream TC 
with multiple antennas. For the balance of the chapter we neglect noise 
when it makes the exposition more complicated, and include it when 
the exposition is not made significantly more difficult. It is generally 
reasonable to neglect noise since in attempting to maximize TC, we 
are by definition packing transmissions as tightly as possible, render- 
ing thermal noise quite inconsequential compared to the background 
interference level, even after strong interferers are cancelled. 

6.3.1 Diversity 

We first consider diversity techniques, which attempt only to improve 
the desired signal, while ignoring interference. 
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Definition 6.4. Single stream MIMO optimal linear diversity 
filters. In a single stream x MIMO channel Hq, the optimum 
linear diversity Tx filter vq and Rx filter wq are the left and right 
singular vectors corresponding to the maximum eigenvalue of Hq, which 
is denoted as (pQ. 

Because these filters ignore interference, they do not maximize SINR 
or TC. They do however maximize the strength of the desired received 
signal since the optimum eigenmode (pQ is used for transmission. Any 
other Tx filter(s) would not concentrate the entire energy on (pQ, and 
hence have a strictly inferior value of S, the SIR numerator. The main 
challenge to deriving OP and TC for the diversity case, relative to the 
single antenna results of Ch. HI is that now multiple iid channels (both 
desired and interference) are combined at the Rx by the Rx filter. We 
begin with the simpler special case of 1 x MRC, where only the filter 
Wq is required, before preceding to general x eigenbeamforming. 
We adapt our channel notation to the assumed 1 xn^ case in the natural 
way — the channel matrices Hq, Hj are denoted by the n^-vectors ho, hj. 

Theorem 6.1. The Maximal Ratio Combiner (MRC) for a vec- 
tor channel Hq which maximizes the desired received energy is wq = 
hQ/||ho||. For the n^-branch MRC with only Gaussian noise, the post- 
combining SNR is the sum of the per-branch SNRs. 

This well-known result {e.g., [HU]) means that one simply weights each 
branch in proportion to the gain on that branch. Since the interference 
is not Gaussian, the branch SIRs cannot simply be summed. The post- 
combining SIR in our network model is 

Mo) = ^ 1„ i-»iu.u 1-^ = 2- (6-4) 



llhol 



Remark 6.1. Interference distribution is unchanged. Since hg 
and hj are iid and is unit norm, p^^j is simply a linear combina- 
tion of complex Gaussian RVs. As proven in [66], a linear combination 
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of Gaussian RVs is again Gaussian, and due to the normalization of 
ho the mean and variance are maintained as and 1, respectively. 
Therefore, I n^^hj^ is exponentially distributed, and the interference is 
unchanged from the single antenna case in Cor. 14. li 

As a consequence, the SIR can simply be expressed as 

allu, ||2 

sir(o) = ' (6-5) 

2, A 

where ^Ig'^ is the standard aggregate single antenna interference. 

Remark 6.2. Signal distribution. Since u is deterministic and So = 
||ho|p is a KV with 2nj^ DoF, the numerator is distributed as well. 
Clearly for = 1 antenna this reduces to an exponential distribution 
and the result in Cor. l4.1l holds. For rip > 1 the CCDF of So is -Fso(^) = 
e"" Efc=o ' and the associated PDF is fs,{x) = e-"(g^. 

We now generalize §4.11 to 1 x with MRC to give the asymptotic 
(A ^ 0) OP. 



Proposition 6.1. OP vi^ith MRC gO]. The OP for 1 x MRC is 



k=l 1=0 



where k = \T^v?Ca and = 7r^(5 csc(7r5). 



Proof. The proof requires significant algebra and is given in |40j . A 
brief version with key steps is as follows. We begin as in the proof of 
Prop. 14. 2i Denote ^Ig'^ by 51 for this proof with PDF f-£{t). Then 

l-q{X) = P(So>™"Z)= / Fs,{st)Mt)dt 

Jo 

r°° / 1 \ Id*' 

(6.7) 
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where the second equaUty in (16.71) uses the LT property t"^ f{t) < — > 
(-l)"£^/:[/(i)](s) and s = tu°'. For Rayleigh fading, as in gHHH 
/^[5I](s) = ex.p{—XCaS^)- Forming the first order Taylor expansion for 
the jth. derivative around k = Xs^Ca = 0, note that any term with 
for A; > 1 is o(k) and can be discarded. Eventuahy it can be shown that 

l{-sf^£[ms) = -l^Xs'C^ll{l-6) + Q{^') (6.8) 

as A — >• 0, which after straightforward manipulation, yields the result. 

□ 



Note that Prop. 16. H is equivalent to Cor. 14.11 for = 1, no noise, and 
given a Taylor series expansion around n. 



Proposition 6.2. TC with MRC [10]. When each Tx transmits on 
a single antenna and each Rx performs MRC with antennas; or 
equivalently each Tx performs MRT with = antennas and each 
Rx uses a single antenna; the asymptotic (as q* — )• 0) TC under iid 
Rayleigh fading for all links, with no noise, is 

X{q*) = . ^ ^ + e{q* f, q* ^ 0. (6.9) 



Further, X{q*) is Q{nji^) (as — )• oo) and ignoring Q{q*)^ terms can 
be bounded as 

l<%^^A(,*)<r(l-5). (6.10) 



The "exact" result (|6.9p follows immediately by solving the result of 
Prop. 16.11 for A. The bounds are more involved to show, but it can be 
shown that ^1 + Yl]!^=i^ H nf=o^(^ ~ be bounded by nR^T{l—S) 

(lower) and n^^ (upper) [3D]. The analysis for Tx antennas and 1 Rx 
antenna follows precisely the same development. However, we do not 
account for the cost of acquiring (non-causally) the necessary channel 
state information hg at the Tx. 
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Remark 6.3. Sublinear gain of MRC in tir. The TOr antennas 

A — 

utilized in MRC provides a gain of = whicli is sublinear since 
a > 2. The same scaling holds if noise is included [20]. Interestingly, the 
gain is larger for small path loss exponents which may seem counter- 
intuitive since it would seem that diversity is increasingly desirable 
when the signal propagation is poor. 

The results on MRC can be extended to the n-r x tt-r eigenbeamforming, 
where the Tx uses what is sometimes called Maximal Ratio Transmis- 
sion (MRT), and the Rx uses MRC. The Tx and Rx filters simply put 
all energy onto the dominant eigenvalue of the channel matrix ho, i.e., 
vq and Wo are set equal to the input and output singular vectors of ho 
corresponding to the maximum singular value of ho- 

Proposition 6.3. OCD of x eigenbeamforming [IQ]. The 

asymptotic (as q* — )• 0) OCD can be bounded (neglecting @{q*)'^ terms) 
as 

(max{nT,nR})^g* ^ „-l^„*^ ^ Tjl - S){nrrnnYq* 



Remark 6.4. TC scaling with n-^ x tir eigenbeamforming. The 

scaling of max{nT,?T.R} observed in the LB follows immediately from 
Prop. 16.21 since simply using MRT at the Tx with tir = 1 gives n-Y^ 
scaling, and similarly using MRC at the Rx with = 1 gives tt-r"^ . Thus 
the scaling with n-^ x tir must be at least the larger of those two, since 
the smaller array must be at least as good as using a single antenna. 
The UB indicates the possibility of superlinear TC scaling for tir = tit 
and a < 4. Although a tighter UB is not available, we conjecture based 
on other results in random matrix theory and simulations that the LB 
is more accurate and the scaling is max{?7-T,nR}. 



6.3.2 Interference cancellation (IC) 

One can instead design the Rx beamforming filter to suppress inter- 
ference from a subset of the nearby interfering nodes, rather than for 
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enhancing the desired signal strength. As in the last section, we as- 
sume TT-T = 1. To illuminate the tradeoff between these two opposing 
beamforming design philosophies, we adopt a general but suboptimal 
structure for the Rx beamformer wq which we term partial zero-forcing 
(PZF). We briefly visit the optimum but less illustrative MMSE Rx 
towards the end of this section. 

Definition 6.5. Partial zero forcing (PZF) Rx. The PZF-z Rx 
uses a beamforming vector wq that is orthogonal to the channel vectors 
of the z strongest interferers, where z < — 1. That is, hj _L wq i = 
1,2, ... z. Furthermore, ||wop = 1 and the tir — z remaining DoF are 
used to maximize the desired received power. Formally, if the columns of 
the X (riR— z) matrix Q form an orthonormal basis for the nullspace of 
(hi, . . . , h^), which can be found by performing a full QR decomposition 
of matrix [hi • • • h^], then the Rx filter is chosen as: 

Wo = 77——; [7- (6.12) 



Note that if z = 0, Q = I and wq is the MRC beamformer of Thm. 
16.11 If z = tih — 1 then we have conventional ("full") zero forcing (ZF); 
note that "full" ZF is only on the closest z interferers and the rest of 
the interferers are treated as background noise. 

Proposition 6.4. PZF SINR. The SINR for PZF-z is 

s^^^''^'-'io) = -—^ 1 ^ _i (6.13) 

where So is X2(nR-2)' = l^o'^jP terms are iid unit-mean exponen- 

tial RVs and also independent of So, the quantities jx^+ip, |x^+2p, . . . 
are the z + 1, z + 2, . . . ordered points of a 1-dim. PPP (i.e., the closest) 
with intensity vrA, and the ordered points are independent of the signal 
and interference terms. Finally, snr = Pu^"/N as in Def. 



We note that Hi, . . . , terms are not in the expression because those 
interfering nodes have been (perfectly) cancelled. Unlike in the MRC 
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case, we will need to concern ourselves with the interference distribu- 
tion, and so define the aggregate interference power for PZF-z as: 

oo 

Z, = n" J2 l>^ir"Hi. (6.14) 

i=z+l 

To understand how PZF-z (and the choice of z) affects the TC, it is 
sufficient to consider upper and lower bounds, since we shall see that 
the scaling in is the same for both bounds, and also for the MMSE 
Rx. 

Proposition 6.5. OP for PZF-z [16]. The OP with PZF-z has UB: 

r((7rn2A)^(f -l)-'(z- [1])'-^ + ^) 



for I"!] < z < - 1. 



Proof. First, rewrite the OP as the tail probability of the RV 1/sinr 
and then apply Markov's inequality as follows: 



pzf— z 



(A) 



sinr 



1 ^ 1 

pzf — Z — J- 



\ (a) 


1 


I < T-E 


.sinrP^f-^ 





(6) 









" 1 " 


T-E 




E 






snr 




% 



(6.16) 



where (a) is due to Markov's inequality, (b) is from ()6.13p and the 
independence of So and Hj. We note that the first expectation term 
Hz + —1 = E[5I^1 + — and corresponds to the effect of IC, whereas 

^ snr J I- snr ^ ' 

^ = — V~r since So is yL„ ,w and this term corresponds to the 
signal power boost from the remaining DoF. The remaining task is to 
find E[5l2], and an upper bound on it (suppressing the leading term) 
can be found as 



E 
E 



E 



.i=z+l 



oo 

E [|x, 

i=z+l 



(6.17) 



1=2 + 1 



from the independence of |xj| and Hj and unit mean fading. Because 
l^iP) |x2pj ■ ■ ■ are a 1-dim. PPP with intensity vrA, RV 7rA|xjp is xii 
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and thus has PDF /x(x) = ^-(jiiiy 
e[(|x,|2)-"/' 



see Rem. 16. 2|) . Therefore, 



T i) 



(6.18) 



This quantity is finite only for i > ^, and thus the expected power 
from the nearest uncancelled interferer is finite only if z + 1 > ^ • As in 
[37] . using Kershaw's inequality, 



r(.-f) 



T(i) 



< i 



.19) 



where [•] is the ceiling function and i > . Therefore 



i=z+l 



E 

i=z+l 



< 



2 
2 



^-1 
2 



-1 



dx 



(6.20) 



where the inequality in the second line holds because x~2 is a decreas- 
ing function. Inserting this expression for ^[Z^] along with E 
(|6.16p gives the desired result. 



into 
□ 



This result for OP is easy to invert for A which yields our next propo- 
sition, and a useful interpretation of it. 



Proposition 6.6. TC for PZF LB. The TC of PZF is lower bounded 
by: 

(6.21) 

for any z satisfying \^~\ < z < — 1 — —r^- Furthermore, if 2; = 
for < ^ < 1, then X{q*)^^^^^ = 0{n^) (as — )■ 00), and the optimum 
value of 6* is 61* = 1 - ^. 
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The linear scaling of PZF can be observed by simply plugging in z = 
Oriji and allowing nj^ to grow large. Similarly, the optimum value of 
can be found taking the derivative of (|6.21|) w.r.t. 6 after plugging in 
z = 6n^, setting equal to zero, and solving for 6. 

Remark 6.5. PZF linear scaling in n^. The linear scaling can be 
interpreted as adaptively combining the MRC portion of the array (to 
get n^^) with the IC portion (to get ), for a total scaling of n^. 

In particular, the amount used for each must grow with n^, any fixed 
value of z does not achieve linear scaling in n^ . 

To be certain that PZF cannot achieve superlinear TC scaling in n^, 
we turn our attention to appropriate upper bounds. One can consider 
two different approaches. The first and most direct is simply to lower 
bound the PZF OP and invert it to observe a TC UB. The second is 
to consider the MMSE Rx, which is by definition strictly better than 
PZF, and hence upper bounds its TC. For completeness, and because 
the proof techniques are very similar for each, we provide both upper 
bounds. First, however, we must define the MMSE Rx beamformer. 

Definition 6.6. The MMSE Rx filter is given as 

where A is the random spatial covariance of the interference plus noise 
and can be expressed as 



A=— I + n° y Ix^r^h^h*, (6.23) 
snr ^-^ 

and I is the x identity matrix. 



Remark 6.6. MMSE Rx filter. The MMSE filter maximizes SINR, 
which is sinr = HqA^^Ho, but is not easy to express in a way amenable 
to analysis. Furthermore, it should be noted that although it appears 
superficially from ()6.23p that a great deal of information (namely the 
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distances and fading values of every single interferer) is needed to com- 
pute A, in fact this covariance matrix can be estimated over time and 
is generally more robust to imperfections than the ZF matrix compu- 
tation, which requires the BF vector wq to be precisely orthogonal to 
each of the z interfering channels. 

With this background on the MMSE filter (see [SH HH [Ml E] for more 
detailed discussion) , we can now state the two TC upper bounds in the 
following proposition. 

Proposition 6.7. TC UB for PZF and MMSE [46]. The TC for 
an MMSE Rx with antennas at high SNR is upper bounded by 

A— (g*) < + (6.24) 

while for PZF the corresponding bound with I uncancelled interferers 
is 

TTU^Ta[l — q*)a V ' ^ / 

The PZF upper bound holds for any < z < — 1 and any integer 
/ > 2. These bounds both scale as ©(n^) (as — )• oo). 

Proof. The complete proof of the MMSE TC result is given in [46] . 
The first step is showing that an UB on success probability (neglecting 
noise) can be written as 

1 _ ^mmse(^) < p [ ^I^M^^ > \ ^ (g 36) 

which follows from an OP lower bound for MMSE receivers developed 
in [27] . Using Markov's inequality and a similar approach to the proof 
of Prop. 16.61 - namely using Kershaw's inequality, dropping various 
terms (while preserving the bound), and exploiting the independence 
of various terms - the desired result can be achieved. The proof of 
the PZF UB follows the same method, starting with the PZF SINK 
expression. Both results scale linearly with tt-r which can be observed 
by inspection. □ 
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This proposition, paired with Prop. 16.6^ estabhshes that by appropri- 
ately combining diversity (array) gains with IC, the TC of a wireless 
network can in fact be increased linearly with just the Rx antenna ar- 
ray size. In contrast, this is not possible without the combination. This 
can also be observed from the above upper bound on PZF TC, which 
reduces to MRC for z = and to full ZF for z = n^ — 1. Plugging those 
values in yields the following Corollary. 

Corollary 6.1. TC UB for MRC. An upper bound on the TC with 
MRC is given by 2: = and is 

KQ*r'' < (6.27) 
ttu^tq (1 — q*) a 

while for full ZF the TC is giving with z = n^ — 1 and 

MgYin. - 1) < (6.28) 

TTW^Ta (1 — q*)a 

Both of these expressions use / = 2 in Prop. 16.71 

These expressions make plain the n^^ scaling for MRC and ^ 
scaling for full ZF, respectively. These same scalings were observed 
above for MRC in Prop. [6^2] and for full ZF in [37j, respectively. Notably, 
in [37], it was found that using multiple Tx antennas did not change 
the scaling if the Tx beamforming vector is not adapted to the channel 
or interference. 

6.4 Main results on multiple stream TC 

We now consider the case where nodes communicate K > 1 streams 
simultaneously. In the case of spatial multiplexing, these K streams all 
originate and terminate between a single Tx-Rx pair. For (Tx) SDMA, 
the streams all originate at a single Tx but terminate at K different 
receivers. 

6.4.1 Spatial multiplexing 

The fundamental tradeoff in a spatial multiplexing (SM) system is one 
of rate versus reliability. Famously formalized in , the basic idea is 
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Fig. 6.2 The OCD vs. ?1r for PZF, MMSE, ZF, and MRC, for a = 4. The superior scaling 
of MMSE and PZF is readily observed. 



that transmitting multiple streams, e.g., K = n^, tends to maximize 
the data rate while K = \ tends to maximize reliability. We will see that 
a related tradeoff also occurs in decentralized wireless networks follow- 
ing the TC model. In our case the tradeoff is typically more complex to 
quantify because the TC of each Tx-Rx pair is affected by all the other 
Tx-Rx pairs in the network, who also are adjusting their number of Tx 
streams. However, as a metric the TC already inherently includes both 
throughput and reliability, so maximizing the TC inherently trades off 
between rate (higher SINR or A) and reliability (outage). 

As noted previously, the TC framework operates on a specific sys- 
tem model, which renders a SINR expression specific to that model. 
In the case of SM there is a fairly large library of possible Tx and Rx 
structures, each of which yield a different SINR, and hence possibly a 
different tradeoff. For example, a SM Tx could use multimode beam- 
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forming, a singular value decomposition (SVD) pre-filter, or it could 
simply transmit each of the K streams on K < rir^ different antennas 
— either in a fixed pattern or rotating over the antennas (as is done 
in BLAST). At the Rx, the number of options is even larger, including 
the Maximum Likelihood (ML) detector or approximations to it like 
the sphere decoder; various forms of linear interference cancellation 
(including ZF, PZF, and MMSE); nonlinear interference cancellation 
(BLAST or SIC); or other approaches like MRC or SVD post-filtering. 
Explaining the details of all these different well-known Rx structures is 
outside the scope of this chapter, and readers lacking this background 
are referred to any of (BOj [72l [30l [28] for details. We will focus on an 
instructive subset of these approaches (mainly limited to Tx and Rx 
very similar to those considered in ^6.3p . and see that several broad 
trends, as well as a few interesting differences, hold across the different 
approaches. 

Simple i^-stream transmission 

First, we consider a class of results where the Tx has the simplest 
possible structure: it simply transmits a single unique data stream from 
each of one or more of the Tx antennas. So, for K = 1 these results 
reduce to the results of §6.31 For 2 < K < min(nT, ?1r), the results will 
generalize our prior results. For this simple transmit structure, we will 
consider four possible different receivers, all of which are linear: 

(1) The MRC Rx [56]. 

(2) The ZF Rx, where only self-interference from the interfering 
K — 1 streams is cancelled [681 [56] . 

(3) The PZF Rx, where in general both self-interference and in- 
terference from other users is cancelled |74j . 

(4) The MMSE Rx [55]. 

We now catalogue and discuss the main results of interest for these 
cases. The proofs are not given but can be found in the referenced 
papers. To express the results simply, we assume in all cases that n-r < 
and that the number of transmitted streams K < n-Y, and so of 
course K* < also. Due to the transmit structure assumed in this 



section, the other n. 
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— K Tx antennas are simply not used. 



Proposition 6.8. SM with MRC [56]. With a MRC Rx and large 
rirYjUji, the optimal number of streams to transmit is 



1 2 



K*=n^, (6.29) 
r(l + snr 

and when K = K* the resulting TC is 0(71^) (as — t- oo). 



Proposition 6.9. SM with ZF Rx [MIES]. With a ZF Rx applied to 
the K — 1 interfering streams, other-user interference treated as noise, 
and large nxj^a, the optimal number of streams to transmit is 

1-1 

^*=^RT- (6-30) 

1 + Tsnr ^ 

and when K = K* the resulting TC is 0(?1r) (as — ^ oo). Further- 
more, if K = the scaling falls to 0(nR/ q) (as — oo). 

Several interesting observations can be made already from these two 
results. At high SNR, the optimum number of streams is quite different 
for MRC and ZF, namely ZF uses r more streams than MRC. They 
both prefer more streams at high path loss exponents: in this case the 
other-user interference is attenuated more rapidly. And we see that with 
ZF, there is a significant penalty attached to sending streams from the 
full antenna array. 

Proposition 6.10. SM with PZF at high SNR [74J. When a PZF 
Rx is applied to cancel the K — 1 interfering streams from the same 
user, and the K streams from each of z closest interfering users, the 
optimal number of streams is K* = 1, the optimal number of users to 
cancel is z* = n^{l — §), and the TC is ©(n^) (as — oo). 



This proposition indicates that the optimal PZF setup is the single 
stream Tx of [36] that was already extensively discussed in §6.3.21 In 
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short, there is in principle no gain to doing SM if one is able to cancel 
interference instead. An important caveat here is that one must be able 
to increase the contention density arbitrarily to get the linear scaling 
gains {e.g., by allowing a higher contention probability in the Aloha-like 
random access protocol). SM may therefore still be desirable since it 
does not require an increased network density to achieve its TC gains. 

Corollary 6.2. PZF Rx without cancellation [74J. If the PZF Rx 
operates only on the streams of the desired user, i.e., z = 0, then 
K* = nR(l — ^) and the TC is still 0(71-^) (as — )• oc). 



This corollary follows from Prop. [6^91 and from [71], because the PZF Rx 
is in this case exactly the ZF Rx. The TC scaling is the same whether 
the Rx cancels interference from the other streams of the same Tx; 
or cancels interference from a number of nearby single stream trans- 
mitters [56]. Furthermore, the optimum number of streams/users to 
cancel turns out to be the same, which is not obvious since the sta- 
tistical properties of the interference are quite different in each case. 



Corollary 6.3. PZF Rx cancelling strongest interferer [71]. If 

the PZF Rx instead cancels the strongest interferers, then K* = 1, 

z* = na — 1, and the TC is still @{n^) (as n^, — )■ oo) but can be further 

1 

characterized as @{n^{q*)"R) (as — )• oo). 

This tells us that there is an additional gain in terms of the outage 
constraint from measuring and cancelling the strongest users, which 
would typically also be more practical since the Rx would not know 
where the interferers are located, but would be able to measure the 
receive signal strengths. 

Remark 6.7. SM with the MMSE Rx. Exact results for the TC 
in this case are given in |55j but are quite complex, and not amenable 
to a simple conclusion on the orderwise scaling or optimum number 
of streams K. In view of the previous results it is safe to conjecture 
that the scaling will be at least linear in if the optimum number of 
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streams are used. In fact from Prop. 16.71 we know that the scahng is 
hnear even K = 1. 



Spatial multiplexing enhancements: multimode BF and 
DBLAST 

The prior section used the simplest possible Tx structure, and four 
different linear receivers. In this section we consider two other repre- 
sentative setups that introduce logical enhancements at both the Tx 
and Rx. These are: 

(1) Multimode (eigen)-beamforming transmission [74j . 

(2) DBLAST, which includes a diversity Tx and non-linear in- 
terference cancelling Rx |68j . 

The multimode beamformer transmits K streams on all Ut antennas, 
where the K streams are placed on the K dominant eigenmodes of 
the Tx-Rx channel. For simplicity we now let and so for 

the balance of the section and can be used interchangeably. 
Obviously, the Tx must acquire the Tx-Rx channel to implement such 
a precoder, which typically requires the singular value decomposition 
(SVD) of the matrix channel. Such an approach achieves the capacity 
(with appropriate power allocation across the eigenmodes) for a point- 
to-point MIMO channel with an optimal Rx [71J. This approach was 
the first considered approach in the TC framework for ad hoc networks 
as well |38) and is difficult to analyze; Vaze and Heath made progress 
by instead using a PZF Rx. Their main result can be summarized by 
the following proposition. 



Proposition 6.11. SM with multimode beamforming and PZF 

receivers [74j. For a multimode beamformer with a PZF Rx, the opti- 
mum number of Tx streams is K* = 1, the optimal number of cancelled 
interferers is z = — 1, and the TC scales as Q{n^) (as — )• oo). 



In short, applying an optimum Tx pre-filter does not change the PZF 
results we have previously seen; however it does increase the optimum 
number of cancelled interferers since the additional Rx DoF are no 
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longer needed for diversity and array gain, since that is accomplished 
at the Tx. There is a fixed (independent of nT,n^) gain of about a 
factor of 4 with multimode beamforming. With an optimum single- 
user Rx (computed on the SVD of the Tx-Rx channel), but no IC, it 
was observed numerically in [38] that the optimal number of streams 
reverts to be compatible with the recurring nT(l — ^) expression. 

Turning to nonlinear Rx structures, Stamatiou et al. considered 
both Vertical and Diagonal BLAST architectures in [68j. BLAST re- 
ceivers successively decode and cancel the K transmitted streams. V- 
BLAST is more bandwidth efficient but does not achieve as much di- 
versity as D-BLAST, which rotates the symbols across the antennas. 



Proposition 6.12. SM with the BLAST architecture [M]. For 

a D-BLAST architecture the optimum number of streams is K* = 
2(nT -M)(l - I) and for both D-BLAST and V-BLAST the TC is stih 
©(tt-t) (as n-Y — )• oo). Furthermore, the TC of V-BLAST is higher than 
D-BLAST by a factor of 2^'^ for a < 4 and by a factor 2-^5-\l-5f-^ 
for a > 4, recalling 5 = —. 



It is notable that BLAST appears to make SM more robust, which 
results in more streams being used especially for high path loss (lower 
interference). For example, if a > 4, K* = tt-t. 

6.4.2 Space division multiple access (SDMA) 

The main distinction between SDMA, also often called multiuser 
MIMO, and spatial multiplexing is that in SDMA, the multiple si- 
multaneous streams originate or terminate at different users. We will 
focus on the Tx SDMA case, and to maximally distinguish between 
SDMA and SM, we will assume that each SDMA stream is sent to a 
different Rx. We will now formalize key aspects of the model, before 
providing key TC results on SDMA, and discussing their implications, 
and how they differ and/or agree with cellular SDMA results. 
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Each Tx has antennas and communicates simultaneously with |/C| = 
K < Ut receivers, which are each equipped with Rx antennas. Each 
of the K streams contains a separate message destined to a different 
Rx, and each Tx and its set of intended receivers /Cj form a broadcast 
cluster. The closest Rx in the cluster is a distance Umin from the Tx, 
and the farthest is Umax, with the rest of the receivers proportionally 
located in a random direction in between Umin and Umax- Although Umin 
and Umax show up in the TC bounds, they can take on any arbitrary 
value and do not affect the scaling results we will present. 

Dirty paper coding (DPC) transmission 

It is well-known that the optimum Tx SDMA strategy in a Gaus- 
sian broadcast channel {i.e., a downlink channel with only Gaussian 
noise/interference) is dirty paper coding (DPC) |20j . which can be 
viewed as a form of successive interference pre-cancellation, and pro- 
vides an effectively interference free channel for each stream |13^ \7E[ IST] . 
Therefore, the SDMA sum capacity (achieved using DPC) with single 
antenna receivers at high SNR is ©(n^ logsnr), and adding Rx an- 
tennas to each node does not increase the scaling by more than a fairly 
inconsequential log log term, which can be achieved with antenna 
selection (which increases the snr by a lognjj factor asymptotically, 
and hence the capacity by log log n^). However, until recently, almost 
nothing was known about SDMA's performance in the presence of un- 
controlled (spurious) interference. 

The next several results, due to Kountouris |52l[5l], provide math- 
ematical expressions that show SDMA's effect on throughput in a de- 
centralized network. The proofs can be found in the sources referenced 
when not provided here. 



Proposition 6.13. Bounds on multistream DPC TC with 

MRC The muhi-stream TC of DPC transmission with MRC 

receivers at high SNR subject to small OP constraint q* is bounded 
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as 



Kq*{l-q*)Td 



^ pTOXC ^ 

^ L-DPC — 



K[4nT(nR,-K + l)]^(l-g*)log 



1 



) 




(6.31) 



where for diversity order d = nT(^^R — K + 1) we have 



1 -1 




(6.32) 



and the constant Jk is 



Jk 




) 



(6.33) 



which depends only on the number of streams K and a. 

Proof. A sketch of the proof is provided, which follows the standard 
TC framework. First, derive an upper and lower bound on the OP, 
which results in bounds that are given by 



where £[5I](s) is the LT of the aggregate interference term with multi- 
stream transmission (given interference power marks that are dis- 
tributed as chi-squared with 2K DoF), and can be given by |38j 




The closest Rx in the cluster gives the lowest OP (almost surely in 
probability), and the farthest similarly gives the highest. These expres- 
sions are then set equal to q* and inverted for A to find corresponding 
lower and upper bounds on the contention density X{q*), which mul- 
tiplied by the the number of transmitted streams K and the success 
probability I — q* gives the result. □ 

These somewhat involved expressions are amenable to the following 
scaling laws, which clearly show the dependence on the number of 
streams and antennas. 






(6.35) 
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Proposition 6.14. SDMA scaling laws with a DPC Tx and 
MRC Rx 



CEpc = ^{K^-^iinr- K + l){n^,- K + (6.36) 
Cg'fc = 0{K'^-i[n^{n^, - K + (6.37) 

where both asymptotic order expressions hold as r^T,?^R — oo. 

Proof. Given Prop. 16.131 for asymptotically large number of antennas 
TT-T, n-R, and streams K, as in the proof Prop. 16. 2| that 

lim ^ = ^r(l - 2/a) (6.38) 



and for asymptotically large d we have Td = S(d") [SOj. Thus, the 
lower bound divided by [(ut — K + l){n^, — K + 1)]^ converges to 

a constant as n^, — )■ oo, and similarly for the upper bound divided 
hy K^^^{n^{n^ - K + 1))^ . □ 



Corollary 6.4. TC scaling and K* for SDMA with DPC Tx 
and MRC Rx [52|. Letting and the number of streams 

K = ^tit, for < 6* < 1 in Prop. 16. 14^ then 

CJ^f(. = S7(nT), Cf^^-c = 0(nT^+-), riT ^ oo (6.39) 
and it can be shown that 

^ib = = 1 - — Att^t, K^ = l--{nT + l). (6.40) 



Considering just the UB, we observe that the optimum number of 
streams ^C*^^ again follows the same theme as the PZF and SM re- 
sults, with a similar but slightly different trend for K*,^. But unique 
thus far to SDMA, we see the possibility of superlinear scaling in the 
number of antennas. 

If we consider a trivial Rx that has a single Rx antenna, the following 
scaling law can be determined for such a multi-input single output 
(MISO) scenario. 
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Proposition 6.15. TC with DPC Tx and single antenna Rx 
[52]. At low q* and high SNR, The TC with a DPC precoder and single 
antenna receivers scales as 

CJJig = n{K^~i{nT, -K + l)i) and Cg-^fg = 0{K^-inJ), 

(6.41) 

(as Ut — )• oo), and ii K = On^ for any 9 G (0, 1), then Cg^pg = 0(nT). 
Finally, the asymptotically optimum number of streams is K* = (1 — 

— I'T-T- 

a' ^ 



Prop. 16.151 — where K streams are transmitted to K different single 
antenna receivers in space — forms an agreeable symmetry with Prop. 
I6.6l where z streams are transmitted by z different single antenna trans- 
mitters to z different riR-antenna receivers in space. In both cases the 
scaling laws are linear in nx and tt-r, respectively, and optimum number 
of streams is the fraction 1 — ^. Furthermore in both cases, fixing the 
number of streams to some constant that does not depend on n-r of 
loses the linear scaling. 

6.5 Practical issues and further research 

In this chapter, we have considered numerous ways to use multiple 
antennas at the Tx, Rx, or both. The results provided here follow a 
few broad themes, which we now briefly summarize. 

6.5.1 Summary of main design insights 

Diversity techniques typically provide linear SINK gains coming from 
the array gain of the antenna array. This then results in a TC gain of 

2_ 

La m most cases, where L is the amount of diversity, e.g., L = nxn^. 
For example, for a pathloss exponent of 4, the TC scales as \/L- This 
is considerably better than in a hypothetical centralized system where 
interference does not exist, in which case the throughput scales about 
as log(L) for diversity. 

Linear capacity increases. We observe that there are in effect three 
ways to increase the number of streams per unit area. The first is to 
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increase the density of independent transmissions and tolerate higher 
interference, e.g., using Rx IC. The second is through SM between a 
given Tx and Rx. And the third is via SDMA. Interestingly, we see 
that all three of these can linearly increase the TC with the number of 
antennas. Rx IC and SDMA require only the Rx or Tx, respectively, to 
have multiple antennas, whereas SM requires both the Tx and Rx to 
have multiple antennas. However, it may not be possible to simply in- 
crease the density (Rx interference cancellation) and the SDMA result 
requires perfect CSIT and the feedback for this has not been accounted 
for. 

Balancing diversity and multiplexing. An interesting recurring 
theme in many of the results is that the number of streams sent should 
not be sent to the maximum, but rather a fraction 1 — — , with the 
remaining DoF going to diversity. For example, with a = 4, equal 
antenna resources should be devoted to diversity and spatial pack- 
ing/multistream transmission. No such corresponding result exists, to 
the best of our knowledge, in classical MIMO theory (without spurious 
interference) . 

6.5.2 Caveats and practical issues 

There were numerous simplifications made in this chapter that will be 
easily recognized by those familiar with MIMO systems and/or ad hoc 
network design. We comment on a few of the significant ones. 
Random access MAC. Throughout this chapter, indeed for any re- 
sults using the PPP as a spatial interference model, the implicit MAC is 
an uncoordinated slotted Aloha- type MAC. Such a MAC is pessimistic 
versus CSMA (which is quite practical) or a centralized MAC (which is 
likely not practical for an ad hoc or self- forming network) . Moving away 
from an Aloha MAC (which gives worst-case interference) may change 
some of the relative benefits and make SDMA/SM more desirable ver- 
sus interference suppression [86]. Some initial progress on extending 
the MIMO results to a CSMA setup are given in [IT] . 
Scaling results are focused on in this chapter, to see how the TC slope 
scales for large nrY,n^. In most practical scenarios, regime of small n-j- 
and is the most important one, and in such a regime, the use of the 
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antennas to decrease OP may be more important from a TC point of 
view. For example, [37j indicates that for small all DoF should be 
used for IC, in which case the OP falls as A"^^^, or equivalently the 
TC can tolerate a stricter outage constraint at the same A. Then, as 
the array size grows, some DoF can be converted to use for MRC or 
sending multiple streams. 

Cellular systems do not directly follow from the results given here, 
which are for decentralized networks. The main modeling difference 
in a (downlink) cellular system is that interferers should not be 
closer /stronger than the desired signal, or else one would simply hand- 
off to this stronger signal. A tractable framework for cellular analysis 
which uses tools similar to those in this monograph was recently pro- 
posed in [3] , where the base stations are drawn from a PPP, and hence 
the interference has a shot noise format and many of the results in this 
monograph can be modified and applied. We conjecture that many of 
the key results and scalings shown here may hold in a MIMO-enhanced 
cellular network, but this is an interesting topic for future work. 
Feedback and overhead. We have neglected the cost of obtaining 
the channel state information at the Rx (e.g., using pilot symbols) or 
more importantly for SDMA and possibly SM, at the Tx. SDMA in 
particularly requires frequent feedback. However, [3H E3] have consid- 
ered exactly this problem, and [32\ [53] has shown that the feedback re- 
quirement for preserving the SDMA linear capacity scaling is log(T), 
which follows the well-known broadcast channel result in [H] , but inter- 
estingly includes the effect of interference and thus uses quite different 
analytical techniques. 
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Table 1 Notation 



Chapter 1 



c(o) 


(random) channel capacity for Rx at o 


Ass.fm 


R 


codebook rate 


Def.lTTI 


T 


SINR threshold required to support rate R 


Def.O 


9(A) 


outage probability (OP) 


Dcf.O 


d 


dimension of the wireless network ({1, 2, 3}) 


Def.O 


A 


intensity of attempted transmissions per unit area 


Def.O 




(random, homogeneous) PPP on R"* of intensity A 


Def.O 


Q* 


target OP 


Def.O 


Hi') 


transmission capacity with target OP q* 


Def.O 


Apot 


spatial intensity of potential interferers 


Rem.lTTI 


Ptx 


Aloha transmission probability 


Rcm.fm 


Chapter 2 


bd(c, r) 


ball in K'' centered at c with radius r 


Def. 12.11 


ad(c, ri,r2) 


annulus in M'' centered at c with radius ri, 


Def.O 




coefHicent on volume of a ball: |bd(o, r)| = c^t'^ 


Prop. O 




Gamma and incomplete Gamma functions 


Def.O 


Fx, Fx, fx 


CDF, CCDF, PDF for RV x 


Def.O 


X ~ N(p, a) 


a normal RV with mean ^ and standard dev. a 


Def.O 


z ~ N{0, 1) = Fz 


standard normal RV 


Def.O 




inverse CDF and CCDF for RV x 


Def.O 


£[x](s) 


Laplace transform (LT) for x with s £ C 


Def.O 


?i[x](t) 


characteristic function (CF) for x with t £ R 


Def.O 


>![x](0) 


moment generating function (MGF) for x with d £ R-|- 


Def.O 


W[x](x) 


hazard rate function (HRF) for x with a; S R 


Def.O 


«(r) 


a generic impulse response function 


Def.O 




(random, cumulative) shot noise (SN) RV 


Def.O 




(random) max shot noise (SN) RV 


Def.O 


a 


pathloss exponent 


Rem. [Til 


e 


nulling radius around Rx 


Rem.O 




pathloss attenuation function over distance r £ 


Rem.O 




(random) SN at o under PPP 11^ \ and la^t 


Def.O 


a 


characteristic exponent 


Def.O 


7 


dispersion coefficient of a stable RV/CDF 


Def.O 


g[U,u] 


point processs probability generating functional 


Dcf.l2Jil| 


Chapter 3 


sinr(o) 


(random) SINR at o 


Def.O 


P 


transmission power 


Def.O 


N 


noise power 


Def.O 


S 


Rx signal power 


Def.O 


u 


Tx-Rx pair separation distance 


Def.O 




a constant 


Def.O 


snr 


Rx SNR 


Def.O 




(random) dominant interferers and interference 


Def.O 




(random) non-dominant interferers and interference 


Def.O 


A{A) 


spatial throughput (TP) 


Def.O 
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Table 2 Notation (continued) 



Chapter 4 





(random, cumulative) interference at o under fading (h) 


Def.liH 




(random) channel fading coef. from i £ XI^j x to Rx at o 


Rem.liAl 


ho 


(random) channel fading from ref. Tx to ref. Rx at o 


Rem.O 


X 


(random) marked PPP 


Rem. 112] 


£0,9(0), g(0) 


event/probability of outage w/o interference 


Rem. 113] 


*^.A,^^:i(o) 


(random) dominant interferers and interference 


Def.lO 




(random) non-dominant interferers and interference 


Def.lO 




(random) link distance under VLD 


Def.lO 


9(A), A(g*) 


OP and TC under FLD 


Prop. HJ] 


g(A),A(g-) 


OP and TC under VLD 


Prop. Hj] 


[/ 


end-to-end (source-destination) distance 




M 


number of hops 


S3] 


A 


maximum allowable end-to-end transmission attempts 


Def.l431 


T 


(random) total number of transmission attempts 


Def.l431 


Amh 


multihop transmission capacity 


Def.[431 


\ ub 


upper bound on Aj^ij 


Prop. |4J2] 




7r^(5 csc(7r(5) 


Prop. 14.131 


Chapter 5 


W 


total bandwidth (Hz) available to network 


Ass.[3Tn 


B 


number of bands employed, each of BW W/ B (Hz) 


Ass.O 


V 


noise power spectral density (W/Hz) 


Def.O 


N{B),N 


noise power over a band, full spectrum 


Def.O 


V 


spectral efficiency requirement 


Def.O 


qiX/B,B) 


OP for intensity A and B bands 


Def.lOl 


\{g*,B*) 


TC for target OP q* using optimal bands B 


Def.O 


K{q*),LoiB*) 


spatial component and spectral component of the TC 


Rem.lOl 


Ik 
V 


energy per bit 


Def.O 


K 


interference cancellation effectiveness 


Def.O 


K 


maximum number of cancellable interferers 


Def.O 


p ■ 

^ min 


minimum receive power required for cancellation 


Def.O 


zp^(o),np^,(o) 


(random) partially cancellable (pc) int. /interferers at o 


Def.O 


z-(o),n-;,(o) 


(random) uncancellable (uc) int. /interferers at o 


Def.O 


' dom ' dom 


mapped constraint thresholds 


Lem.O 


ho,o 


(random) fading channel coef. from reference Tx to o 


Def.O 




(random) fading channel coef. from interferer i to its Rx 


Def.O 


hi,o 


(random) fading channel coef. from interferer i to o 


Def.O 


h 


fading threshold for transmission 


Def.O 


A 


intensity of attempted transmissions under threshold h 


Def.O 


A . 


(random) MPPP of attempted Tx's under threshold h 


Def.O 


q(h),A(h) 


OP and TP under threshold h 


Def.O 


f 


FPC exponent 


Def.O 




(random) Tx power by reference Tx and interferer i 


Def.O 


^0,/, 9/(0), 5/(0) 


event/probability of outage w/o interference 


Rem.O 


4^,A>t^;;;(o) 


(random) dominant interferers and interference 


Def.O 


*d,A,t^:'(o) 


(random) non-dominant interferers and interference 


Def.O 
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Table 3 Notation (continued) 



Chapter 6 


Vi 


(random) Tx beamforming vector for Tx i 


Def.O 


Wi 


(random) Rx beamforming vector for Rx i 


Dcf.lCT 


riT 


number of Tx antennas 


Def.lCT 


riR 


number of receive antennas 


Def.O 


H,h 


(random) matrix, vector channel coefficients 


^6.3.11 


So 


(random) signal power, So = ||ho|P 


Rem. EH 


K 


number of independent streams of data 


Dcf.lOl 


z 


number of interferers cancelled in PZF Rx, z < jir . 


Def.[631 


e G [0, 1] 


fraction of antenna array used for a given purpose 


Prop. ESI 




min. and max. distances to a desired Rx in a SDMA cluster 


mx2\ 
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Table 4 Acronyms 



ASE 


Area spectral efficiency 


BF 


Beamforming 


BLAST 


Bell Labs space time (a MIMO Rx) 


BPP 


Binomial point process 


iid 


Independent and identically distributed 


CCDF 


Complementary cumulative distribution function 


CDF 


Cumulative distribution function 


CF 


Characteristic function 


CSMA 


Carrier sense multiple access 


DoF 


Degrees of freedom 


DPC 


Dirty paper coding 


FLD 


Fixed link distances 


FPC 


Fractional power control 


FTS 


Fading threshold scheduling 


HRF 


Hazard rate function 


IC 


Interference cancellation 


LB 


Lower bound 


LT 


Laplace transfrom 


MAC 


Medium access control 


MGF 


Moment generating function 


MIMO 


Multiple-input multiple-output (multiple antennas) 


MISO 


Multiple-input single-output 


MMSE 


Minimum mean square error 


MPPP 


Marked Poisson point process 


MRC 


Maximal ratio combining 


MRT 


Maximal ratio transmission 


OCD 


Optimal contention density 


OP 


Outage probability 


PC 


Power control 


PDF 


Probability density function 


PGFL 


Probability generating functional 


PZF 


Partial zero forcing 


PPP 


Poisson point process 


RV 


Random variable 


Rx 


Receiver 


SDMA 


Space division multiple access 


SINR 


Signal to interference plus noise ratio 


SIR 


Signal to interference ratio 


SM 


Spatial multiplexing 


SN 


Shot noise 


SNR 


Signal to noise ratio 


TC 


Transmission capacity 


TP 


Throughput 


Tx 


Transmitter 


UB 


Upper bound 


VLD 


Variable link distances 


ZF 


Zero forcing 



A 
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Table A.l Results in Ch. [T] Introduction and preliminaries. 



§ 


o 


Motivation and assumptions 


Ass. 


o 


Key assumptions 


Fig. 


o 


Reference Rx and Tx and PPP of intorferers 






Key definitions: PPP, OP, and TC 


Def. 


o 


Outage probability 


Dcf. 


o 


Homogeneous Poisson point process (PPP) 


Fig. 


o 


Instance of a PPP 


Fact 


o 


OP g(A) is continuous, strictly increasing, onto (0, 1) 


Def. 


roi 


Transmission capacity (TC) 


Rem. 


o 


Potential vs. actual transmitters and Aloha 


§ 




Overview of the results 
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Table A. 2 Results in Ch.[2] Mathematical preliminaries. 



Def. 




JJdll dillj. OiliilUlLlO 




EI] 


Rail anrl anmilnc; "vrr^lnmo 


Def. 




G3ymm.ci fiiiictioii 


s 


lo 1 1 


Probability: notations, definitions, key inequalities 


Rem. 


lo 1 1 


RV notation 


JJer. 




Standard probability definitions 


Prop. 


lo ol 


Jensen's inequality 


Prop. 


jo q| 


IVlarkov's inequality 


Prop. 


lo /1 1 


Chebychev's inequality 


Prop. 


lo f^l 


Chernoff' s inequality 




12.21 


PPP void probabilities and distance mappings 


Ass. 


12.11 


Labeling convention for PPP 


Prop. 


jo fij 
1.^.01 


Void probability 


Thm. 


lo 1 1 

12.1! 


Mapping theorem 


Prop. 


12.71 


Distance mapping 




12.31 


Shot noise (SN) processes 


Fig. 


12.11 


SN process 


Def. 


12.41 


SN process 


Rem. 


12.21 


SN index convention 


Rem. 


12.31 


Radial symmetry 


Ass. 


12.21 


Power law impulse response 


Def. 


12.51 


Power law bN and characteristic exponent 


Rem. 


12.41 


Pathloss attenuation and the singularity at the origin 


Def. 


12.61 


Frechet distribution 


Cor. 


12.11 


Max SN RV CDF 




lo q| 


Interference mapping 


Trim. 


lo ol 
12. Zl 


Integration of radially symmetric functions 


Thm. 


lo qI 
12.31 


Campbell-JVIecke 


Prop. 


lo qI 


Shot noise mean and variance 


Prop. 


IZ.IUI 


Shot noise series expansion 


oor. 


|o o| 


jT.syrnpLOLic i ur ano o u r oi luc oin ixv 




I2.4I 


Stable distributions, Laplace transforms, and PGFL 


Def. 


12.71 


btable KV and distribution 


Def. 


lo «| 


Stable CF 


Def. 


[H 


Levy distribution 


Prop. 


12.111 


Stable moments 


Fig. 


m\ 


Levy PDFs and CDFs 


Def. 


I2.101 


Point process PGFL 


Prop. 




PPP PGFL 


Cor. 


[231 


Point process SN LT 


Cor. 




PPP SN LT 


Cor. 


[231 


Pathloss SN MGF 


Cor. 


EH 


Pathloss SN CF 


Cor. 


Em 


Pathloss SN for (5 = i 


§ 


12.51 


Maximums and sums of RVs 


Prop. 


|2.13| 


Sum and max SN CCDF ratio 


Fig. 


[231 


CCDFs for sum and max SN 


Def. 


|2.11| 


Subexponential distribution 


Prop. 


|2.14| 


Sufficient subexponential condition) 


Def. 


|2.12| 


binomial point process (BPP) 


Lem. 


EH 


BPP distances and interference 


Cor. 


[231 


Subexponential BPP interferences 
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Table A. 3 Results in Ch. |3] Basic model. 



Def. 


\M\ 


Basic model SINR 


Def. 




Rjc SNR 


Ass. 




SNR LB 


§ 


m 


Exact OP and TC 


Prop. 


[311 


OP is SN CCDF 


Cor. 




Explicit OP for <5 = i 


Prop. 


[32] 


TC (e = 0) 


Cor. 


[Q 


TC (e = and 5 = |) 


Fig- 


[Q 


Exact OP and TC 




[3:2] 


Asymptotic OP and TC 


Prop. 


[Q 


Asymptotic OP and TC 


Rem. 


[311 


TC as sphere packing 


Fig. 




Taylor series expansions used in asymptotic TC 


§ 




Upper bound on TC and lower bound on OP 


Def. 


[33] 


Dominant interferers and interference 


Prop. 


[3:4] 


OP LB and TC UB 


Rem. 


[3:2] 


Dominant and maximum interferers 


Cor. 


[33] 


OP and TC bounds (e = 0, N = 0, S = ^) 


Fig. 


[33] 


OP and TC exact results vs. bounds 


Fig. 


EM 


Lower bound on normal CDF 




m 


Throughput (TP) and TC 


Def. 


[Ml 


MAC layer TP 


Prop. 


[Ml 


Slotted Aloha TP and OP 


Fig. 


[331 


Slotted Aloha TP and OP 


Prop. 


EU 


MAC layer TP UB 


Prop. 


[TTl 


TC is constrained TP maximization 


Prop. 


EU 


Maximum TP equals maximum TC 


Cor. 


[3:41 


Maximizing TP and TC UBs 


Fig. 


EE 


Throughput and TC UBs 




[331 


Lower bounds on TC and upper bounds on OP 


Prop. 


[331 


Exact OP in terms of OP LB 


Prop. 


|3.10| 


Markov inequality OP UB 


Prop. 


[31T1 


Chebychev inequality OP UB 


Prop. 


|3.12| 


Chernoff inequality OP UB 


Fig. 




The three OP upper and TC LBs 
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Tabic A. 4 Results in Ch. |4] Extensions of the basic model 







Clianiiel fading 


Def. 




SINR. undGr fading 


Ftem. 


ITTl 


Signal and interference fading coefficients 




k.i.ii 


Exact OP and TC with fading 


Prop. 




LT of the interference 






OP and TC under Ftayleigh signal fading 


Leni. 


ITTl 


Moments of exponential RV 


Cor. 




OP and TC under Ftayleigh fading 


Cor. 


[4:21 


OP and TC under Rayleigh fading {6 = ^, N = 0) 


Fig. 




OP and TC under Rayleigh fading 


SS 


14 1 21 


Asymptotic OP and TC with fading 






Marked PPP fMPPP^ 


Thm. 




PPP marking theorem 


Prop. 




Void probability of the non-homogeneous M!PPP 




l4~4l 


A^IPPP distance and interference mapping 


Prop 


[431 


Shot noise series expansion 


Thm. 


[4:21 


Interference with fading is stable 


Rem. 


1131 


Interference and fading moments 


Rem. 


14:41 


Fading and outage with no interference 


Prop. 


11:61 


Asymptotic OP and TC under fading 


Rem. 


1131 


OP and fading moments 


Cor. 


1131 


Fading degrades performance 


Fig. 


11:21 


Fading degrades performance 




|4.1.3| 


Lower (upper) bound on OP (TC) with fading 


Def. 


[U 


Dominant intcrfcrers and interference 


Prop. 


[Ml 


OP LB 


Fig. 


m 


MGF for the RV -h^* 


Fig. 


11:11 


Exact, asymptotic, bound OP and TC under fading 




14.21 


Variable link distances (VLD) 


Def. 


m 


SINR for VLD 


Def. 


11:11 


OP for VLD 


Prop. 


HI 


Asymptotic OP and TC (e = 0, Af = 0) 


Prop. 


11:91 


OP LB as an MGF 


Prop. 


|4.10| 


Nearest neighbor RV characteristics 


Cor. 


HI 


Exact OP (e = 0, AT = 0, (5 = i) 


Fig. 


1131 


Exact, LB, and asymptotic OP 


§ 


1131 


Multihop TC 


Fig. 


14.61 


The multihop TC model with M = 3 


Def. 


1131 


Multihop TC 


Prop. 




Multihop TC inequality 


Prop. 


|4.12| 


Multihop TC UB 


Fig. 


Kt\ 


Multihop TC and its UB vs. allowed number of Tx attempts A 


Def. 


1131 


Optimal number of hops 


Prop. 




Optimal number of hops 


Cor. 


1131 


Optimal number of hops for a = 3 


Cor. 


1131 


Optimal number of hops for a = A 


Fig. 


[131 


Amh(A, U, M, A) vs. M for A e {6, 12} 
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Table A. 5 Results in Ch. [5] Design techniques for wireless networks 





|M A 1 

15.11 


Spectrum management 


Ass. 


15.11 


Random band selection 


Rem. 


IP HI 

15.11 


rim ■ 1 ■ J r i r t~i 

Inmned mterierence seen by reference Kx 


Def. 


15.11 


Noise, SINR, SNR, capacity, spectral efficiency 


Def. 


15.21 


OP and 1 C under multiple bands 


Prop. 


1 r -1 1 
15.11 


IC under multiple bands 


Rem. 


15.21 


Optimal number of bands independent of target OP 


Def. 


15.31 


Energy per bit anad receive SNR 


Rem. 


15.31 


Relaxation of integrality constraint 


Cor. 


15711 


TC under multiple bands 


Prop. 


15.21 


Minimum energy per bit required for solution 


Rem. 


15.41 


Low SNR regime 


Fig. 


1 r- 1 1 

15.11 


Spectral component of performance 


Prop. 


lOol 


Optimal spectral efficiency 


t Ig. 


10. 


Optimal spectral efficiency 






Asymptotic optimal spectral efficiency (high SNR) 


Fig. 


1 r: ol 


Optimal high SNR spectral efficiency 


Prop. 


Ik /1 1 
|0-4| 


Maximum possible spectral efficiency (low SNR) 


t ig. 


Ic /I 1 
10.41 


Optimal spectral efficiency 


§ 


15.21 


Interference cancellation (iC) 


Def. 


15.41 


Ine (k, -K, i^niin) 1^ model 


Def. 


Ic cl 

15.51 


bllNK at a (z^, K , i^jnin) 1^ capable reterence Kx 


Lem. 


15.11 


Constraint mapping 


Lem. 


15.21 


T~»j."n nil / nil j 

Partially canceiiabie / uncancellable nodes 


^ Ig. 


Ir: cl 
10.01 


Illustration of the thresholds 


1 iim. 


Ic 1 1 
10.11 


Ordered distances marginal distributions 


Oor. 


10. ol 


Ordered distances in Hi i marginal distributions 


Fig. 


Ic c| 
10. Ol 


PDFs for the ordered distances in TTi.i 


Prop. 


1 r: c 1 

10. ol 


L;i3 under A , i^T^i^j IC K,x model 


l^lg. 


15.71 


OP vs. A 


l^lg. 


Ic q| 
10. Ol 


OP vs. K and -Pmin 


Ftem. 


10.01 


Winer oounos on wj: ano i unoer 


§ 


\rr ol 

15.31 


lading threshold scheduling (1 lb) 


Def. 


[M] 


Fading coefficients, signal interference, and SINR 


Def. 


15.71 


OP and TP 


Rem. 


15.61 


Quantity vs. quality of transmissions through lib 




|r: c| 
10. Ol 


A cnrrv^ 1 tit' nn/^/ZiV Tv I Q 

/\aynipLOLie Kjr mioei jt ± o 


Prop. 


EH 


Asymptotic TP under FTS 


Fig. 


[5:91 


Asymptotic TP and optimal FTS threshold 


Prop. 


15:81 


FTS exploits fading to improve performance 


Fig. 


|5.10| 


Throughput under FTS 


Prop. 


EH 


OP LB and TP UB 


Fig. 


15.111 


Asymptotic TP and UB 






Fractional power control (FPC) 


Def. 


[5:81 


Transmission powers and SINR 


Lem. 


01 


Transmitted power RV moments under FPC 


Fig. 


|5.12| 


Variance of the transmitted power under FPC 


Rem. 


[5:71 


Fading and outage with no interference 


Prop. 




Asymptotic OP and TC under FPC 


Prop. 


15.111 


Asymptotic optimality of / = 1/2 


Rem. 


EH 


Optimal FPC may incur large power variance 


Def. 


[5:91 


Dominant and maximum interferers 


Prop. 




The OP LB under FPC 


Fig. 


|5.13| 


OP vs. FPC exponent / 


Fig. 




OP vs. intensity of attempted transmissions A 



163 



Table A. 6 Results in Ch. |6] Multiple antennas 



§ 




MIMO with interference 


§ 


16:21 


Categorizing MIMO in decentralized networks 


Fig. 




Receive diversity, SM, and SDMA 


§§ 


|6.2.1| 


Single stream techniques 


Def. 


Em 


MIMO single stream SINR 


§§ 


|6.2.2| 


Multi-stream models: spatial multiplexing and SDMA 


Def. 


El 


Multi-stream OP and OCD 


Def. 




Multi-stream TC 




le.sl 


Single stream MIMO TC results 


§§ 


16.3.11 


Diversity 


Def. 




Single stream MIMO optimal linear diversity filters 


Thm. 


leH 


Maximal Ratio Combiner (MRC) 


Rem. 


lem 


Interference distribution is unchanged 


Rem. 


16:21 


Signal distribution 


Prop. 


lem 


OP with MRC 


Prop. 


16:21 


TC with MRC 


Rem. 


lOl 


Sublinear gain of MRC in riR 


Prop. 


lOl 


OCD of riT X riR eigenbeamforming 


Rem. 


lOl 


TC scaling with Ut x Ur eigenbeamforming 


§§ 


|6.3.2| 


Interference cancellation (IC) 


Def. 


l631 


Partial zero forcing (PZF) Rx 


Prop. 


El 


PZF SINR 


Prop. 


EH 


OP for PZF-2 


Prop. 


EH 


TC for PZF LB 


Rem. 


EH 


PZF linear scaling in riR 


Def. 


EH 


MMSE Rx filter 


Rem. 


EH 


MMSE Rx filter 


Prop. 


E3 


TC UB for PZF and MMSE 


Cor. 


ED 


TC UB for MRC 


Fig. 


EH 


OCD vs. riu for PZF, MMSE, ZF, and MRC, for a = 4 


§ 


El 


Main results on multiple stream TC 


§§ 


|6.4.1| 


Spatial multiplexing 


Prop. 


Esl 


SM with MRC 


Prop. 


EH 


SM with ZF Rx 


Prop. 


|6.10| 


SM with PZF at high SNR 


Cor. 


El 


PZF Rx without cancellation 


Cor. 


lOl 


PZF Rx cancelling strongest interferer 


Rem. 


E3 


SM with the MMSE Rx 


Prop. 


|6.11| 


SM with multimode beamforming and PZF receivers 


Prop. 


ET2| 


SM with the BLAST architecture 


§§ 


|6.4.2| 


Space division multiple access (SDMA) 


Prop. 


ET3| 


Bounds on multistream DPC TC with MRC 


Prop. 


|6.14| 


SDMA scaling laws with a DPC Tx and MRC Rx 


Cor. 


El 


TC scaling and K* for SDMA with DPC Tx and MRC Rx 


Prop. 


16.151 


TC with DPC Tx and single antenna Rx 




16.51 


Practical issues and further research 


§^§ 


Ifi.B.ll 


Summary of main design insight 


§§ 


16.5.21 


Caveats and practical issue 



Acknowledgements 



The authors wish to acknowledge several individuals with whom we 
have collaborated in the development of the transmission capacity 
framework over the last eight years. First, many of the ideas and re- 
sults in this volume would not exist without our active collaboration 
with Dr. Nihar Jindal (now at Broadcom). In particular, he was lead 
author on several journal papers ([l5l HTl HB]) on the topics discussed 
in Chapters 5 and 6. 

The authors wish to also acknowledge the seminal contributions 
made by Prof, de Veciana (UT Austin) and his former Ph.D. student 
Dr. Xiangying Yang (Apple) in two early papers on transmission capac- 
ity. The authors learned stochastic geometry largely via a course (and 
the accompanying course notes) taught by Dr. de Veciana in Spring 
2003 — this course led directly to the research collaboration on trans- 
mission capacity [501 [77] . Other collaborators and colleagues who pro- 
vided valuable input and/or feedback to this monograph are as follows, 
in alphabetical order: 

• Prof. Frangois Baccelh (ENS and INRIA) 

• Prof. Radha Krishna Ganti (IIT Madras) 



164 



Acknowledgements 165 



• Prof. Martin Haenggi (University of Notre Dame) 

• Prof. Robert W. Heath Jr. (The University of Texas at 
Austin) 

• Prof. Kaibin Huang (Yonsei University) 

• Mr. Andrew Hunter (The University of Texas at Austin) 

• Prof. Marios Kountouris (Supelec) 

• Dr. Chun-Hung Liu (The University of Texas at Austin) 

• Dr. Raymond Louie (Univ. of Sydney) 

• Prof. Rahul Vaze (TIFR Mumbai) 

The support of a 2006-2009 NSF grant to Weber, Andrews, and Jindal 
is gratefully acknowledged, as is NSF CAREER CCF-0643508 (An- 
drews). This vohimc would not have come about without the intel- 
lectual environment and financial support provided by the DARPA 
IT-MANET project (#W911NF-07-l-0028). 



References 



[1] N. Abramson, "The throughput of packet broadcasting channels," IEEE Trans- 
actions on Communications, vol. COM-25, pp. 117-128, January 1977. 

[2] O. Ali, C. Cardinal, and F. Gagnon, "Performance of optimum combining in a 
Poisson field of interferers and Raylcigh fading channels," IEEE Transactions 
on Wireless Communicaiions, vol. 9, pp. 2461-2467, August 2010. 

[3] J. G. Andrews, F. Baccelli, and R. K. Ganti, "A tractable approach to coverage 
and rate in cellular networks," IEEE Transactions on Communications, vol. 59, 
pp. 3122-3134, November 2011. 

[4] J. G. Andrews, W. Choi, and R. W. Heath, "Overcoming interference in spatial 
multiplxing MIMO cellular networks," IEEE Wireless Communications Maga- 
zine, vol. 14, pp. 95-104, December 2007. 

[5] J. G. Andrews, S. Weber, M. Kountouris, and M. Haenggi, "Random access 
transport capacity," IEEE Transactions on Wireless Communications, vol. 9, 
pp. 2101-2111, June 2010. 

[6] F. Baccelli, B. Blaszczyszyn, and P. Muhlcthalcr, "An ALOHA protocol for 
multihop mobile wireless networks," IEEE Transactions on Information The- 
ory, vol. 52, pp. 421-436, February 2006. 

[7] F. Baccelli, J. Li, T. Richardson, S. Shakkottai, S. Subramanian, and X. Wu, 
"On optimizing CSMA for wide area ad-hoc networks," in Modeling and Opti- 
mization in Mobile, Ad Hoc and Wireless Networks (WiOpt), (Princeton, NJ), 
pp. 354-359, May 2011. 

[8] F. Baccelli and B. Blaszczyszyn, "Stochastic geometry and wireless networks 
— volume I: Theory," Foundations and TYends@ in Networking, vol. 3, no. 3-4, 
pp. 249-449, 2009. 



166 



References 167 



[9] F. Baccelli and B. Blaszczyszyn, "Stochastic geometry and wireless networks 
— volume II: Applications," Foundations and Trends® in Networking, vol. 4, 
no. 1-2, pp. 1-312, 2009. 

[10] J. A. Baker, "Integration of radial functions," MAA Mathematics Magazine, 
vol. 72, pp. 392-395, December 1999. 

[11] R. Blum, "MIMO capacity with interference," IEEE Journal on Selected Areas 
in Communications, vol. 21, pp. 793-801, June 2003. 

[12] R. Blum, J. Winters, and N. Sollenberger, "On the capacity of cellular systems 
with MIMO," IEEE Communications Letters, vol. 6, pp. 242-244, June 2002. 

[13] G. Cairo and S. Shamai, "On the achievable throughput of a multi-antenna 
Gaussian broadcast channel," IEEE Transactions on Information Theory, 
vol. 49, pp. 1691-1706, July 2003. 

[14] S. Catreux, P. Driessen, and L. Greenstein, "Attainable throughput of an 
interference-limited multiple-input multiple-output (MIMO) cellular system," 
IEEE Transactions on Communications, vol. 49, pp. 1307-1311, August 2001. 

[15] B. Chen and M. J. Cans, "MIMO communications in ad hoc networks," IEEE 
Transactions on Signal Processing, vol. 54, pp. 2773-2783, July 2006. 

[16] Y. Chen and J. G. Andrews, "An upper bound on multi-hop transmission 
capacity with dynamic routing selection," in IEEE International Symposium 
on Information Theory (ISIT), (Austin, TX), pp. 1718-1722, June 2010. 

[17] V. P. Chistyakov, "A theorem on sums of independent positive random variables 
and its applications to branching random processes," Theory of Probability and 
Its Applications, vol. 9, pp. 640-648, 1964. 

[18] H. Cho and J. Andrews, "Resource-redistributive opportunistic scheduling for 
wireless systems," IEEE Transactions on Wireless Communications, vol. 8, 
pp. 3510-3522, July 2009. 

[19] T. L. Chow and J. L. Teugels, "The sum and the maximum of I.I.D. random 
variables," in Proceedings of the 2nd Prague Symposium on Asymptotic Statis- 
tics, pp. 81-92, August 1978. 

[20] M. Costa, "Writing on dirty paper," IEEE Transactions on Information The- 
ory, vol. 29, pp. 439-441, May 1983. 

[21] D. J. Daley and D. Vere-Jones, An introduction to the theory of point processes. 
Springer, 2nd od., 2003. 

[22] D. A. Darling, "The influence of the maximum term in the addition of indepen- 
dent random variables," Transactions of The American Mathematical Society, 
vol. 73, no. 1, pp. 95-107, 1952. 

[23] G. J. Foschini, "Layered space-time architecture for wireless communication 
in a fading environment when using multiple antennas," Bell Labs Technical 
Journal, vol. 1, no. 2, pp. 41-59, 1996. 

[24] A. E. Gamal and Y.-H. Kim, Network Information Theory. Cambridge Uni- 
versity Press, 2012. 

[25] R. K. Ganti, J. G. Andrews, and M. Haenggi, "High-SIR transmission capacity 
of wireless networks with general fading and node distribution," IEEE Trans- 
actions on Information Theory, vol. 57, pp. 3100-3116, May 2011. 



168 References 

[26] R. K. Ganti and M. Haenggi, "Interference and outage in clustered wireless 
ad hoc networks," IEEE transactions on Information Theory, vol. 55, no. 9, 
pp. 4067-4086, 2009. 

[27] H. Gao and P. Smith, "Some bounds on the distribution of certain quadratic 
forms in normal random variables," Australian and New Zealand Journal of 
Statistics, vol. 40, no. 1, pp. 73-81, 1998. 

[28] A. Ghosh, J. Zhang, J. G. Andrews, and R. Muhamed, Fundamentals of LTE. 
Prentice-Hall, 2010. 

[29] C. Goldie and C. Kliippelberg, "Subexponential distributions," in A practical 
guide to heavy tails: statistical techniques for analysing heavy tails, (R. Adler, 
R. Feldman, and M. Taqqu, eds.), Birkhauser, 1997. 

[30] A. J. Goldsmith, Wireless Communications. Cambridge University Press, 2005. 

[31] S. Govindasamy, D. W. Bliss, and D. H. Staelin, "Spectral efficiency in single- 
hop ad-hoc wireless networks with interference using adaptive antenna arrays," 
IEEE Journal on Selected Areas in Communications, vol. 25, pp. 1358-1369, 
September 2007. 

[32] S. Govindasamy, D. W. Bliss, and D. H. Staelin, "Asymptotic spectral efficiency 
of multi-antenna links in wireless networks with limited Tx CSI," Computing 
Research Repository (CORE), 2010. 

[33] J. A. Gubner, "Computation of shot-noise probability distributions and densi- 
ties," SIAM Journal of Scientific Computing, vol. 17, pp. 750-761, May 1996. 

[34] P. Gupta and P. R. Kumar, "The capacity of wireless networks," IEEE Trans- 
actions on Information Theory, vol. 46, pp. 388-404, March 2000. 

[35] M. Haenggi, "On distances in uniformly random networks," IEEE Transactions 
on Information Theory, vol. 51, pp. 3584-3586, October 2005. 

[36] M. Haenggi and R. K. Ganti, "Interference in large wireless networks," Foun- 
dations and Trends® in Networking, vol. 3, no. 2, pp. 127-248, 2008. 

[37] K. Huang, J. G. Andrews, R. W. Heath, D. Guo, and R. Berry, "Spatial inter- 
ference cancellation for multi-antenna mobile ad hoc networks," IEEE Transac- 
tions on Information Theory, to appear, available at arxiv.org/abs/0804.0813. 

[38] A. Hunter and J. G. Andrews, "Adaptive rate control over multiple spatial 
channels in ad hoc networks," in Workshop on Spatial Stochastic Models for 
Wireless Networks (SpaSWiN), (Bcrhn, Germany), April 2008. 

[39] A. M. Hunter, J. G. Andrews, and S. Weber, "Capacity scaling laws for ad 
hoc networks with spatial diversity," in Proceedings of the IEEE International 
Symposium on Information Theory (ISIT), (Nice, France), ,Iune 2007. 

[40] A. M. Hunter, J. G. Andrews, and S. Weber, "Capacity scaling of ad hoc net- 
works with spatial diversity," IEEE Transactions on Wireless Communications, 
vol. 7, pp. 5058-5071, December 2008. 

[41] A. M. Hunter, R. K. Ganti, and J. G. Andrews, "Transmission capacity of 
multi-antenna ad hoc networks with CSMA," in IEEE Asilomar Conference on 
Signals, Systems and Computers, (Pacific Grove, CA), pp. 1577-1581, Novem- 
ber 2010. 

[42] J. How and D. Hatzinakos, "Analytic alpha-stable noise modeling in a Pois- 
son field of interferers or scatterers," IEEE Transactions on Signal Processing, 
vol. 46, pp. 1601-1611, June 1998. 



References 169 



[43] H. Inaltekin, S. B. Wicker, M. Chiang, and H. V. Poor, "On unbounded path- 
loss models: effects of singularity on wireless network performance," IEEE 
Journal on Selected Areas in Communications (JSAC), vol. 27, pp. 1078-1092, 
September 2009. 

[44] N. .Jindal, "MIMO broadcast channels with finite rate feedback," IEEE Trans- 
actions on Information Theory, vol. 52, pp. 5045-5060, November 2006. 

[45] N. Jindal, J. G. Andrews, and S. Weber, "Bandwidth partitioning in decen- 
tralized wireless networks," IEEE Transactions on Wireless Communications, 
vol. 7, pp. 5408-5419, December 2008. 

[46] N. Jindal, J. G. Andrews, and S. Weber, "Multi-antenna communication in ad 
hoc networks: achieving MIMO gains with SIMO transmission," IEEE Trans- 
actions on Communications, vol. 59, pp. 529-540, February 2011. 

[47] N. Jindal, S. Weber, and J. G. Andrews, "Fractional power control for decen- 
tralized wireless networks," IEEE Transactions on Wireless Communications, 
vol. 7, pp. 5482-5492, December 2008. 

[48] Y. Kim, F. Baccelli, and G. de Veciana, "Spatial reuse and fairness in mo- 
bile ad-hoc networks with channel-aware CSMA protocols," submitted to IEEE 
Transactions on Information Theory, May 2011. 

[49] J. F. C. Kingman, Poisson Processes. Oxford Studies in Probability, Oxford 
University Press, 1993. 

[50] S. Kotz and S. Nadarajah, Extreme value distributions: theory and applications. 
World Scientific, 2001. 

[51] M. Kountouris and J. G. Andrews, "Transmission capacity scaling of SDMA 
in wireless ad hoc networks," in IEEE Information Theory Workshop (ITW), 
(Taormina, Italy), pp. 534-538, October 2009. 

[52] M. Kountouris and J. G. Andrews, "Capacity bounds on multiuser MIMO 
transmission in random wireless networks," submitted to IEEE Transactions 
on Information Theory, 2010. 

[53] M. Kountouris and J. G. Andrews, "Downlink SDMA with limited feedback 
in interference-limited wireless networks," submitted to IEEE Transactions on 
Wireless Communications, 2011. 

[54] P. Levy, "Proprictcs asymptotiques des sommes de variables aleatoires 
indcpendentes en enchaines," Journal de Mathematiques, vol. 14, pp. 347-402, 
1935. 

[55] R. Louie, M. McKay, N. Jindal, and I. CoUings, "Spatial multiplexing with 
MMSE receivers in ad hoc networks," in IEEE International Conference on 
Communications (ICC), (Kyoto, Japan), June 2011. 

[56] R. H. Y. Louie, M. R. McKay, and I. B. CoUings, "Open-loop spatial multi- 
plexing and diversity communications in ad hoc networks," IEEE Transactions 
on Information Theory, vol. 57, pp. 317-344, January 2011. 

[57] S. B. Lowen and M. C. Teich, "Power-law shot noise," IEEE Transactions on 
Information Theory, vol. 36, pp. 1302-1318, November 1990. 

[58] V. Mordachev and S. Loyka, "On node density — outage probability tradeoff 
in wireless networks," IEEE Journal on Selected Areas in Communications 
(JSAC), vol. 27, pp. 1120-1131, September 2009. 



170 References 



[59] J. P. Nolan, Stable Distributions - Models for Heavy Tailed Data. Birkhauser, 
2012. 

[60] A. Paulraj, D. Gore, and R. Nabar, Introduction to Space-time Wireless Com- 
munications. Cambridge University Press, 2003. 

[61] R. Ramanathan, "On the performance of ad hoc networks with bcamforming 
antennas," in ACM International Symposium on Mobile Ad Hoc Networking 
and Computing (MobiHoc), (Long Beach, CA), pp. 95-105, October 2001. 

[62] R. Ramanathan, J. Redi, C. Santivanez, D. Wiggins, and S. Polit, "Ad hoc net- 
working with directional antennas: a complete system solution," IEEE Journal 
on Selected Areas in Communications (JSAC), vol. 23, pp. 496-506, March 
2005. 

[63] S. O. Rice, "Mathematical analysis of random noise," Bell Systems Technical 
Journal, vol. 23, pp. 282-332, 1944. 

[64] G. Samorodnitsky and M. S. Taqqu, Stable Non-Gaussian Random Processes: 
Stochastic Models with Infinite Variance. Chapman and Hall, 1994. 

[65] W. Schottky, "Uber spontane stromschwankungen in verschiedenen elektriz- 
itatsleitern," Annalen der Physik, vol. 57, pp. 541-567, 1918. 

[66] A. Shah and A. M. Haimovich, "Performance analysis of maximal ratio com- 
bining and comparison with optimum combining for mobile radio communica- 
tions with co-channel interference," IEEE Transactions on Vehicular Technol- 
ogy, vol. 49, pp. 1454-1463, July 2000. 

[67] I. M. Slivnyak, "Some properties of stationary flows of homogeneous random 
events," Theory of Probability and its Applications, vol. 7, no. 3, pp. 336-341, 
1962. 

[68] K. Stamatiou, J. G. Proakis, and J. R. Zeidler, "Spatial multiplexing in random 
wireless networks," Advances in Electronics and Telecommunications, vol. 1, 
pp. 5-12, April 2010. 

[69] K. Stamatiou, F. Rossctto, M. Haenggi, T. Javidi, J. R. Zeidler, and M. Zorzi, 
"A delay-minimizing routing strategy for wireless multihop networks," in Work- 
shop on Spatial Stochastic Models for Wireless Networks (SpaSWiN), (Seoul, 
Korea), June 2009. 

[70] D. Stoyan, W. S. Kendall, and J. Mecke, Stochastic geometry and its applica- 
tions. John Wiley and Sons, 2nd ed., 1995. 

[71] E. Teletar, "Capacity of multi-antenna Gaussian channels," European Transac- 
tions on Telecommunications, vol. 6, pp. 585-595, November-December 1999. 

[72] D. Tse and P. Viswanath, Fundamentals of Wireless Communication. Cam- 
bridge University Press, 2005. 

[73] R. Vazc, "Throughput-dclay-rcliability tradeoff with ARQ in wireless ad hoc 
networks," IEEE Transactions on Wireless Communications, vol. 10, pp. 2142- 
2149, July 2011. 

[74] R. Vaze and R. W. H. Jr., "Transmission capacity of ad-hoc networks with mul- 
tiple antennas using transmit stream adaptation and interference cancelation," 
submitted to IEEE Transactions on Information Theory, December 2009. 

[75] S. Verdu, "Spectral efficiency in the wideband regime," IEEE Transactions on 
Information Theory, vol. 48, pp. 1319-1343, June 2002. 



References 171 



S. Vishwanath, N. Jindal, and A. Goldsmith, "Duality, achievable rates and sum 
rate capacity of the Gaussian MIMO broadcast channel," IEEE Transactions 
on Information Theory, vol. 49, pp. 2658-2668, October 2003. 
S. Weber, J. G. Andrews, X. Yang, and G. de Veciana, "Transmission capacity 
of wireless ad hoc networks with successive interference cancellation," IEEE 
Transactions on Information Theory, vol. 53, pp. 2799-2814, August 2007. 
S. Weber, J. G. Andrews, and N. Jindal, "The effect of fading, channel in- 
version, and threshold scheduling on ad hoc networks," IEEE Transactions on 
Information Theory, vol. 53, pp. 4127-4149, November 2007. 
S. Weber, J. G. Andrews, and N. Jindal, "An overview of the transmission 
capacity of wireless networks," IEEE Transactions on Communications, vol. 58, 
pp. 3593-3604, December 2010. 

S. Weber, X. Yang, J. G. Andrews, and G. de Veciana, "Transmission capacity 
of wireless ad hoc networks with outage constraints," IEEE Transactions on 
Information Theory, vol. 51, pp. 4091-4102, December 2005. 
H. Weingarten, Y. Steinberg, and S. Shamai, "The capacity region of the Gaus- 
sian multiple-input multiple-output broadcast channel," IEEE Transactions on 
Information Theory, vol. 52, pp. 3936-3964, September 2006. 
F. Xue and P. R. Kumar, "Scaling laws for ad hoc wireless networks: an infor- 
mation theoretic approach," Foundations and Trends® in Networking, vol. 1, 
no. 2, pp. 145-270, 2006. 

S. Ye and R. S. Blum, "On the rate regions for wireless MIMO ad hoc net- 
works," in IEEE Fall Vehicular Technology Conference (VTC), (Los Angeles, 
CA), pp. 1648-1652, September 2004. 

W. Yu and J. Cioffi, "Sum capacity of Gaussian vector broadcast cliauiiols," 
IEEE Transactions on Information Theory, vol. 50, pp. 1875-1892, September 
2004. 

L. Zheng and D. Tse, "Diversity and multiplexing: A fundamental tradeoff in 
multiple antenna channels," IEEE Transactions on Information Theory, vol. 49, 
pp. 1073 1096, May 2003. 

M. Zorzi, J. Zeidler, A. Anderson, B. Rao, J. Proakis, A. L. Swindlehurst, 
M. Jensen, and S. Krishnamurthy, "Cross-layer issues in MAC protocol design 
for MIMO ad hoc networks," IEEE Transactions on Wireless Communications, 
vol. 13, pp. 62-76, August 2006. 



